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Abstract 

This document will describe the current status on the search for 
factors of Smarandache consecutive numbers and their reverse. A 
complete list up to index 200 will be given, with all known factors. 
Smarandache numbers are the concatenation of the natural numbers 
from one up to the given index, and reverse Smarandache numbers are 
the concatenation of the natural numbers from the given index down 
to 1. 



1 Introduction 

As a followup to Ralf Stephan’s article in this journal [St], I decided to extend 
his factorizations to index 200. The Smarandache consecutive sequence, 
as well as their reverse is described in [Sm]. In this article Smll denotes 
1234567891011 for example, and Rsmll denotes 1110987654321. 

Most of the factors that have been found by me and others, have been 
found by using the elliptic curve method (ECM) [Le], some have been found 
using the Multiple Polynomial Quadratic Sieve (MPQS) [Si]. 

All factors and remaining cofactors have been proven prime or composite 
by means of Elliptic Curve Primality Proving (ECPP) [At], or the Adleman- 
Pomerance-Rumely test [Ad], which has been simplified in 1984 by Cohen 
and Lenstra [Co]. 



2 Software used 

The main factoring program used has been GMP-ECM by Paul Zimmermann 
[Zi, Le, Gr]. The first small factors were filtered out quickly by ECMX, a 
program of the UBASIC package [Ki, Le]. 
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The factors which were probably prime were then tested with Frangois 
Morain’s ECPP [Mo, At]. Some factors have been proven prime by use of 
APRT-CLE [Ad] from the UBASIC package [Ki], 

All these fine pieces of software are freely available from the internet. 
The appropriate adresses are enlisted in the references. 

3 Progress of calculations 

All numbers have been factored using GMP-ECM up to 20 digits. First 25 
runs with B 1 = 2000 were run, and if the factorization wasn’t complete, 90 
runs with B 1 = 11000 were run. 

Work is in progress to extend this to 25 digits. Some factors have already 
been tried to 25 digits (300 curves with Bl = 50000). For more detail on 
the progress check the following URL: 

http : //ww.sci .kun.nl/ sigma/Persoonli jk/michaf /ecm/ecmtries .html 
Currently the lowest not-completely factored numbers are Sm63 and Rsm59. 



4 Factorization results 

The lists presented here are an up to date representation of the factors 
known so far. When more factors are found they will be added to the list, 
which can be found on the internet at the following URL: 
http : / / www . sci . kun . nl/ sigma/Persoonli jk/michaf /ecm/ 

Most of the factors up to Sm80 and Rsm80 should be creditted to Ralf 
Stephan, (unless otherwise stated). All contributors, together with their 
email- adresses can be found in tables 1 and 3. 

A denotes an un complete factorization, pxx denotes a prime of xx dig- 
its and cxx denotes a composite number of xx digits. 

4.1 Smarandache Factors 



Contributors of Smarandache factors 



RB 


Robert Backstrom 


TC 


Tim Charron 


BD 


Bruce Dodson 


MF 


Micha Fleuren 



bobbOatinet.com.au I 
tcharronOinterlog.com j 
bad001ehigh.edu 
michaf0sci.kun.nl 
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AM 


Allan MacLeod 


MACL-MSOfflwpmail . paisley .ac.uk 


RS 


Ralf Stephan 


stephanQtmt . de 


EW 


Egon Willighagen 


egonvfcsci . kun . nl 


PZ 


Paul Zimmermann 


zimmerma01oria.fr 
(LORIA, Nancy, France) 



Table 1: Contributors of Smarandache factors 



n 


Factors of Sm(n) 


2 


2^-3 


3 


3-41 


4 


2-617 


5 


3-5-823 


6 


2 6 • 3 - 643 


7 


127 • 9721 


8 


2 - 3 2 - 47 - 14593 


9 


3 2 - 3607 - 3803 


10 


2 - 5 • 1234567891 


11 


3 - 7 • 13 ■ 67 • 107 • 630803 


12 


2 s - 3 • 2437 • 2110805449 


13 


113 • 125693 • 869211457 


14 


2-3 

pl8 : 205761315168520219 


15 


3-5 

pl9 : 8230452606740808761 


16 


2 2 

plO : 2507191691 
pl3 : 1231026625769 


17 


3 2 - 47 • 4993 

pl8 : 584538396786764503 


18 


2 - 3 2 • 97 • 88241 

pl8 : 801309546900123763 


19 


13-43-79-281-1193 
pl8 : 833929457045867563 


20 


2 s - 3 • 5 - 323339 • 3347983 
pl6 : 2375923237887317 


21 


3 • 17 - 37 - 43 - 103 - 131 - 140453 
pl8 : 802851238177109689 


22 


2 • 7 - 1427 • 3169 - 85829 




continued... 
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n 


Factors of Sm (n) 




23 


p22 : 2271991367799686681549 
3 - 41 • 769 




24 


p32 : 13052194181136110820214375991629 
2 2 • 3 - 7 




25 


pl8 : 978770977394515241 
pl9 : 1501601205715706321 
5 2 15461 




26 


pll : 31309647077 

p25 : 1020138683879280489689401 

2 3 4 - 21347 • 2345807 






pl2 : 982658598563 




27 


pl8 : 154870313069150249 
3 3 • 19 2 • 4547 ■ 68891 




28 


p32 : 40434918154163992944412000742833 
2 3 • 47 • 409 




29 


pl5 : 416603295903037 

p27 : 192699737522238137890605091 

3-859 




30 


p20 : 24526282862310130729 
p26 : 19532994432886141889218213 
2 3 ■ 5 - 13 - 49269439 




31 


pl8 : 370677592383442753 
p23 : 17333107067824345178861 
29 




32 


plO : 2597152967 

p42 : 163915283880121143989433769727058554332117 
2 2 • 3 - 7 




33 


p23 : 45068391478912519182079 

p30 : 326109637274901966196516045637 

3 23 - 269 - 7547 




34 


pl8 : 116620853190351161 
p31 : 7557237004029029700530634132859 
2 






p50 : 6172839455055606570758085909601061116212631364146515661667 


35 


3 2 5 • 139 - 151 - 64279903 




36 


plO : 4462548227 

p37 : 4556722495899317991381926119681186927 
2 4 • 3 2 • 103 • 211 






p56 








continued ... 
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n Factors of Sm(n) 

37“ 71 • 12379 • 4616929 
p 52 

38 2-3 

p23 : 86893956354189878775643 

p43 : 2367958875411463048104007458352976869124861 

39 3 - 67 • 311 • 1039 

p25 : 6216157781332031799688469 

p36 : 305788363093026251381516836994235539 

40 2 2 - 5 - 3169 • 60757 • 579779 
p 10 : 4362289433 

p20 : 79501124416220680469 
p26 : 15944694111943672435829023 
i 41 3 • 487 • 493127 • 32002651 

p56 

42 2 • 3 - 127 - 421 
pll : 22555732187 

p25 : 4562371492227327125110177 

p34 : 3739644646350764691998599898592229 

43 7-17-449 
p 72 

44 2 3 • 3 2 

p26 : 12797571009458074720816277 
p 52 

45 3 2 - 5 • 7 - 41 - 727 - 1291 
pl3 : 2634831682519 

pl8 : 379655178169650473 

p41 : 10181639342830457495311038751840866580037 

46 2 • 31 • 103 • 270408101 
pl8 : 374332796208406291 

p25 : 3890951821355123413169209 
p28 : 4908543378923330485082351119 

47 3 • 4813 ■ 679751 

p22 : 4626659581180187993501 
p 53 

48 2 2 - 3 179 • 1493 • 1894439 

p29 : 15771940624188426710323588657 

p46 : 1288413105003100659990273192963354903752853409 

49 23 • 109 • 3251653 

1 plO : 2191196713 

continued... 
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n 


Factors of Sm(n) 


50 


p23 : 53481597817014258108937 

p47 : 12923219128084505550382930974691083231834648599 
2 - 3 • 5 2 • 13 - 211 - 20479 


51 


pl8 : 160189818494829241 
p20 : 46218039785302111919 

p44 : 19789860528346995527543912534464764790909391 
3 


52 


p20 : 17708093685609923339 
p 73 
2 7 


53 


pl7 : 43090793230759613 
p 76 
3 3 • 7 3 


54 


pl8 : 127534541853151177 
p76 

2 • 3 6 • 79 • 389 - 3167 • 13309 


55 


pll : 69526661707 

p22 : 8786705495566261913717 

p51 

5 • 768643901 


56 


pl5 : 641559846437453 

p22 : 1187847380143694126117 

p 55 

2 2 -3 


57 


p25 : 4324751743617631024407823 (BD) 
p 77 

3 - 17 • 36769067 


58 


pl3 : 2205251248721 

p37 : 2128126623795388466914401931224151279 (RB) 
p47 : 14028351843196901173601082244449305344230057319 
2-13 


59 


p31 : 1448595612076564044790098185437 (BD) 

p75 

3 


60 


pl8 : 340038104073949513 

p36 : 324621819487091567830636828971096713 (RB) 
p 55 

2 3 • 3 • 5 - 97 - 157 




pl03 




continued ... 
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n 


Factors of Sm(n) 


61 


10386763 

pl4 : 35280457769357 
p 92 


62 


2 • 3 2 • 1709 ■ 329167 • 1830733 

p34 : 970395623292182 1 22640122334854 1 28 1 (TC) 

p64 


63* 


3 2 

pll : 17028095263 
cl05 


64 


2 2 • 7 • 17 - 19 - 197 - 522673 
pl9 : 1072389445090071307 
p29 : 20203723083803464811983788589 (PW) 

p60 


65* 


3 • 5 - 31 • 83719 
cll3 


66* 


2 ■ 3 • 7 • 20143 • 971077 
clll 


67 


397 

pl8 : 183783139772372071 
pl05 


68* 


2 4 • 3 • 23 - 764558869 
plO : 1811890921 
cl 05 


69 


3-13-23 

p22 : 8684576204660284317187 

p24 : 281259608597535749175083 

p32 : 15490495288652004091050327089107 (RB) 

p49 : 3637485176043309178386946614318767365372143115591 


70 


2-5-2411111 

p24 : 109315518091391293936799 

p41 : 11555516101313335177332236222295571524323 

p60 


71 


3 2 • 83 • 2281 

p31 : 7484379467407391660418419352839 (AM) 
p95 


72 


2 2 -3 2 -5119 

p27 : 596176870295201674946617769 (BD) 
pi 03 


73* 


37907 




continued ... 
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71 


Factors of Sm(n) 


74 


"cl32 

2 • 3 - 7 - 1788313 • 21565573 


75* 


p20 : 99014155049267797799 

p25 : 1634187291640507800518363 (PW) 

p31 : 1981231397449722872290863561307 

p49 : 2377534541508613492655260491688014802698908815817 

3 • 5 2 • 193283 


76 


cl 33 
2 3 


77 


p 18 : 828699354354766183 

p27 : 213643895352490047310058981 

p9 7 

3 


78* 


p24 : 383481022289718079599637 (PW) 

p24 : 874911832937988998935021 

P39 : 164811751226239402858361187055939797929 (RB) 

p58 

2 • 3 • 31 185897 


79* 


cl39 
73 • 137 


80 


p20 : 22683534613064519783 
p24 : 132316335833889742191773 
cl02 

2 2 • 3 3 • 5 - 101 • 10263751 


81 


p25 : 1295331340195453366408489 
pi 15 
3 3 • 509 


82* 


p30 : 152873624211113444108313548197 (AM) 
pll9 

2 ■ 29 ■ 4703 • 10091 


83* 


p35 : 12295349967251726424104854676730107 (AVI) 

clll 

3 - 53 • 503 


84 


pl8 : 177918442980303859 (MF) 
cl 34 
2 5 • 3 


85* 


pl57 
5 ■ 7 2 




cl58 




continued ... 
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n 


Factors of Sm(n) 


86* 


2 3 23 1056149 ~ 

cl55 


87* 


3 - 7 - 231330259 
plO : 4275444601 (MF) 
cl 45 


88* 


2 2 

pl4 : 12414068351873 (MF) 
cl53 


89* 


3 • 3 • 13 • 31 • 97 • 163060459 
pl8 : 789841356493369879 (MF) 
cl37 


90* 


2 - 3 - 3 - 5 • 1987 • 179827 • 2166457 
cl54 


91* 


37 • 607 

pl6 : 5713601747802353 (MF) 

p24 : 100397446615566314002487 (MF) 

cl30 


92* 

! 


2 3 - 3 - 75503 
cl68 


93* 


3 - 73 • 1051 

plO : 3298142203 (MF) 
cl62 


94* 


2 • 12871181 

pll : 98250285823 (MF) 
cl60 


95* 


3 - 5 • 7 • 401 
cl76 


96 


2 - 2 • 3 - 23 • 60331 
pi 75 


97 


13 

pi 83 


98* 


2 • 3 2 • 23 • 37 - 199 

pl6 : 1495444452918817 (MF) 

cl 65 


99* 


3 2 • 31601 

pl2 : 786576340181 (MF) 
cl71 


100* 


2 2 • 5 2 - 7 3 - 8171 • 1065829 
plO : 2824782749 (AM) 




continued . . . 
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n 


Factors of Sm(n) 


101* 


p20 : 20317177407273276661 (MF) 

cl49 

3-8377 


102 


p21 : 799917088062980754649 (AM) 
cl69 

2 3 • 19 • 89 • 3607 • 15887 • 32993 


103* 


plO : 2865523753 (MF) 
pi 72 

131 • 1231 


104* 


pl6 : 1713675826579469 (MF) 
cl80 

2® - 3 • 59 • 773 


105* 


p20 : 19601852982312892289 (AM) 
cl 77 

3-5-193 


106* 


pl3 : 6942508281251 (MF) 
cl 90 

2 11 -127 *827 


107 


c203 

3 3 


108* 


pl2 : 536288185369 (MF) 

pl99 

2 2 - 3 3 


109* 


pl8 : 128451681010379681 (AM) 
cl96 

7 - 1559 - 78176687 


110 


p20 : 73024355266099724939 (AM) 
cl87 

2 • 3 - 5 • 4517 


111 


p20 : 18443752916913621413 (AM) 
pl97 

3 - 293 - 431 - 230273 • 209071 - 241423723 




plO : 3182306131 (MF) 

pl2 : 171974155987 (MF) 

pl3 : 1532064083461 (MF) 

pl7 : 59183601887848987 (MF) 

pl9 : 8526805649394145853 (AM) 

p23 : 27151072184008709784271 (AM) 

pl09 




continued... 
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n 


Factors of Sm(n) 


112 


2 J • 16619 • 449797 • 894009023 
pl7 : 74225338554790133 (MF) 
p23 : 10021106769497255963093 (MF) 
pl69 


113* 


3 • 11 • 13 • 5653 • 1016453 • 16784357 

pl8 : 116507891014281007 (AM) 

p37 : 6844495453726387858061775603297883751 (AM) 

cl57 


114* 


2 • 3 • 7 • 178333 
c227 


115* 


5 - 17 • 19 • 41 • 36606 • 71518987 
pl8 : 283858194594979819 (AM) 
c202 


116* 


2 2 - 3 2 - 2239 
c235 


117* 


3 2 • 31883 

pl2 : 333699561211 (MF) 

p20 : 28437086452217952631 (MF) 

c206 


118* 


2-83 

pll : 33352084523 (MF) 

p20 : 20481677004050305811 (MF) 

c214 


119* 


3 • 59 • 101 - 139 • 2801 
c239 


120* 


2 4 • 3 • 5 • 13 • 16693063 
c241 


121* 


278240783 

c246 


122 


2 • 3 • 23 - 618029123 

pl4 : 31949422933783 (MF) 

p233 


123* 


3 • 7 • 37 • 413923 

plO : 1565875469 (MF) 

pl6 : 5500432543504219 (MF) 

c227 


124* 

| 

j 


2 2 739393 

pl6 : 1958521545734977 (MF) 
c242 




continued ... j 
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n 


Factors of Sm(n) 


125* 


3 2 • 5 3 • 4019 

pl3 : 7715697265127 (MF) 
c247 


126 


2 • 3 2 - 29 • 103 • 70271 

p20 : 11513388742821485203 (MF) 

p241 


127* 


53 • 269 • 4547 

p20 : 56560310643009044407 (AM) 
c245 


128* 


2 3 • 3 • 7 - 11 - 59 - 215329 

p22 : 8154316249498591416487 (MF) 

c243 


129* 


3 - 19 
c277 


130* 


2-5 

pl2 : 166817332889 (MF) 
c269 


131* 


3 - 19 • 83 • 1693 

pll : 23210501651 (MF) 

pl2 : 575587270441 (MF) 

c256 


132* 


2 2 • 3 ■ 79 

pl3 : 2312656324607 (MF) 
c272 


133* 


pl9 : 8223519074965787731 (AM) 
c272 


134* 


2 - 3 3 • 73 - 6173 

pl6 : 5527048386371021 (AM) 

p28 : 1417349652747970442615118133 (AM) 

c243 


135* 


3 3 - 5 • 11 • 37 • 647 
plO : 1480867981 (MF) 
pl2 : 174496625453 (MF) 
pl5 : 151994480112757 (MF) 
c255 


136* 


2 5 • 1259 -4111 

pl3 : 9485286634381 (MF) 

p26 : 10151962417972135624157641 (AM) 

c253 




continued ... 
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n 


Factors of Sm(n) 


137* 


3 - 7 2 

pl3 : 7459866979837 (MF) 
c288 


138* 


2 • 3 • 181 • 78311 • 914569 
pl5 : 413202386279227 (MF) 
c 277 


139* 


13 

pll : 62814588973 (MF) 

pl2 : 115754581759 (MF) 

pl2 : 964458587927 (MF) 

p22 : 9196988352200440482601 (MF) 

c252 


140* 


2 2 • 3 - 5 • 23 • 761 - 1873 • 12841 

pll : 34690415939 (MF) 

pl8 : 226556543956403897 (AM) 

p23 : 10856300652094466205709 (AM) 

c248 


141 


3 • 107171 
p309 


142* 


2 • 7 - 4523 • 14303 • 76079 

p22 : 2244048237264532856611 (AM) 

c282 


143* 


3 2 - 859 
c317 


144 


2 s • 3 2 - 6361 

pl3 : 6585181700551 (MF) 

pl4 : 81557411089043 (MF) 

p21 : 165684233831183308123 (MF) 

p271 


145* 


5 • 96151639 
c326 


146* 


2-3-13-83 

pl2 : 720716898227 (MF) 

pl9 : 1122016187632880261 (MF) 

c296 


147* 


3-59113 

p22 : 1833894252004152212837 (AM) 

p31 : 1519080701040059055565669511153 (MF) 

c276 




continued... 
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n 


Factors of Sm(n) 


148* 


2 s • 197 • 11927 • 17377 • 273131 • 623321 
pl3 : 3417425341307 (AM) 
pl3 : 4614988413949 (MF) 
c288 


149* 


3 - 103 • 131 • 1399 
c331 


150* 


2 • 3 • 5 2 • 11 • 23 

pl6 : 2315007810082921 (MF) 

p26 : 92477662071402284092009799 (MF) 

c296 


151* 


7 • 53 - 1801 • 3323 
c335 


152* 


2 4 • 3 2 • 131 • 10613 

p20 : 29354379044409991753 (AM) 

p22 : 2587833772662908004979 (MF) 

c298 


153* 


3 2 ■ 29 - 7237 • 6987053 • 8237263 • 389365981 
c322 


154* 


2-17-19-43 

pi 8 : 444802312089588077 (MF) 
p21 : 855286987917657769927 (EW) 
c311 


155 


3 • 5 • 66500999 

p24 : 223237752082537677918401 (EW) 
p323 


156* 


2 2 - 3 - 7 - 3307 
c354 


157* 


11 • 53 - 492601 • 43169527 
pl2 : 645865664923 (MF) 
pl8 : 125176035875938771 (MF) 
c318 


158* 


2 - 3 ■ 17 - 29 ■ 53854663 

p21 : 164031369541076815133 (EW) 

c334 


159* 


3 - 71 • 647 

plO : 3175105177 (AM) 

p25 : 1957802969152764074566129 (EW) 

c330 


160* 


2 3 • 5 • 37 - 130547 • 859933 21274133 




continued... 
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n 


Factors of Sm (n) 


161 


p27 : 122800249349203273846720291 (EW) 
c324 

3 4 • 59 • 491 • 81705851 


162* 


p360 

2 - 3 5 • 2999 


163* 


p21 : 393803780657062026421 (AM) 

c351 

2381 


164* 


pll : 72549525869 (AM) 

pl2 : 666733067809 (AM) 

p25 : 1550529016982764630292633 (AM) 

c330 

2 2 -3 


165* 


c383 

3 - 5 ■ 7 - 13 • 31 - 247007767 


166 


pl5 : 490242053931613 (MF) 

c359 

2-89 


167* 


p23 : 55566524959746113370037 (AM) 
p365 
3 • 3313 


168 


c389 

2 7 - 3 - 532709 


169* 


p 387 

2671 • 5233 


170* 


c392 

2 • 3 2 • 5 • 7 • 701 


171* 


pl4 : 73406007054077 (MF) 

c382 

3 2 • 1237 


172* 


pl9 : 6017588157881558471 (AM) 
c382 

2 2 • 11 - 13 - 37 


173* 


c403 

3 - 17 • 53 - 101 - 153 • 11633 • 228673 


174* 


c394 

2 • 3 • 59 • 277 • 2522957 




p22 : 2928995151034569627547 (AM) 
c381 




continued ... 
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n 


Factors of Sm(n) 


175* 


5 2 

pl3 : 2606426254567 (MF) 
c403 


176* 


2 3 • 3 • 19 • 1051 

pl9 : 1031835687651103571 (AM) 
c396 


177* 


3 • 109 - 153277 • 6690569 
pll : 32545700623 (MF) 
pl6 : 2984807754776597 (MF) 
c382 


178 


2 

pl3 : 3144036216187 (MF) 
pi 7 : 11409535046513339 (MF) 
p397 


179* 


3 2 - 7 • 11 • 359 
c423 


180* 


2 2 - 3 2 • 5 - 43 • 89 - 7121 
c422 


181* 


31 • 197-70999 

p20 : 46096011552749697739 (AM) 
c406 


182* 


2 - 3 • 123529391 
c429 


183* 


3 - 29 • 661 • 1723 

pl6 : 3346484052265661 (AM) 

c417 


184* 


2 4 - 7 - 59 - 191 - 1093 - 1223 
pll : 22521973429 (MF) 
pl7 : 15219125459582087 (MF) 
pl8 : 158906425126963139 (MF) 
pl9 : 2513521443592870099 (MF) 
c369 


185* 


3 - 5 - 94050577 

pl3 : 4716042857821 (MF) 

pl6 : 3479131875325867 (MF) 

c409 


186* 


2 - 3 - 1201 

p21 : 574850252802945786301 (MF) 
c425 




continued... 
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n 


Factors of Sm(n) 


187* 


349 • 506442073 
c442 


188* 


2 2 -3 3 

c454 


189* 


3 3 - 47 - 1515169 

plO : 1550882611 (MF) 

plO : 1687056803 (MF) 

p21 : 348528133548561476953 (AM) 

c410 


190 


2-5-379 

p23 : 46645758388308293907739 (AM) 
p 435 


191* 


3 13 • 5233 

pl2 : 164130096629 (MF) 

p20 : 13806214882775315521 (MF) 

c429 


192* 


2 3 • 3 • 29 • 41 
c463 


193* 


7-419 

c467 


194* 


2-3-11-31-491-34188439 
pl4 : 28739332991401 (MF) 
pl6 : 8203347603076921 (MF) 
pi 9 : 1507421050431503839 (MF) 
p20 : 22805873052490568609 (MF) 
p21 : 168560953170124281211 (MF) 
c373 


195* 


3 • 5 - 397 - 21728563 • 300856949 - 554551531 

plO : 8174619091 (MF) 

c438 


196 


2 2 • 17 • 73 • 79 

plO : 3834513037 (MF) 

p465 


197* 


3 2 - 37 - 6277 

pl6 : 1368971104990459 (MF) 
c461 


198* 


2 ■ 3 2 ■ 7 2 • 13 
c482 


199* 


151 




continued... 
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n 


Factors of Sm(n) 


200* 


c487 
2 5 - 3 - 5 2 
c488 



Table 2: Factorizations of Sm(n), 1 < n < 200 
4.2 Reverse Smarandache Factors 



Contributors of Reverse Smarandache factors 


RB 

BD 

MF 

AM 

RS 

PZ 


Robert Backstrom 
Bruce Dodson 
Micha Fleuren 
Allan MacLeod 
Ralf Stephan 
Paul Zimmermann 


bobbQat inet . com . au 
badOQlehigh.edu 
michafQsci.kun.nl 
MACL-MSOQvpmail . paisley .ac.uk 
stephanQtmt . destephanQtmt . de 
Paul . ZimmermannQloria . f r 


continued... 





Table 3: Contributors of Reverse Smarandache factors 



n 


Factors of Rsm(n) 


2 


3.7 


3 


3.107 


4 


29.149 


5 


3.19.953 


6 


3.218107 


7 


19.402859 


8 


3 2 . 1997.4877 


9 


3 2 .17 2 .379721 


10 


7.28843.54421 


11 


3 




pl2 : 370329218107 


12 


3.7 




pl3 : 5767189888301 


13 


17.3243967.237927839 


14 


3.11.24769177 




plO : 1728836281 


15 


3.13.19 2 .79 


continued... 
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n 


Factors Rsm(n) 


16 


pl5 : 136133374970881 
23.233.2531 


17 


pl6 : 1190788477118549 

3 2 . 13. 17929.25411.47543.677181889 


18 


3 2 .11 2 . 19.23.281. 397.8577529.399048049 


19 


17.19 


20 


pl3 : 1462095938449 
pl4 : 40617114482123 
3.89.317.37889 


21 


p21 : 629639170774346584751 
3.37 


22 


pl2 : 732962679433 

pl9 : 2605975408790409767 

13.137.178489 


23 


pl3 : 1068857874509 
pl4 : 65372140114441 
3.7.191 


24 


p32 : 578960862423763687712072079528211 
3.107.457.57527 


25 


p28 : 28714434377387227047074286559 
11.31.59.158820811.410201377 


26 


p20 : 19258319708850480997 
3 3 . 929. 1753. 2503.4049. 11171 


27 


p24 : 527360168663641090261567 
3 5 .83 


28 


plO : 3216341629 
pl3 : 7350476679347 
pl8 : 571747168838911343 
23.193.3061 


29 


pl9 : 2150553615963932561 
p21 : 967536566438740710859 
3.11.709.105971.2901761 


30 


plO : 1004030749 

p24 : 405373772791370720522747 

3.73.79.18041.24019.32749 


31 


plO : 5882899163 

p24 : 209731482181889469325577 

7.30331061 


j 


p45 : 147434568678270777660714676905519165947320523 




continued... 
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n 



Factors Rsm(n) 

32 3.17.1231.28409 
pl2 : 103168496413 
p35 : 17560884933793586444909640307424273 

33 3.7.7349 

plO : 9087576403 

p42 : 237602044832357211422193379947758321446883 

34 89.488401 .2480227.63292783.254189857 
plO : 3397595519 

pl9 : 5826028611726606163 

35 3 2 .881.1559.755173.7558043 
plO : 1341824123 

pl6 : 4898857788363449 
pl6 : 7620732563980787 

36 3 2 .11 2 .971 

pl3 : 1114060688051 

p22 : 1110675649582997517457 

p24 : 277844768201513190628337 

37 29.2549993 

p20 : 39692035358805460481 

p38 : 12729390074866695790994160335919964253 

38 3.9833 
p63 

39 3.19.73.709.66877 
p 58 

40 11.41.199 

p27 : 537093776870934671843838337 

p39 : 837983319570695890931247363677891299117 

41 3.29.41.89.3506939 
pl4 : 18697991901857 

p20 : 59610008384758528597 

p28 : 3336615596121104783654504257 

42 3.13249.14159.25073 
plO : 6372186599 
p52 

43 52433 

p20 : 73638227044684393717 
p53 

44 3 2 .7.3067. 114883.245653 

p23 : 6571190708843766 0760939 

continued... 
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n 


Factors Rsm(n) 


45 


p41 : 12400566709419342558189822382901899879241 
3 2 .23.167.15859.25578743 


46 


p6 5 

23.35801 


47 


pl2 : 543124946137 

p23 : 45223810713458070167393 

p43 : 2296875006922250004364885782761014060363847 

3.11.31.59 


48 


pl6 : 1102254985918193 

p28 : 4808421217563961987019820401 

p38 : 14837375734178761287247720129329493021 

3.151.457.990013 


49 


pl5 : 246201595862687 

p24 : 636339569791857481119613 

p39 : 15096613901856713607801144951616772467 

71 


50 


plO : 9777943361 
p 77 

3.157.3307 


51 


pl3 : 3267926640703 

p30 : 771765128032466758284258631297 

p43 : 1285388803256371775298530192200584446319323 

3.11 


52 


p 92 

7.29.670001 


53 


pl2 : 403520574901 
pl4 : 70216544961751 
pi 6 : 1033003489172581 

p47 : 13191839603253798296021585972083396625125257997 
3 4 .499.673.6287.57653. 199236731 


54 


pl6 : 1200017544380023 
p28 : 1101541941540576883505692003 
p31 : 2061265130010645250941617446327 
3 3 .7 4 . 13.1427.632778317 


55 


pll : 57307460723 
pl3 : 7103977527461 
pl5 : 617151073326209 

p43 : 2852320009960390860973654975784742937560247 
357274517.460033621 




continued ... 
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n 



Factors Rsm(n) 
p 84 

56 3.13 2 

pl4 : 85221254605693 
p 87 

57 3.41 

pll : 25251380689 
p93 

58 11.2425477 

pl5 : 178510299010259 
pl8 : 377938364291219561 
p28 : 5465728965823437480371566249 
p40 : 5953809889369952598561290100301076499293 
59* 3 

cl09 

60 3 

plO : 8522287597 
plOl 

61 13.373 

p22 : 6399032721246153065183 

p42 : 214955646066967157613788969151925052620751 (RB) 
p46 : 9236498149999681623847165427334133265556780913 

62 3 2 .1 1.487.6870011 
pl3 : 3921939670009 
pl4 : 11729917979119 

p28 : 9383645385096969812494171823 

P50 : 43792191037915584824808714186111429193335785529359 

63 3 2 .97.26347 

p24 : 338856918508353449187667 

p86 

64 397.653 

pl2 : 459162927787 

pl4 : 27937903937681 

p24 : 386877715040952336040363 

p65 

65* 3.7.23.13219.24371 

cllO 

I 66 3.53.83.2857.1154129.9123787 

pl03 

i 67* 43 

L_ continued... 
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n 


Factors Rsm(n) 




68 


pll : 38505359279 
cl 13 

3.29.277213.68019179.152806439 




69 


pl8 : 295650514394629363 
p20 : 14246700953701310411 
p6 7 

3.11.71.167.1481 




70 


plO : 2326583863 

p23 : 19962002424322006111361 

p 89 

1157237.41847137 




71 


£>22 : 8904924382857569546497 
p96 

3 2 . 17.131. 16871 




72 


plO : 1504047269 

pll : 82122861127 

pl9 : 1187275015543580261 

p 87 

3 2 .449.1279 




73 


£>129 

7.11.21352291 




74 


plO : 1051174717 

pl7 : 92584510595404843 

p20 : 33601392386546341921 

£>83 

3.177337 






plO : 6647068667 
pll : 31386093419 
pl5 : 669035576309897 
pl6 : 4313244765554839 

p32 : 67415094145569534144512937880453 (PW) 


75 


p51 

3.7.230849.7341571.24260351 






plO : 1618133873 

pl4 : 19753258488427 

pl7 : 46752975870227777 

p28 : 7784620088430169828319398031 (PW) 


76* 


£>53 

53 








continued... 
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n 


Factors Hsm(n) 


77 


cl42 

3.919 


78* 


pl5 : 571664356244249 

p22 : 6547011663195178496329 (PW) 

p27 : 591901089382359628031506373 (BD) 

p33 : 335808390273971395786635145251293 (PW) 

p46 : 3791725400705852972336277620397793613760330637 

3.17.47 


79 


pl4 : 17795025122047 

cl31 

160591 


80* 


pl5 : 274591434968167 
pl9 : 1050894390053076193 

pi 12 

3 3 . 11.443291. 1575307 


81 


pl7 : 19851071220406859 
cl21 

3 3 .23 2 . 62273. 22193.352409 


82 


pl5 : 914359181934271 (MF) 
pl20 

PRIME'. (RS) 


83* 


3 


84* 


cl57 

3.11.47.83.447841.18360053 


85 


pl4 : 53294058577163 (MF) 
cl30 

pl2 : 465619934881 (MF) 


86* 


p22 : 5013354844603778080337 (AM) 

p 128 

3.7.3761 .20511 1.16080557. 16505767 


87 


cl39 

3.2423 


00 

00 

* 


p25 : 4433139632126658657934801 (AM) 
p30 : 951802198132419645688492825211 (MF) 
pl07 
73.8747 


89* 


cl 62 

3 2 . 19.7052207 




cl61 




continued... j 
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n 


Factors Rsm(n) 


90* 


3^.157.257.691 

cl40 


91* 


1 1.29. 163.3559.2297.22899893 

pl5 : 350542343218231 (MF) 

p25 : 8365221234379371317434883 (MF) 

cll5 


92* 


3.17.113.376589.3269443.6872137 

cl53 


93* 


3.13.69317.14992267 
cl 64 


94* 


7.593.18307 

pll : 51079607083 (MF) 
cl61 


95* 


3.11.13.53.157.623541439 

cl66 


96* 


3.7.211.14563.2297 
cl 72 


97* 


1553 

cl82 


98 


3 2 . 101.401.5741.375373 
pi 73 


99* 


3 2 . 109.41829209 

pl2 : 174489586693 (MF) 

cl 68 


100* 


13.6779 

pll : 48856332919 (MF) 

p26 : 41858129936073024200781901 (MF) 

cl50 


101* 


3 

pll : 16320902651 (MF) 

pl9 : 3845388775716560041 (MF) 

p33 : 693173763848292948494434792706137 (AM) 

cl32 


102* 


3.101.103.36749 

pll : 10189033219 (MF) 

p20 : 23663501701518727831 (AM) 

p26 : 52648894306108287380398039 (AM) 

cl 33 


103* 


19.29.103.3119.154009291 




continued... 
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n 


Factors Rsm(n) 


90* 


3^.157.257.691 

cl40 


91* 


1 1.29. 163.3559.2297.22899893 

pl5 : 350542343218231 (MF) 

p25 : 8365221234379371317434883 (MF) 

cll5 


92* 


3.17.113.376589.3269443.6872137 

cl53 


93* 


3.13.69317.14992267 
cl 64 


94* 


7.593.18307 

pll : 51079607083 (MF) 
cl61 


95* 


3.11.13.53.157.623541439 

cl66 


96* 


3.7.211.14563.2297 
cl 72 


97* 


1553 

cl82 


98 


3 2 . 101.401.5741.375373 
pi 73 


99* 


3 2 . 109.41829209 

pl2 : 174489586693 (MF) 

cl 68 


100* 


13.6779 

pll : 48856332919 (MF) 

p26 : 41858129936073024200781901 (MF) 

cl50 


101* 


3 

pll : 16320902651 (MF) 

pl9 : 3845388775716560041 (MF) 

p33 : 693173763848292948494434792706137 (AM) 

cl32 


102* 


3.101.103.36749 

pll : 10189033219 (MF) 

p20 : 23663501701518727831 (AM) 

p26 : 52648894306108287380398039 (AM) 

cl 33 


103* 


19.29.103.3119.154009291 




continued... 
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n 


Factors Rsm(n) 


104* 


pl2 : 329279243129 (MF) 
pl3 : 1240336674347 (MF) 
cl61 

3.7.60953.1890719 


105* 


pll : 10446899741 (MF) 
pl5 : 216816630080837 (MF) 
pl9 : 1614245774588631629 (MF) 
cl49 

3.7.859.6047.63601 


106* 


cl 94 

p22 : 1912037972972539041647 (AM) 


107* 


p22 : 3052818746214722908609 (AM) 
cl67 

3 3 . 13.45 19. 114967 


108 


plO : 1425213859 (MF) 
pl4 : 17641437858251 (MF) 
cl79 

3 3 . 23.457.1373 


109 


pl2 : 605434593221 (MF) 
pl2 : 703136513561 (MF) 
pl83 

11.29.31 2 . 1709.30345569 


110* 


pl4 : 42304411918757 (MF) 
pl89 

3.11.19.53.229.24672421 


111* 


p24 : 611592384837948878235019 (AM) 
cl83 

3.61.269.470077.143063.544035253 


112* 


c200 

137 


113* 


pl2 : 262756224547 (MF) 
c214 

3.19.45061.111211 


114* 


c219 

3.19.53.59 


115* 


c228 

137.509.1720003 


116 


c226 

3 2 . 83. 103. 156307. 176089.21769 127 




continued... 



30 







p 217 

117 3 2 
p242 

118 7.4603 
p241 

119* 3.7 

c247 

120* 3.73 

c249 

121* 31.371177 

c248 

122* 3.17 

pll : 91673873887 (MF) 
c245 

123 3.1197997 

pll : 15744706711 (MF) 
p244 

124* 37.1223 

c259 

125 3 2 . 59.83 

plO : 5961006911 (MF) 
pl3 : 1096598255677 (MF) 
p240 

126* 3 2 . 13.68879.135342173 

c255 

127* 97 

pl6 : 1385409249340483 (AM) 
c255 

128* 3.34613.29497667 

c263 

129* 3.23.1213.82507 

pl2 : 420130412231 (MF) 
c 257 

130* 31.263.86969.642520369 

c264 

131* 3.11.4111.852143 

pl2 : 606617222863 (MF) 

p23 : 33247682213571703426139 (AM) 

c239 

continued... 
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n 


Factors Rsm(n) 


132 


3.7.11.41.43.31259.69317.180307.199313 
pl7 : 16995472858509251 (MF) 
p20 : 56602777258539682957 (AM) 

p226 


133 


7.13 

p20 : 22533511116338912411 (AM) 
p269 


134* 


3 3 . 37.29004967 

pl7 : 60164048964096599 (AM) 
c266 


135* 


3 3 .211.5393.98563 

pl2 : 207481965329 (MF) 

p22 : 6789282931372049267693 (AM) 

c251 


136* 




137* 


3.179 

p22 : 6796599525965619205571 (AM) 
c278 


138* 


3.119611.314087617 

c292 


139* 




140* 


3.317.772477 

pl5 : 153629260660723 (AM) 
c289 


141* 


3.631.65831 

c307 


142* 


859.2377.2909.6521 

pl4 : 41190901651547 (MF) 

c291 


143 


3 2 . 93971 

pl2 : 9053448211979 (MF) 
p302 


144* 


3 2 

pi 9 : 5028055908018884749 (MF) 
c304 


145* 


57719.2691841 

p20 : 45690580335973653419 (MF) 
c296 


146* 


3.7 2 . 277. 19319.55807 




continued ... 
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n 


Factors Rsm(n) 


147* 


pl3 : 2454423915989 (MF) 
c304 

3.7 2 . 19.31. 15467623 


148* 


c321 

p20 : 33825333713396366003 (AM) 


149* 


p23 : 25082957895838310384953 (AM) 
c294 

3.109.34442413 


150* 


c329 

3.59.257 


151* 


c337 

plO : 7134941903 (MF) 


152 


c335 

3 2 .13 


153* 


p21 : 412891312089439668533 (MF) 

p325 

3 2 . 67793 


154* 


pl8 : 237333508084627139 (MF) 

c328 

11.53861 


155* 


plO : 1118399729 (MF) 
c339 

3.41.33842293 


156* 


c347 

3.21961 


157* 


c355 

plO : 4136915059 (MF) 


158* 


c353 

3.31.89209 


159* 


plO : 1379633699 (MF) 
pl4 : 54957888020501 (MF) 
c336 

3.13.5669.11213.816229087 


160* 


plO : 50611041883 (MF) 
c340 

7.942037.1223207 


161* 


p21 : 125729584994875519171 (AM) 
c339 

3 7 . 7.37.67.6521. 826811. 6018499 




continued ... 
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n 


Factors Rsm(n) 


162* 


p23 : 77558900444266075256801 (MF) 
c328 

3 4 . 1295 1 13.202557967 


163* 


c361 

pl6 : 1139924663537993 (MF) 


164 


pl7 : 17672171439068059 (MF) 

c350 

3.193 


165* 


p24 : 105444241520715055381519 (AM) 
p 358 
3 


166* 


c386 

pl5 : 396444477663149 (MF) 


167* 


p32 : 15221332593310506150048824812249 (AM) 
c344 

3.17.373.7346281.8927551.194571659 


168* 


p20 : 68277637362521294401 (AM) 
c347 

3.59.35537.68102449 


169* 

170 


pl9 : 7766035514845504007 (AM) 
c362 

3 2 .23. 


171* 


pl6 : 3737994294192383 (MF) 

p384 

3 2 .37 


172 


pl2 : 237089136881 (MF) 

pl9 : 2153684224509566597 (MF) 

p21 : 175530075465216996787 (MF) 

p22 : 8105319358780665120301 (MF) 

c330 

17.29.281 


173* 


plO : 4631571401 (MF) 
pll : 31981073881 (MF) 
pl5 : 119749047957053 (MF) 
p 368 
3.1787 


174* 


c407 

3.7.269.397.156894809 j 




continued... 
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n 


Factors Rsm(n) 


175* 


c399 

7.11 


176* 


plO : 3763462823 (MF) 
c405 

3.11.47.49613 


177 


pl3 : 2800890701267 (MF) 

pl5 : 315698062297249 (MF) 

p27 : 880613122533775176075766757 (MF) 

c358 

3.73.1753 


178 


pl4 : 29988562180903 (MF) 
p404 

13.47.353.644951.487703.1436731 


179* 


pl2 : 728961984851 (MF) 
pl4 : 34686545199997 (MF) 
pl4 : 36329334000803 (MF) 
p364 
3 2 .23.43 


180* 


pl4 : 50981967790529 (MF) 

c411 

3 2 .29 


181* 


pl7 : 33644294710009721 (MF) 
c413 

325251083 


182* 


pl7 : 57421731284347247 (MF) 
c410 

3.107.5568133 


183* 


pl2 : 139065644033 (MF) 

c417 

3.23.89 


184* 


c437 

23.19531 


185* 


pl5 : 196140464783429 (MF) 

c424 

3.13.919 


186* 


pll : 32173266383 (MF) 

c432 

3.23 




c448 




continued ... 
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Factors Rsm(n) 


187* 


61.83.103.523.3187 

pl9 : 1018598504636281577 (MF) 

c423 


188* 


3 3 .7.7681.65141 
c445 


189* 


3 3 .7.2039.3823.9739.212453.105865 19 
c433 


190* 


83.107.1871.25346653 

c447 


191 


3.809 

pl8 : 627089953107590081 (MF) 
p444 


192* 


3.2549 

c464 


193* 


47.503.12049 

c463 


194* 


3.179 

p22 : 8000103240831609636731 (AM) 
p23 : 77947886830169946060329 (MF) 
c426 


195* 


3.79.8219 

c471 


196* 


19 

pl6 : 8982588119304797 (AM) 
c463 


197* 


3 2 . 11. 43. 11743. 125201. 867619 
pll : 61951529111 (MF) 
pl4 : 27090970290157 (MF) 
c440 


198 


3 2 . 11.37.2837 

pl9 : 1245013373736039779 (MF) 
p461 


199* 


103.2377 

c484 


200* 


3.1666421 

c485 



Table 4: Factorizations of Rsm(n), 1 < n < 200 
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a load of factors to these lists, as well as some primality proofs. All other 
contributes earn a warm thank you here too. 
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Smarandache Continued Fractions 



Henry Ibstedt 1 



Abstract: The theory of general continued fractions is developed to the extent required in 
order to calculate Smarandache continued fractions to a given number of decimal places. 
Proof is given for the fact that Smarandache general continued fractions built with positive 
integer Smarandache sequences having only a finite number of terms equal to 1 is 
convergent. A few numerical results are given. 



Introduction 

The definitions of Smarandache continued fractions were given by Jose Castillo in the 
Smarandache Notions Journal, Vol. 9, No 1-2 [1]. 

A Smarandache Simple Continued Fraction is a fraction of the form: 



o(l)+- 



o(2)- 



o(3)- 



o(4) + 



1 



o(5) + - 



where a(n), for n>l, is a Smarandache type Sequence, Sub-Sequence or Function. 
Particular attention is given to the Smarandache General Continued Fraction defined as 



o(l)+- 



b{ 1) 



o(2) + - 



b{ 2) 



o(3) + - 



b( 3) 



o(4) + 



6(4) 
o(5) -i 



where a(n) and b(n), for n>l, are Smarandache type Sequences, Sub-Sequences or 
Functions. 



1 HTbstedtS; swipnet.se 
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As a particular case the following example is quoted 



1 + 



12 +- 



21 



123- 



321 



1234 + - 



4321 



12345 + • 



Here 1, 12, 123, 1234, 12345, ... is the Smarandache Consecutive Sequences and 1, 21, 
321, 4321, 54321, ... is the Smarandache Reverse Sequence. 

The interest in Castillo’s article is focused on the calculation of such fractions and their 
possible convergens when the number of terms approaches infinity. The theory of simple 
continued fractions is well known and given in most standard textbooks in Number Theory. 
A very' comprehensive theory* of continued fractions, including general continued fractions 
is found in Die Lehre von den Kettenbruchen [2]. The symbols used to express facts about 
continued fractions vary a great deal. The symbols which will be used in this article 
correspond to those used in Hardy and Wright An Introduction to the Theory of Numbers 
[3]. However, only simple continued fractions are treated in the text of Hardy* and Wright. 
Following more or less the same lines the theory of general continued fractions will be 
developed in the next section as far as needed to provide the necessary tools for calculating 
Smarandache general continued fractions. 



General Continued Fractions 

We define a finite general continued fraction through 



c n =q 0 +- 



' ~ 9o 






?i + <h+ + 






— (i) 



+ " 



?3 + 



q 4 +- 



where {q 0 , qi, q 2 , .... qj and {ij, r 2 , r 3 ... r n } are integers which we will assume to be 
positive. 

The above definition is an extension of the definition of a simple continued fraction where 
r i=r 2 = ... =r n =l. The theory developed here will apply to simple continued fractions as well 
by replacing r k (k=l, 2, ... ) in formulas by 1 and simply ignoring the reference to r k when 
not relevant. 



The formula (1) will usually be expressed in the form 

C n =[q 0 ,qi,q:,q3, -..qn,ri,r 2r r 3 ...r n ] 



( 2 ) 
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For a simple continued fraction we would write 



C n =[qo.q 1 .q 2 ,q 3 , ... 4 ] ( 2 ’) 

If we break off the calculation for m<n we will write 

Cm=[qo,qi,q 2 ,q 3 . ...qm,ri,r 2 ,r 3 ...rj (3) 

Equation (3) defines a sequence of finite general continued fractions for m=l, m=2, m=3, 
.... Each member of this sequence is called a convergent to the continued fraction 

Working out the general continued fraction in stages, we shall obviously obtain expressions 
for its convergents as quotients of two sums, each sum comprising various products formed 
with q 0 , qi, q?, ...q^and r If r 2 , ...r m . 

If m=l, we obtain the first convergent 

Ci = [qo>qi4i] = <7o +— = (4) 



For m=2 we have 



^ r , <?2 r l ?0?1<72 + 4o>* 2 + ?2>i 

C 2 = , „ = 1 ^ , , C>) 



+r 2 



9l<?2 +r 2 



In the intermediate step the value of q l + — from the previous calculation has been quoted, 
putting q 1? q 2 and r 2 in place of q 0 ,, qi and ri. We can express this by 
C2 = [qo,[qi,q2,t2], r i] 

Proceeding in the same way we obtain for m=3 

0?2<73 + '3>1 



( 6 ) 



C3=[<lo,<h,q2,<h,r u r 2 ,r 3 ] = q 0 +- 

<7l<?2<73 + ?i' 3 + ^3 r 2 

?0<7l<?2<?3 + g0^1^3 + <?og3'2 + g2<?3 r l +^l r 3 
? 1 ? 2?3 + ? l '’ 3 + ? 3 r 2 ' 



( 7 ) 



or generally 

Cm - [q0?qi 5 - - - qm-2> [qm-1 ?qjn?t‘ m ]»r* l ,r 2 , . . . T m _i ] (8) 

which we can extend to 

On - [qo?qi ? * • • qm-2y[qm-l ?qm? - * - qn?^m> - - - fn] >^2^2, . . . r m _}] (9) 
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Theorem 1: 



Let Am and B m be defined through 

Ao=q 0 , A^qoqi+rj, A fn =q m A fn .i + r m A in -2 (2^m<n) 

Bo=l, Bi=qj, B m =q m B m . 1 +r m B m . 2 (2<m<n) (10) 

A A 

then C m -tqQ,qi,...q pi ,ri,...r m ] = — — , i.e. — — is the m th convergent to the general 

B m B m 

continued fraction. 



Q Q Aq 

Proof: The theorem is true for m=0 and m=las is seen from [q 0 ]= — = and [q 0 ,qi,ri]= 

1 Bq 

= — . Let us suppose that it is true for a given m<n. We will induce that it is true 
<7i A 
for m+1 

[qo,qi , * * - q^i ,i*i , - - - Tmf i ] “ [qo,qi,-*qm-i,[qm ? c Wi J riirri],ri,...r m ] 



+r m A m _ 2 



(?«+— )4-i +r m A m _ 2 
. ffm+I 

(9m + m_ ' )Bm-\ + r mB m - 2 



- <lm+tenA*-\ + Vin- 

1 B m — 1 ^m-rl^m-1 B m~ 1 



The recurrence relations (10) provide the basis for an effective computer algorithm for 
successive calculation of the convergents C&. 

Theorem 2: 

AoB^^B^.K-ir 1 nr, (11) 

ifc = 1 

Proof: For m=l we have AiBo-BiAo=qoqi+ii-qoqi = ri. 

AmB m .| “BmAm-1 — ( ^”^mA m -2)B m -i “( q^Bm-l "^r^Bm^) Am-I 

*^m(Am-l Bm-2'B 

m-i Am-2) 



42 




By repeating this calculation with m-1, m-2, 2 in place of m, we arrive at 

A ffl B m . I -B m A n , 1 = ...= (A.Bo-B.AoX-ir' U r k =(-l f 1 n r k 

k =2 *=1 

□ 

Theorem 3: 



A m Bm-2'B m A m .2=(- 1 ) ra qm nr* 



( 12 ) 



Proof: This theorem follows from theorem 3 by inserting expressions for A„ and B m 



A TO B m _ 2 *B In A nl _ 2 — ( Qj n A nl _]+r m Ani. 2 )Bni- 2 _ { tJjnBm-l"*’t m B m . 2 )A m .2 



m— I 

Qm( Aro-i Bm.i-Bjn-i A m .2) = (- 1 ) m <lni fl r. 



*=1 



Using the symbol Cn,= — — we can now express important properties of the number 

B m 

sequence m=l, 2, ... , n. In particular we will be interested in what happens to Q, as n 
approaches infinity. 



From (1 1) we have 



(- 1)- 1 nr, 

^n— 1 k=\ 



B n B n _ x B n _ x B n 



while (12) gives 



71-1 



A„ A„_ 



(-D-^-nr* 



C„-C n _ 2 = — - 



k-\ 



B„ B, 



n— 2 



B n - 2 B n 



(13) 



(14) 



We will now consider infinite positive integer sequences {q 0 , qi, q 2 , and {r t , r 2 , ...} 
where only a finite number of terms are equal to 1. This is generally the case for 
Smarandache sequences. We will therefore prove the following important theorem. 

Theorem 4: 

A general continued fraction for which the sequences qo, qi, q 2 , .... and n, r 2 , .... are 
positive integer sequences with at most a finite number of terms equal to 1 is convergent. 

Proof: We will first show that the product B n .]B n , which is a sum of terms formed by- 
various products of elements from {qi, q 2 , ... q,> ri, r 2 , ... r n .;}, has one term which is a 

n 

multiple of I r, . Looking at the process by which we calculated C-, , C 2 , and C 3 , equations 
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4, 5 and 7, we see how terms with the largest number of factors r k evolve in numerators and 
denominators of This is made explicit in figure 1. 





Ci 


C? 


C3 


C4 


C 5 


C* 


C 7 


Cs 


Num. Am 


n 


Oor? 


rir 3 


aoT2r 4 


ri r 3 r 5 


ao w 6 


ur^r? 


OoW^s 


Den. 5m 


- 


h 


qir 3 


r?u 


qir 3 r 5 




qi r3T5T7 


WsTa 



Figure /. The terms with the largest number of r-factors in numerators and denominators. 

As is seen from figure 1 two consecutive denominators BnBn-i will have a term with r 2 r 3 ...r n 
as factor. This means that the numerator of (13) will not cause Cn-Cn_i to diverge. On the 
other hand Bn_iB n contains the term (q^ which approaches oo as n-xc. it follows 

that \im{C n -C n -\) = 0. 

n— fco 

From (14) we see that 

1. If n is odd say n=2k+l, than C^i <C 2 k-i forming a monotonously decreasing number 
sequence which is bounded below (positive terms). It therefore has limit. 

Cik+\ ~ * 

cc 

2. If n is even, n=2k, than C 2 k > C 2 k -2 forming a monotonously increasing number 
sequence. This sequence has an upper bound because C 2 k < C 2 ask-* oc. It 
therefore has limit. 

lim C lk = C 2 . 

3. Since lim \) = 0 we conclude that C\=C 2 Consequently lim C n -C exists. 

n ~ >x 



Calculations 

A UBASIC program has been developed to implement the theory of Smarandache general 
continued fractions. The same program can be used for classical continued fractions since 
these correspond to the special case of a general continued fraction where r!=r 2 = ... =r n =l. 

The complete program used in the calculations is given below. The program applies eq uall y 
well to simple continued fractions by setting all element of the array R equals to 1. 

10 point 10 

20 dim Q(25),R(25),A25),B25) 

30 input "Number of decimal places of accuracy: ";D 
40 input "Number of input terms for R (one more for QJ ”;N% 

50 cis 

60 for i%=0 to N%:recd Q(!%):next 

70 date ‘The relevant data ao, qu ... 

80 for I%=1 to N%:read R(l%):next 

90 data ‘The relevant data for n. r 2 , ... 

100 print tcb(10}:“Smarandcche Generalized Continued Fraction" 

1 10 print tcb(10):‘‘Sequence G:"; 

120 for l%=0 to 6:print Q(!%);:next:print " ETC" 

130 print tab (10) Sequence R:“; 

14C for !%=1 to 6:print R(!%);:next:print " ETC" 
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'Initiating recurrence algorithm 



150 print tcb(10);”Number of decimal places .of accuracy: rt ;D 

160 A(0)=Q(0):B{0)-1 

170 A(1)=Q(0)*Q(1)+R(1):B(1)=Q[1) 

180 Deltas 1:M=1 'M=loop counter 

190 while abs (Delta) >10 A (*D) ‘Convergens check 

200 inc M 

210 A(M)=Q(M)*A(M-1 )+R(M)*A(M-2) ‘Recurrence 

220 B(M)=Q(M)*B(M-l)+R(M)*B(M-2) 

230 Deita=A(M)/B(M)-A(M-1)/B(M-l) 'Cm-Cm-i 

240 wend 

250 print tab(10);*'An/Bn= M ;:print using(2,20),A(M)/B(M) *Cn in decimaiform 

260 print tab(10); M An/Bn=";^rint A(M);7**;B(M) ‘Cn in fractional form 

270 print tcb(10);“Delta=";:print using(2.20), Delta; ‘Delta=Last difference 

280 print " for n-';M *n=number of iterations 

290 print 

300 end 



To illustrate the behaviour of the convergents Cn have been calculated for qi=q?= ...= q^l 
and ri=r 2 = r n =10. The iteration of Cn is stopped when A n = ICn-C^i |<0.01. Table 1 
shows the result which is illustrated in figure 2. The factor (-l)"' 1 in (13) produces an 
oscillating behaviour with diminishing amplitude approaching lim C n -C 

n—^o 

Table 1. Value of convergents Cn for qs{U,...} and reflOJO,...} 






A number of sequences, given below, will be substituted into the recurrence relations (10) 
and the convergence estimate (13). 



Si={l, 1, 1, } . 

S 2 ={1, 2, 1,2, 1,2, } 

S 3 ={3, 3, 3, 3, 3, 3, } 

S 4 ={k 12, 123, 1234, 12345, 123456, } Smarandache Consecutive Sequence. 

S 5 ={1, 21, 321, 4321, 54321, 654321, } Smarandache Reverse Sequence. 

CS1={ 1,1,2,8,9,10,512,513,514,520,521,522,729,730,731,737,738, ... 

NCS1={1,2,3, 4,5, 6,7, 10,11,12,13,14,15,16,17,18, 19,20,21,22,23,24,25,26,29, 30,... 

The Smarandache CS1 sequence definition: CSl(n) is the smallest number, strictly greater 
than the previous one (for n>3), which is the cubes sum of one or more previous distinct 
terms of the sequence. 

The Smarandache NCS1 sequence definition: NCSl(n) is the smallest number, strictly 
greater than the previous one, which is NOT the cubes sum of one or more previous distinct 
terms of the sequence. 

These sequences have been randomly chosen form a large number of Smarandache 
sequences [5]. 

As expected the last fraction in table 2 converges much slower than the previous one. 
These general continued fractions are, of course, very artificial as are the sequences on 
which they are based As is often the case in empirical number theory it is not the individual 
figures or numbers which are of interest but the general behaviour of numbers and 
sequences under certain operations. In the next section we will carry out some experiments 
with simple continued fractions. 



Experiments with Simple Continued Fractions 

The theory of simple continued fractions is covered in standard textbooks. Without proof we 
will therefore make use of some of this theory to make some more calculations. We will first 
make use of the fact that 

There is a one to one correspondence between irrational numbers and infinite 
simple continued fractions . 

The approximations given in table 2 expressed as simple continued fractions would 
therefore show how these are related to the corresponding general continued fractions. 
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Table 2. Calculation of general continued fractions 



Q 


R 


n 


A, 


Q (dec.form) 


Cn (fraction) 


s, 


Si 


18 


-9-10'® 


1.6180339 


6765 












4181 


S2 


s, 


13 


8- 10"® 


1.3660254 


7953 












5822 


S; 


S 3 


22 


-9-10-® 


1.8228756 


1402652240 












769472267 


s 4 


Si 


2 


-7-10 -6 


1.04761 


7063 






3 


5-10’ 12 


1.04761198457 


6742 

30519245 






4 


-2-10' 20 


1.0476119845794017019 


29132203 

1657835914708 

1582490405905 


s 4 


s 5 


2 


-MO' 5 


1.082 


540 






4 


-7-10' 10 


1.082166760 


499 

8245719435 






6 


-MO' 19 


1.08216676051416702768 


7619638429 

418939686644589 1 50004 












387130433063328840289 


s 5 


Si 


2 


-7-10 -6 


1.04761 


7063 






3 


5-10' 12 


1.04761198457 


6742 

30519245 






4 


-2-10' :o 


1.04761198457940170194 


29132203 

1657835914708 

1582490405905 


S 5 


s 4 


2 


-8-1 O' 5 


1.0475 


2358 






3 


7-10' 9 


1.04753443 


2251 

2547455 






5 


1-10' 20 


1.04753443663236268392 


2431858 

60363763803209222 

57624610411155561 


CS1 


NCS1 


6 


-1-10' 7 


1.540889 


1376250 






7 


3*10' 12 


1.54088941088 


893153 

1412070090 






9 


-MO' 20 


1.540889410887887952 55 


916399373 
377447939426190 
24495459359974 3 


NCS1 


CS1 


16 


- 5 - 10' 5 


0.6419 


562791312666017539 












876693583206100846 
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Table 3. Some general continued fractions converted to simple continued fractions 




These sequences show no special regularities. As can be seen from table 3 the number of 
terms required to reach 20 decimals is much larger than for the corresponding general 
continued fractions. 

A number of Smarandache periodic sequences were explored in the author’s book 
Computer Analysis of Number Sequences [6], An interesting property’ of simple continued 
fractions is that 



A periodic continued fraction is a quadratic surd ' i.e. an irrational root of a 
quadratic equation with integral coefficients. 

In terms of An and which for simple continued fractions are defined through 

Ao=q 07 A] =q 0 qi + 1 , A n =q p An .1 +. A^ 

B 0 =l, Bi=qi, B n =qJB n _i+B n _2 (15) 

the quadratic surd is found from the quadratic equation 

BnX^CBn.! - A n )X- An.! =0 (16) 



where n is the index of the last term in the periodic sequence. The relevant quadratic surd is 
A + J7l TbIx -2 A U B^ -^A n _ l B n 



An example has been chosen from each of the following types of Smarandache periodic 
sequences: 

1. The Smarandache two-digit periodic sequence: 

Definition: Let N k be an integer of at most two digits. N k ’ is defined through 
fthe reverse of N k if N k is a two digit integer 

NY H 

l NY 10 if N k is a one digit integer 
N k ~! is then determined by' 

Nm^nvn/I 
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The sequence is initiated by an arbitrary' two digit integer N, with unequal digits. 



One such sequence is Q={9, 81, 63, 27, 45}. The corresponding quadratic equation is 
6210109x 2 -55829745x-1242703=0 

2. The Smarandache Multiplication Periodic Sequence: 

Definition: Let Ol be a fixed integer and N 0 and arbitrary- positive integer. N^, is derived 
from N k by multiplying each digit x of N k by c retaining only the last digit of the product 
cx to become the corresponding digit of N^. 

For c=3 we have the sequence Q={ 1, 3, 9, 7} with the corresponding quadratic equation 
199x 2 -235x-37=0 

3. The Smarandache Mixed Composition Periodic Sequence: 

Definition. Let N 0 be a two-digit integer a r 10+ao. If a,+ao<10 then bi= ai+ao otherwise 
b,= a,+ao+l. bo=|ai-ao| . We define Ni=b r 10+bo. N w is derived from N k in the same way. 

One of these sequences is Q={18, 97, 72, 95, 54, 91} with the quadratic equation 
326258351 5x 2 -58724288064x-645584400=0 

and the relevant quadratic surd 

58724288064 + V f 3456967100707577532096 
6525167030 

The above experiments were carried out with UBASIC programs. An interesting aspect of 
this was to check the correctness by converting the quadratic surd back to the periodic 
sequence. 

There are many interesting calculations to carry- out in this area. However, this study will 
finish by this equality between a general continued fraction convergent and a simple 
continued fraction convergent. 

[ 1 , 12 , 123 , 1234 , 12345 , 123456 , 1234567 , 1 , 21 , 321 ,- 1321 , 54321 , 654321 ]= 

[ 1 , 12 , 5 , 1 , 6 , 1 , 1 , 1 , 48 , 7 , 2 , 1 , 20 , 2 , 1 * 5 , 1 , 2 , 1 , 1 , 9 , 1 , 1 , 10 , 1 , 1 , 7 , 1 , 3 , 1 , 7 , 2 , 1 , 3 , 31 , 1 , 2 , 6 , 38 * 2 ] 
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SOME CONNECTIONS BETWEEN THE SMARANDACHE 
FUNCTION AND THE FIBONACCI SEQUENCE 



"onstantin Durrutrescu I and Carmen Rocsoreanu 



University of Craiova, Dept, of Mathematics 
Craiova 1 100, Romania 



I. INTRODUCTION 

The Smarandache function S:N* ->N* is defined [9] by the condition that S(n) 
is the smallest positive integer k such that k\ is divisible by n. 

If 

n = PT Pz 1X P? (1) 

is the decomposition of the positive integer n into primes, then it is easy to verify that 
S(«)=max(s(p“‘)) (2) 

One of the most important properties of this function is that a positive integer p is 
a fixed point of S if and only if/? is a prime or p = 4. 

This paper is aimed to provide generalizations of the Smarandache function. They 
will be constructed by means of sequences more general than the sequence of the 
factorials. Such sequences are monotonously convergent to zero sequences and divisibility 
sequences (in particular the Fibonacci sequence). 

Our main result states that the Smarandache generalized function associated with 
every strong divisibility sequence (sequence satisfying the condition (15) from bellow) is a 
dual strong divisibility sequence (i.e. it satisfies the condition (26), the dual of (15)). 

Note that the Smarandache function S is not monotonous. Indeed, n x < n 2 does 
not imply S(n x ) < S(n 2 ). For instance 5 < 12 and S( 5) = 5, 5(12) = 4. 

d 

Let us denote by v the least common multiple, by a the greatest common divisor 

d 

and let a = min , v = max . It is known that 

are lattices. The order on N* corresponding to the lattice N 0 is the usual order: 
n x < n 2 o ttj a n 2 -n x 

and it is a total order. On the contrary, the order < corresponding to the lattice N d , 

d 

defined as 

n x <n 2 <=> n, a n 2 — n x 

d d 

( the divisibility relation) is only a partial order. 

More precisely we have 



N 0 = (Af, a, v) and N d = 

V 
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n x <n 2 <=> n x divides n 2 . 

d 

For n x <n 2 we shall also write n 2 >n x . We notice that N d has zero as the greatest 

d d 

element, N 0 does not possess a greatest element and both lattices have 1 as the smallest 
element. Then it is convenient to consider in N 0 the convergence to infinity and in N d , 

the convergence to zero. 

Let 

= n p? and /, 2=n p< i 

be the decompositions into primes of «, and n 2 . Then we have 

-ma x(a it fr) 



n , v fi n 



Up: 



The definition of the Smarandache function implies that 

vn 2 j = 5(/7i)v5(« 2 ) (3) 

Also we have 

n x ^n 2 S{n x )<S{n 2 ). (4) 

d 

In order to make explicit the lattice (so, the order) on the set N * , we shall write 
N 0 instead of N ‘ , if the order on the set of the positive integers is the usual order and 

N d instead of N ’ , if we consider the order < respectively. 



Then (4) shows that the Smarandache function, considered as a function 



S '.N d — » JV 0 , 



( 5 ) 



is an order preserving map. 

From (2) it follows that the determination of S(n) reduces to the computation of 
s{p a ). In addition, it is proved [1] that if the sequence 

( P ): P, P\ P\ - (6) 

is the standard p - scale and the sequence 

[p]: a x {p\ a 2 {p\ ..., a,(p), ... 

is the generalized numerical scale determined by the sequence 

a,(p)- P '~ l 



P~ 1 



then 



S(p a )= p(a b] \ p) . (7) 

In other words, ^(/> £r ) can be obtained by multiplying by p the number obtained 
writing the exponent a in the generalized scale \p\ and "reading" it in the scale (p) . 

For instance, in order to calculate ^(s 99 ) let us consider the scale 
[3] 1,4,13,40,121,... 

Then, for a - 99 , we have 

«[ 3 ] = 2fir 4 (3) + a 3 (3) + a 2 (3) + 2a x (3) = 2 1 1 2 [3] 
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and "reading" this number in the usual scale 
(3) 1, 3, 3 2 , 3 3 , ... 

we get 5 ( 3 ")= 3 ( 2 *3 3 + 3 2 + 3 + 2)=204. So, 204 is the smallest positive integer whose 
factorial is divisible by 3". 

We quote also the following formula used to compute s(p a ): 

s(p a )=(p~ 1 ) a + a [ P ]( a )> ( 8 ) 

where cr^(a) stands for the sum of the digits of the integer a written in the scale \p]. 



2. GENERALIZED SMARANDACHE FUNCTIONS 

A sequence of positive integers is a mapping a : N* — > N* and it is usualy denoted 
by (a„) n€W - (i.e. the set of its values). Since in the sequel an essential point is to 

make evident the structure (the lattice) on the domain and on the range of this function 
respectively, we adopt the notation from (5). 

Then 

& ■ N 0 -*-N d (9) 

shows that a is a sequence of positive integers defined on the set N' . This set was 
structured as a lattice by a and v and its range has also a structure of lattice, induced by 

d 

a and v . 

d 

Definition 2.1. [3] The sequence (9) is a multiplicatively convergent to zero 
sequence ( mez ) if 

(v)neA r * (3) m n ^N" (y)m>m n =>n<<j{m). (10) 

d 

In other words, a (mez) sequence is a sequence defined as in (9), which is 
convergent to zero. 

These sequences, satisfying in addition the condition 

cr(n)<(j{n + \) (11) 

d 

(that is cr(«)divides cr(/? + l)) were considered by G. Christol [3] in order to obtain a 
generalization of p-adic numbers. 

As an example of a {mez) sequence we may consider the sequence defined by 
cr{n) = n\ . This sequence also satisfies the condition (1 1). 

Remark 2.1. We find that the value S{n ) of tbe Smarandache function at the point 
n is the smallest integer m n provided by (10), whenever cr(«)=«!. This enables us to 
define a Smarandache type function for each {mez) sequence. Indeed, for an arbitrary 
{mez) sequence a , we may define S a {n) as the smallest integer m n given by (10). 

The {mez) sequences satisfying the extra-condition (11) generalize the factorial. 
Indeed, if 
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then 



cr(n + \) = k n ^a(n) 



( 12 ) 



<j(n) = k : -k 2 ■ ... ■ k n , with k ] = 1 and k t e N* for / > 1 . 

Starting with the lattices N 0 and N d , we can construct sequences 
a:N d ^N d . (13) 

Definition 2.2. A sequence (13) is called a divisibility sequence (ds) if 
n<m => a(n)<<y(m) (14) 

d a 

(that is if the mapping a from (13) ia monotonous). The sequence (13) is called a strong 
divisibility, sequence (sds) if 

ayi /\m)= a(n)/\a(m) for every n, m e N * . (15) 

Strong divisibility sequences are considered, for instance, by N. Jensen in [5], 

It is known that the Fibonacci sequence is also (sals’) . 

For a sequence cr of positive integers, concepts as (usual) monotonicity, 
multiplicatively convergence to zero, divisibility, have been independently studied by many 
authors. A unifying treatement of these concepts can be achieved if we remark that they 
are monotonicity or convergence conditions of a given sequence o : N' — » N ' , for 
adequate lattices on A r * . 

We shall consider now all the possibilities to define a sequence of positive integers, 
with respect to the lattices N 0 and N d . To make briefly evident the kind of the lattice 
considered on the domain and on the range of a , we shall use the following notation: 

(a) a sequence a ^ : N o — » N o is an ( oo ) - sequence 

( b ) a sequence a od : N a — » N d is an {od ) - sequence 

(c) a sequence a do :N d — » N a is an (do) - sequence 

id) a sequence : N d — » N d is a (dd)- sequence 

We have already seen (Remark 2.1) that, considering ( mez ) sequences, the 
Smarandache function may be generalized. 

In order to generalize tbe Smarandache function for each type of the above 
sequences, it is necessary to consider the monotonicity and the existence of a limit 
corresponding to each of the cases (a) - (d). 

Of course, the limit is infmit for N o -valued sequence and it is zero for the others. 

We have four kinds of monotonicity. 

For a (do)- squence cr do , the monotonicity reads: 

( m do ) (v)/i„7i 2 eW\ n x ^n 2 )<a do (n 2 ) 

d 

and the condition of convergence to infinity is: 

i c do) (v)/!Gtf* (3)m n eN' (vWm„ => (/w) > w . 

d 

Similarly, for a (dd)- sequence a dd , the monotonicity reads: 
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( m &) ( v )w„n 2 eN\ n x <n 2 => cypher 

d d 

and the convergence to zero is: 

i c dd) (V) n<=N' (3)/n„ sN* if)m>m n z^o^m^n. 

d d 

Definition 2.3. The generalized Smarandache function associated to a sequence 
cr y satisfying tbe condition (c (> . ), with i, j e (p, d], is 

Sy (n) = min {m n \ m n given by the condition M (16) 

Remark that (oo)- sequences are the classical sequences of positive integers. As 
examples of (od)- sequences we quote the ( mcz ) sequences. Examples of 
(dd)- sequences are (ds) and (sds)- sequences . Finally, the generalized Smarandache 
functions S od associated with (od)- sequences satisfying the condition (c^) are 
(do) - sequences . 

The functions S I} have the following properties: 

Theorem 2.1. Every function S ^ satisfies : 

(/) (v)/7,, t? 2 eN', », <n 2 ^ S^n^, 

that is satisfies (m oo ) . 

(") S oo (”. V «2 ) = («1 ) V S oo («2 ) 

("0 S oo («. A n 2 ) ' = S oo («1 ) ■ A S oo (»2 ) ' 

Proof: (/) The definition of S^n) implies that: 

Soo («, ) = min { m n t l(V>w > m„ t => (ot) > «, }, for i = 1, 2 
Therefore 

( v ) / ”^^ 00 (« 2 )=> o ' 00 ( m )^«2 ^«1 

and so (/ 7 , ) < S^n^. The equalities («) and (///) are consequences of (/'). 

Theorem 2.2. Every function S od has the following properties : 

(iv) (v)/7,, «, e N', /?, <n 2 ^ )< 

d 

that is S od satisfies (m od ) . 

( v ) ^^iV« 2 J = ^(t7 1 )v^(« 2 ). 

( V/ ) S od (l, A/7, )< («, ) A 5^ (77 2 ) . 

Proof The equality (v) may be proved in the same manner as the equality (3) for 
the function S. Then from (v) it follows (iv ) . 

For (v/) let us note u = S od (//, ) a S od (n 2 ) . From 

77, A/7, <77,, 77, A 77, <77, 
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and from (zy) , it follows that 

S^AZzJ^fo), 

SO ^(z, AH 2 )<S^(njA^(/7 2 ). 

Theorem 2.3. The functions satisfy : 

(v//) (V)w,,w 2 <=N\ n, <n 2 => ^(zz,)^ 

(viii) S do (», v» 2 )<5^ (zz, ) v (n 2 ) . 
fa) S do {n x v« 2 )=5 ito (« ] )v^(» 2 ). 

(*) ^fa 

Proof: Let us note that fa) and (x) are consequences of (vzz). In our terms (vzz) 
is just the fact that the Smarandache generalized function S do associated with a 
{do)- sequence is {po)~ monotonous . To prove this assertion, let zz, <zz 2 . Then for 
every m > m n , we have 

and so S^fafa S do (n 2 ). 

(vzzz) For z = 1, 2 we have: . 

S do fa ) = min \ m nt I ( v ) m \ m n, => fa) ^ ", } 

Let us suppose that zz, < n 2 , so zz, v n 2 = zz 2 and fa v n 2 ) = {n 2 ) . If we take 

d r\ 

m 0 - S do (fh)v S do {n 2 ) , then for every m>m Q it follows that cr do {m) > n s , for z = 1, 2, so 

d 

c do{ m ) - n \ v n , whence the desired inequality. 

Consequence 2.1. S*(w,)aS*(h 2 )< S do (n l )AS do (n 1 ) = S^fa azz 2 )< 

a 

S d o( n ihS do (n 2 ) = S do {n l vn 2 )<S do {n l )vS do (n 2 ). 

Theorem 2.4. The functions S dd satisfy : 

fa) S Jfh v » 2 ] ^fazfavS^fa). 

V J 

(xzz) If n, <n 2 or n 2 <zz, then 

d d 

v « 2 j = ^(«,)v^(n 2 ). 
fa«) fa, (f*l A ^2 )~ ^ dd (P] ) A ^ dd (d 2 ) ‘ 

Proof : The proof of fa) is similar to the proof of (vzzz) and the other assertions 
may be easily obtained by using the definition of S ^ from (17) (for z =j = d). 
Consequence 2.2. For all n x ,n 2 s N ’ we have 
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s dd (»1 ) V S* i n 2 ) £ [n x vn 2 < S„ (w, )v^ ( w 2 ). 



This follows from the fact that 

«, <«, vn 2 for i = 1, 2 => S* («, ) < S^f/7, vw 



If cr^ is a divisibility sequence, the above theorem implies that the associated 
Smarandache function satisfies the inequality (x/) . In the following we shall see that, if the 
sequence a ^ is a divisibility sequence with additional properties, namely if it is a strong 
divisibility sequence, then the inequality ( xi ) becomes equality. 

Theorem 2.5: If <j dd is a (sds) satisfying the condition (c^), then : 

S<id { n \ v 1 = S dd («, ) v S dd (n 2 ) (17) 



( v )«,,« 2 gN \ n ] <n 2 => S& 

d d 

(i.e. S dd satisfies the monotonicity condition (m dd )) . 

Proof: In order to prove the equality (17), it is sufficient to show that 

S d d(n,)<Sdd( n \ vi^, for i =1,2. 



But if, for instance, the above inequality does not hold for n, and we denote 



d o = S d d { n x)^S dd \n x wn 2 y 

it follows that d 0 < S dd (n, ) and taking into account that 

°dd{Sdd( n i))>«i and «, <«, v « 2 ^aJsJn, vn 



1 » 



we have 



f f d X) 

Mo)=°dd S dd{ n x)f}S d d «, V «2 = 

V d K JJ 

f f d f] 

T dd\ S dd{ n x))^dd $dd «1 V «2 >»,A »,=«, 

" k v )) d d 



Thus, we obtain the contradiction 
S dd {n l )<d Q <S tU (n I ). 

So, if the sequence a dd is a (sz&), that is if the equality (15) holds, then the 
corresponding Smarandache function S dd satisfies the dual equality (17). 

Example. The Fibonacci sequence (F n ) neN - is a (sds). Therefore, the generalized 
Smarandache function S F associated with this sequence satisfy: 
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f a \ a 

S F n,vn 2 =S F {n l )vS F (n 2 ) 



By means of this equality, the computation of S F (n) reduces to the determination 
of S F (p a ), where p is a prime number. For instance 

^ (52) = min I m n | (y)m>m n => 52<.F(/w)j= 

' d d ' 

= S f {2 2 )vS f {\3)=6v7 = 42. 

So, 42 is the smallest positive integer m such that F(m) is divisible by 52. 

Also, we have 

^f(12)- *SV(2 2 '3)= *S f (2 2 )vS f ( 3)= 6v4 = 12, (20) 

therefore n =12 is a fixed point of S F . 

The values of S F {p a ) may be obtained by writing all F n in the scale (p) given by 
(6), which is a difficult operation. At the time being, we are not able to provide a closed 
formula for the computation of S F (p a ) • However, we shall present some partial results in 
this direction. In [8] it is stated that 



3* <f.o4-3 H 



2 ‘<F„o 3-2 h <», for k >3 . 

d d 

It is known (see for instance [6], [7]) that if cr is a non-degenerate second-order 
linear recurrence sequence defined by 

<j{n)= Ao{ti-\)- Ba{n-2^) (21) 

where A and B are fixed non-zero coprime integers and ct(i)= 1, <r(2 )= A , then 

neZ*, nr\B = 1 => (3)/n e N' n<<j(m). (22) 

d d 

The least index of these terms is called the rank of appearance of n in the sequence 
and is denoted by r(ri). 

If D = A 2 -4B and {Din) stands for the Jacobi symbol, then for mn/\BD = 1 and 

d 

p a prime we have ([6]) 

n <a(m) <^> r{n)<nr, r(p)<p-{DI p) 

d d d 

U^ PziP'p) o (Blp) = 1; r [ mini = r{m)vr(n). ^ 



Let us denote N* B = jn e N * | n a5 = 1 j. Obviously, if r is considered as a function 

r : N" b — > N’ , then we can write: 
r(n) = min vm \ n<a{m)\. 



Whence an evident parallel between the above methods described for the construction of 
the generalized Smarandache functions and the definition of the function r. 

For the Fibonacci sequence (F n ) we have A = 1, B = - 1 and so D = 5. 

This implies 

/? = 5£±1=>(5//>)=1 (24) 
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p = 5k±2=>(5lp)=-\ (25) 

and it follows that if (24) holds, then p divides . Thus S F (p) is a divisor of p - 1 . In 
the second case p divides F p+X and S F (p) is a divisor of p+ 1. 



From (23) we deduce 

Sf(p)^P~(5/p) 

for any prime number p. 

Lemma 2 from [6] implies that the fraction (p-( 5/ p))/S F (p) is unbounded. We 
also have 

P l <F. <=> S r {p l ) r . 

a a 

Example. For p = 1 1 it follows (5/p) = 1, so S F ( 1 1)<10 . In fact, we have precisely 

d 

S F (l l) = 1 1 - (5 / 1 1) = 1 0 , but there exist prime numbers such that S F (p) < p - (5 / p) . For 
instance,/? =17, for which p - (5/p) =18 and 5 F (l7) = 9 . 

Definition 2.4. The sequence cr is a dual strong divisibility sequence (dsds) if 



nvm 
V J 



d 

= a(n)va(m) for all n,m<=N * . 



(26) 



It may be easily seen that every strong divisibility sequence is a divisibility 
sequence. We also have: 



Proposition 2.1 Every dual strong divisibility sequence is a divisibility sequence. 
Proof. We have to prove that (26) implies (14). But if n<m, it follows 

d 

d 

n\zm-m and then 



( a \ d 

< 7 ( 171 )= <y nvm = a(n)vo(m) 



so, o(n)<<j(m) . 



K 



(27) 



Then Theorem 2.5 asserts that the Smarandache generalized function S a 
associated with any strong divisibility sequence cr is a dual strong divisibility sequence. Of 
course, in this case, both sequences cr and S a are divisibility sequences. 

It would be very interesting to prove whether the converse assertion holds. That is 
if S dd is the generalized Smarandache function associated with a ( divisibility) sequence 

< 7 dd satisfying the condition (c^), then the equality (17) implies the strong divisibility. 



Remarks. (1) It is known that the Smarandache function S is onto. But given a 
(dd)- sequence <7 dd , even if it is a (sds), it does not follow that the associated function 

S dd is onto. Indeed, the function S F associated with the Fibonacci sequence is not onto , 
because n = 2 is not a value of S F . 
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(2) One of the most interesting diophantine equations associated with a 
Smarandache type function is that which provides its fixed points. We remember that the 
fixed points for the Smarandarche function are all the primes and the composit number 
n = 4. For the functions the equation providing the fixed points reads S dd ( x )= x and 

for S F we have as solutions, for instance, n - 5, n =12. 

At the end of this paper we quote the following question on the Smarandache 
function, also related to the Fibonacci sequence: 

T. Yau [10] wondered if there exist triplets of positive integers («, n- 1, n- 2) such 
that the corresponding values of the Smarandache function satisfy the Fibonacci 
recurrence relation S(n) = S(n - 1) + S(n - 2). 

He found two such triplets, namely for n =1 1 and for n =121. Indeed, we have 
5(9) + 5(10) = 5(11) and 5(119) + 5(120) = 5(121). 

Using a computer, Charles Ashbacher [2] found additional values. These are for 
n = 4902, n = 26245, n = 32112, n = 64010, n = 368139, n = 415664. 

Recently H. Ibsent [4] proposed an algorithm permitting to find, by means of a 
computer, much more values. But the question posed by T. Yau "How many other 
triplets with the same property exist?"is still unsolved. 

Acknowledgements 

We wish to express our gratitude to Professor Adelina Georgescu for her 
important contribution to the English translation of this paper. 

REFERENCES 

1. M. Andrei, C. Dumitrescu, V. Seleacu, L. Tutescu, St. Zanfir: "Some Remarks on the 

Smarandache Function", Smarandache Function Journal, Vol. 4-5, No. 1, (1994), pp. 1-5. 

2. Ashbacher, Ch.: " Some Problems on Smarandache Function ", Smarandache 

Function Journal, Vol. 6, No. 1, (1995), pp. 21-36. 

3. Christol, G.: "p-adic Numbers and Ultrametricity", From Number Theory to Phisics, M. 

Woldschmidt et all Editors, Springer-Verlag, 1992, pp. 440-475. 

4. Ibsent, H: "Bases Solutions (The Smarandache Function)", Smarandache Notion Journal, Vol. 

7, No. 1, (1996),86-89. 

5. Jensen, N.: "Some Results on Divisibility Sequences", Appl. of Fibonacci Numbers, 

Vol. 3, G. E. Bergum et al. Editors, Kluwer Academic Press, 1990, pp. 181-189. 

6. Kiss P., Phong B. M., Lieuwens E.: "On Lucas Pseudoprimes Which Are Products of s 

Primes", Fibonacci Numbers and Their Applications, Vol. 1, A. N. Philippou et al. 
Editors, D. Reidel Publishing Company, 1986, pp. 131-139. 

7. Kiss, P.: "Primitive Divisors of Lucas Numbers" Appl. of Fibonacci Numbers, Vol. 2, A. N. 

Philippou et al. Editors, Kluwer Academic Prees, 1988, pp. 29-38. 

8. Robbins N.: "Fibonacci Numbers of the Forms PX~ ±1, PX 3 ±1, Where P is a Prime", 

Appl. of Fibonacci Numbers, Vol. 2, A. N. Philippou et al. Editors, Kluwer Academic 
Press, 1988, pp. 77-88. 

9. Smarandache, F.: "A Function in the Number Theory", An. Univ. Timisoara, Ser. St. Mat. 

Vol. XVIII, fasc. 1, (1980), pp. 79-88. 



59 




10. Yau, T.: "A Problem Concerning the Fibonacci Series", Smarandache Function Journal, Vol. 
4-5, No. 1, (1984), pp, 42. 

AMS Classification Numbers: 11 A 25, 11B39. 



60 




COMPUTATIVE PARADOXES IN 
MODERN DATA ANALYSIS 



Y. V. CHEBRAKOV AND V.V.SHMAGIN 
Department of Mathematics, Technical University, 
Nevsky 3-1 1, 191186, St-Petersburg, Russia 
E-mail: chebra@phdeg.hop.stu.neva.ru 



By developing F. Smarandache thema on paradoxes in mathematics it is stated, 
firstly, if in measurement (natural science) experiments the best solutions are found 
by using methods of modem data analysis theory, then some difficulties with the 
interpretation of the computation results are liable to occur; secondly, one is not 
capable to overcome these difficulties without a data analysis theory modification, 
consisted in the translation of this theory from Aristotelian “binary logic” into more 
progressive “fuzzy logic”. 

Key words: data analysis, revealing outliers, confidence interval, fuzzy logic. 



I Introduction 

As generally known from history of science, a scientific theory may have crisis in 
process of its development, when it disjoints in a set of fragment theories, that weak- 
coordinate each other and, as a whole, form a collection of various non-integrated 
conceptions. For instance, as we assume, F. Smarandache mathematical notions and 
questions 1 ~ 2 help us to understand quite well that a stable equilibrium, observed in 
mathematics at the present time, is no more than fantasy. Thus, it falls in exactly 
with F. Smarandache views that the finding and investigating paradoxes in 
mathematics is a very effective way of approximating to the truth and so at present 
each of scientific researches, continuing F. Smarandache thema 2 , should be 
considered as very actual one. 

Let us assume that computative paradoxes in mathematics are mainly such 
computation results, obtained by using mathematical methods, which are 
contradicted some mathematical statements. The main goal of this paper is to 
demonstrate that the mentioned crisis, demanding practical action instead of debate, 
occurs in modem data analysis, which formally has its own developed mathematical 
theory, but does not capable “to cope worthily” with a large number of practical 
problems of quantitative processing results of measurement experiments. 

Another goal of this paper is to equip the mathematicians and software 
designers, working in the data analysis field, with a set of examples, demonstrating 
dramatically that, if, for solving some problems on analysing data arrays, one uses 
the standard computer programmes and/or time-tested methods of modem data 
analysis theory, then a set of the paradoxical computative results may be obtained. 
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2 Approximative problems of data analysis 

2.1 The main problems of regression analysis theory and 
standard solution methods 

As generally known 3-7 , for found experimental dependence {y„, x„} (n = 1, 2, 

N) and given approximative function F( A, x), in the measurement (natural science) 
experiments the main problems of regression analysis theory are finding estimates of 
A' and y ' and variances of 5A' and 5(y-y') , where A' is an estimate of vector 
parameter A of the function F(A,x) and {y„'}= {F(A',x„)}. In particular, if 
F(A, x) = X^a/fyfx) {F(A, x) is a linear model}, where hfe) are some functions 

on x, then in received regression analysis theory standard solution of discussed 
problems has form 

A' = (H T H) -1 H t Y , (8A') 2 = s/(N-L) diag(H T H) -1 (1) 

8 'p(y-y')=y , ±tpsyl i + h i t (h t h) -1 h i - , 

where H is a matrix LxN in size with n-th row (hfx„), /j 2 (x„), ... , /zz.(x„)); H T is the 
transposed matrix H; Y= {y„}; s = l% =1 (y n ~ yhf KN - L ) ; H,= (hfx,), h 2 (x,), 

... , h L (x,))] the value of t p is determined by t-Student distribution table and 
generally depends on the assigned value of the significance level of p and the value 
of N-L (a number of freedom degree); at the assigned value of the significance 
level of p the notation of h p (y -y ') means confidence interval for possible 
deviations of experimental values of y from computed values y'= F(A', x,). 
According to Gauss - Markov theorem 4 ’ 5 , for classical data analysis model 

y„ = F( A, x„) + e„ (2) 

the solution (1) is the best (gives minimum value of s), if the following conditions 
are fulfilled: 

all values of { x„ } are not random, mathematical expectation of random value 
{e„} is equal to zero and random values of {e„} are non-correlated and have the 
same dispersions o 2 . 

Example 1. In table 1 we adduce an experimental data array, obtained by 
Russian chemist D.I.Mendeleev in 1881, when he investigated the solvability (y, 
relative units) of sodium nitrate (NaNCb) on the water temperature (x, °C). 



n 




/n 


Yn-Yn 


i 


0 


l 


-0.80 


2 


4 


71.0 


0.02 


3 


10 


76.3 


0.10 



D.I.Mendeleev data array 



n 




Yn 


1 

> 


4 






0.05 


5 


21 


85.7 


-0.07 


6 


29 


92.9 


0.17 



Table 1. 



n 




Yn 


Yn ~ Yn 


7 


36 


99.4 


0.58 


8 






1.73 


9 


68 


125.1 


-1.56 
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By analysing the data array {y„, x„}, presented in table 1, Y.V.Linnik 3 states that 
these data, as it was noted by D.LMendeleev, are well-fitted by linear model 
y' = 67.5 + 0.871 x (5A' = (0.5; 0.2)), although the correspondence between 
experimental and computed on linear model values ofy is slightly getting worse at 
the beginning and end of investigated temperature region (see the values {y„ -y'„} 
adduced in table 1). We add that for discussed data array Y.V.Linnik 3 computes the 
confidence interval of 5^ (y -y') from (1) at the significance level of p = 0.9: 

§ 0.9 (y~y') = ±0.593 Vl + U-26) 2 /451J . (3) 



AV 




Figure 1 . The plots of confidence interval of the deviation of y from y' (heavy lines) and 
residuals y-y' (circles) for D.LMendeleev data array. 



We show the plots of 5o.9 (y-y') °n x by heavy lines in figure 1 and {y„ -y„', 
x„ } by the circles. Since the plot of {y„ -y„'x„} steps over the heavy lines in figure 
1, some computative difficulty is revealed: 

the standard way (1), used by Y.V.Linnik 3 for determining the confidence 
interval of the deviations of y from y' is out of character with the discussed 
experimental data array. 

It follows from results presented in table 1 and/or figure 1, if one assumes that 
5 (y-y') ^ max ly„ —yf I = 1.73 then the broken connections of the confidence 
interval S(y -y') with D.LMendeleev data array will be pieced up. But values of 
8A', calculated by Y.V.Linnik from (1), disagree with the values 8(y- v') > 1.73, 
and, consequently, 

standard values of & 4 'is out of character with D.LMendeleev data array also. 

2.2 Alternative methods of regression analysis theory 

P. Huber noted that, as the rule, 5 - 10% of all observations in the majority of 
analysing experimental arrays are anomalous or, in other words, the conditions of 
Gauss - Markov theorem, adduced above, are not fulfilled. Consequently, in practice 
instead of the standard solution (1), found by “least squares (LS) method”, 
alternative methods, developed in the frames of received regression analysis theory, 
should be used. In particular, if the data array (y„, x„} contains a set of 
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outliers, then for finding the best solution of discussed problem it is necessary 6, 7 or 
to remove all outliers from the analysing data array {strategy 1 ), or to compute the 
values of A' on the initial data array by means of M-robust estimators ( strategy 2). 
For revealing outliers in the data array P.J.Rousseeuw and A.M.Leroy 9 suggest to 
use one of two combined statistical procedures, in which parameter estimates, 
minimising the median of the array {(y„ -y n r ) 2 } ( the first procedure) or the sum of K 
first elements of the same array ( the second procedure), are considered as the best 
ones. If F{ A, x) is a linear function (see above), then the robust M-estimates of A' 
are obtained as result of the solving of one from two minimisation problems 6-9 
N N 

5 9 (A) = I<p()> n -^)^mi n or dSy/da[= '£\\r(y n -y' n ) h l (x n ) = 0, (4) 

n - 1 n-l 

where function cp(r) is symmetric concerning Y-axis, continuously differentiable 
with a minimum at zero and cp(0) = 0; \jr(r) is a derivative of cp(r) with respect to r. 

Continued example 1. Since D.I.Mendeleev data array from table 1 contains 
outliers, we adduce results of quantitative processing this data by alternative 
methods, defined above. 

1. Let in (4) Andrews function 10 be applied: cp(r) = d(\—cos(r/d)) if \r\<dK 
and tp(r) =0 if I r \ > dn . It is articulate in figure 2 that in this case the values of the 
linear model parameters ao and a\ depend on 

a) the values of parameter d of Andrews function cp(r); 

b) the type of the minimisation robust regression problem (solutions of the first 
and second minimisation problem of (4) are marked respectively by triangles and 
circles in figure 2). 

Thus, in this case a computative paradox declares itself in the fact, that 

in actual practice the robust estimates are not robust 

and so, as K.R.Draper and H. Smith 11 wrote already, 

“unreasoning application of robust estimators looks like reckless application of 
ridge-estimators: they can be useful, but can be improper also. The main problem is 
such one, that we do not know, which robust estimators and at which types of 
supposes about errors are effectual to applicate; but some investigations in this 
direction have been done...” 





Figure 2. Dependences of parameters values of linear model a 0 +a^x on values internal 
parameter of robust Andrews estimator and the type of the minimisation problems (4). 
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2. Let us reveal outliers in D.I.Mendeleev data array by both combined statistical 
procedures 9 , mentioned above. 

Our computation results show 

a) both procedures could not find the all four outliers (1,7, 8 and 9) but the 
only three ones with numbers 1, 8 and 9; 

b) if a set of readings with numbers 1 , 8 and 9 is deleted from D.I.Mendeleev 
data array, then for the truncated data array the first procedure will not find a new 
outlier, but the second procedure will find two outliers yet that have numbers 2 and 
6 in the initial data array. 

Thus, in this case the main computative paradox is exhibited in the fact, that 

revealing outliers problems solutions depend on a type of the used statistical 
procedures. 

It remains for us to add, if one 

a) computes y by formula 6> 7 

/£, x) = g a ( 68. 12 + 0.02^ + (0.85652 - 0.00046^) x ), (5) 

then for each n the difference I y„ - y„ I will keep within the limit of the chosen 
above value for the confidence interval 5(y — y'}, where a=0.94; 0 < £, < 35; g a (y) — 
2a[y/(2a)] + 2a at I y - 2a [y/( 2a)] I > a, otherwise g a (y) - 2a [y/( 2a)], [ b ] means 
integer part of b. Thus, another computative paradox occurs: 

although for each contaminated data array a family of analytical solutions 
exists, the only single solution of the estimation problems is found in modem 
regression analysis theory. 

b) puts the mentioned above extremal values of £, in (5), one will be able to 
determine the exact limit of the variation for the linear model parameters ao and a\. 
a 0 = 67.77 ± 0.35 and a x = 0.865 ± 0.008; 

c) deletes a set of readings with numbers 1, 7, 8 and 9 from D.I.Mendeleev data 
array, one will obtain that in the truncated data array {y„,x„}* the difference of 
I y„ -y» I for each n keeps within the limit of the error e, where £ is the measuring 
error for readings {y„}*: £ = 0.1. Since in this case 5 (y -y') < £, the complete family 
of analytical solutions has form 6 ’ 7 

/(£, x) = g a ( 67.566 + 0,002q + (0.870047 - 0.000097^) x ), (6) 

where a=0.07; 0<q<45 and, consequently, ao = 67.521 ±0.045 and a\ = 
0.872 ±0.002; 

d) compares solutions (5) and (6) with the standard LS-solution, one can 
conclude that LS-estimations of parameter ao and a\ (A' = (67.5±0.5; 0.87±0.2)} 
are pretty near equal of the mean values of these parameters in the general analytical 
solutions (6) and (7). However, 

values of variances 5a o' and 5a/, computed by standard method, disagree with 
exact values determined by (5). 
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2.3 The main paradox of regression analysis theory 

As it emerges from analysis of information presented in Sect. 2.2, the main paradox 
of modem regression analysis theory is exhibited in a contradiction between this 
theory statements, which guarantee uniqueness of data analysis problems solution, 
and multivarious solutions in actual practice. In this section we adduce yet several 
computative manifestations of this paradox. 

Example 2. In table 2 a two-factors simulative data array is presented. 



Table 2. 



Simulative data array 



n 


Xn 


y n 


n 


Xn 


Yn 


n 


Xn 


Yn 


1 


- 1.0 


0.50 


8 


- 0.3 


0.75 


15 


0.4 


0.87 


2 


- 0.9 


0.55 


9 


- 0.2 


0.77 


16 


0.5 


0.89 


3 


- 0.8 


0.59 


10 


- 0.1 


0.79 


17 


0.6 


0.90 


4 


- 0.7 


0.63 


11 


0.0 


0.81 


18 


0.7 


0.91 


5 


- 0.6 


0.66 


12 


0.1 


0.83 


19 


0.8 


0.92 


6 


- 0.5 


0.69 


13 


0.2 


0.84 


20 


0.9 


0.93 


7 


- 0.4 


0.72 


14 


0.3 


0.86 


21 


1.0 


0.94 



Let the approximative model have form 

y— (<2o + a] x + a 2 x 2 )/( 1 + aye + ayx 2 ) . ( 7 ) 

To find vector parameter A estimates of the model (7) on the data array from 
table 2 we use two different estimation methods. As the first method we choose the 
estimation one, involved in the software CURVE-2.0, designed AISN. In this case 
we obtain, that 

A'= {0.81; 0.008; -0.31; -0.22; -0.24). 

As the second estimation method we select Marquardt method 12 . Using the value 
A', found above by the first estimation method, as initial value of A we obtain that 
in the second case 

A'={0.81; 0.55; 0.035; 0.45; 0.34}. 

Thus, 

values of A', obtained by two different estimation methods, differ from each 
other. 

Example 3. In table 3 yet one two-factors data array is presented. Let us select 
the model y = a\x + e as approximative one and assume, that v is the random 
variable with the known density function p: 

p= exp{-0>-a 1 x)/(2/(a 2 ))}/V 27r /( a 2) » ( 8 ) 

where f(a 2 )= a) ay, b) a 2 x 2 ; c) aye. 
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We will find estimates of parameters a i and <22 by method of maximum likelihood. 
L= nexp{-(>’„ -a^n)/ (2/ (a 2 ))l Jhtfia^j => max (9) 

n - 1 

or din L ! eb,= 0, where symbol “In” means the natural logarithm. 



Table 3. 

Two-factors data array 



n 


*n 




n 


X, 


Yn 


n 




Yp 


1 


5 


21 


6 


11 


53 


11 


16 


81 


2 


6 


31 


7 


11 


56 


12 


19 


97 


3 


8 


38 


8 


12 


60 


13 


20 


98 


4 


8 


37 


9 


14 


68 


14 


22 


107 


5 


10 


53 


10 


15 


72 


- 




- 




a) b) c) 

Figure 3. Dependences {y„ - a/x n } for different hypothesis about Law for the random 

variable y variance. 



Computation results of V.I.Mudrov and V.P.Kushko 13 show, that in the 
discussed case the estimates values of parameters a\ depend on hypothesis about 
Law for the random variable y variance: for case (a) in (8) a\ = 4.938 (£' = 

4.995- 1(T 14 ); for case (b) a{ = 4.896 ( L ' = 4.421- 1(T 16 ) and for case (c) < 21 ' = 4.927 
( L' = 9.217- 10~ 15 ). By analysing obtained results the authors 13 conclude, that, since 
likelihood function (9) has maximum values for case (a), the more likelihood 
hypothesis about Law for the random variable y variance is the hypothesis (a): 
variance of y is the constant value. 

We demonstrate in figure 3 that for cases (a), (b) and (c) dependences 
A y- {e„} = {y„- a i'x„} have practically the same form and, consequently, 

the strong distinction of values L' for all mentioned cases does not tread on 
infirm ground. 

It should be noted that 

- the very apparent expression of the discussed main computative paradox of 
regression analysis theory one may find also in books 6, 7 ’ n , where, for the problem 
on finding the best linear multiple model, fitting Hald data array, a set of solutions. 



67 




found by various procedures and statistical tests of modem regression analysis 
theory, is adduced; 

- the most impressive formulation of the main paradox of regression analysis 
theory is contained in Y.P.Adler introduction 14 : 

“When the computation had arisen, the development of regression analysis 
algorithms went directly «up the stairs, being a descending road». Computer was 
improving and simultaneously new more advanced algorithms were yielded: whole 
regression method, step-by-step procedure, stepped method, etc., - it is impossible 
to name all methods. But again and again it appeared that all these tricks did not 
allow to obtain a correct solution. At least it became clear that in majority cases the 
regression problems belonged to a type of incorrect stated problems. Therefore 
either they can be regularised by exogenous information, or one must put up with 
ambiguous, multivarious solutions. So the regression analysis degraded ingloriously 
to the level of a heuristic method, in which the residual analysis and common sense 
of interpreter play the leading role. Automation of regression analysis problems 
came to a dead-lock”. 



3 Data analysis problems at unknown theoretical models 

Let us assume, that a researcher is to carry out a quantitative analysis of a data array 
{X„} in the absence of theoretical models. Further consideration will be based on the 
fact 6, 7 that the described situation demands a solution of following problems 

- verification of the presence (or absence) of interconnections between analysed 
properties or phenomena; 

- determining (in the case when the interconnection is obvious, a priori and 
logically plausible) in what force this interconnection is exhibited in comparison 
with other factors affecting the discussed phenomena; 

- drawing a conclusion about the presence of a reliable difference between the 
selected groups of analysed objects; 

- revealing object’s characteristics irrelevant to analysed property or 
phenomenon; 

- constructing a regression model describing interconnections between analysed 
properties or phenomena. 

In following sections we consider some methods allowing to solve foregoing 
problems. 

3.1 Correlation analysis 

When one is to carry out a quantitative analysis of the data array {X„} in the absence 
of theoretical models, it is usual to apply correlation analysis at the earlier 
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investigation stage, allowing to determine the structure and force of the connections 
between analysed variables 15 “ 11 . 

Let, for instance, in an experiment each n-th state of the object be characterised 
by a pair of its parameters y and x. If relationship between y and x is unknown, it is 
sometimes possible to establish the existence and nature of their connection by 
means of such simple way as graphical. Indeed, for realising this way, it is sufficient 
to construct a plot of the dependence {y„, x„} in rectangular coordinates y-x. In this 
case the plotted points determine a certain correlation field, demonstrating 
dependences x = x(y) and/or y = y(x) in a visual form. 

To characterise the connection between y and x quantitatively one may use the 
correlation coefficient R, determined by the equation 15 - 17 

s.».,(,.-yx*-») . (10) 

where y and x are the mean values of parameters y and x computed on all N 
readings of the array {y„,x n }. It can be demonstrated that absolute value ofi?^ does 
not exceed a unit: —1 < Ryx< 1. 

If variables y and x are connected by a strict linear dependence y = a 0 + a x x, then 
R >x = ±1, where sign of Ryx is the same as that of the a\ parameter. This can follow, 
for instance, from the fact that, using Ryx, one can rewrite the equation for the 
regression line in the following form 15 - 17 

y= y +Ryx(Sy/S x ) (x-x), (11) 

where S y and S x are mean-square deviations of variables y and x respectively. 

In a general case, when -1 <Ryx< 1, points {y„, x„} will tend to approach the line 
(11) more closely with increasing of \Ryx lvalue. Thus, correlation coefficient (10) 
characterises a linear dependence of y and x rather than an arbitrary one. To 
illustrate this statement we present in table 4 the values Ryx = Ryfo.) for the 
functional dependence y = x a , determined on x-interval [0.5; 5.5] in 11 points 
uniformly. 



Table 4. 

The values Ryx = Ryx(o) for the functional dependence y=x“, 
determined on interval [0.5; 5.5] 



a 


^y*( a a ) 


Ryx(a = a) 


a 


Ryx( a = -a) 


Ryx(a = a) 


3.0 


-0.570 


0.927 


1.0 


-0.795 


1.000 


2.5 


-0.603 


0.951 


0.5 


-0.880 


0.989 


2.0 


-0.650 


0.974 


0.0 


0.0 


0.0 


1.5 


-0.715 


0.992 




- 


- 
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Let us clear up a question what influence has the presence of outliers in the data 
array {y n ,x„} on the value of correlation coefficient (10). To perform it let us 
analyse a data array 

{x„}=(-4;-3;-2;-l;0;10), (12) 

{y n ) = (2.48; 0.73; -0.04; -1.44; -1.32; 0), 

where, on simulation conditions 8 , the reading with number 6 is an extremal outlier 
(such reading that contrasts sharply from others); approximative function F(A, X) = 
ao + a\x and A^e = (-2; -1). 

By computing the values of Ryx of (10) and s of (1), we determine the number i 
of a reading, which elimination from this data array leads to the maximum absolute 
value of and, consequently, to the minimum value of s (the most simple 
combinatoric-parametric procedure Ps, allowing to find one outlier ’ in a data 
array). Our calculations show that the desirable value Ryx = — 0.979 and i— 1. We 
note, if extremal outlier ye is removed from the array (12), Ryx = —0.960, but s of (1) 
takes the minimum value. Presented results enable us to state that 

procedure Ps loses its effectiveness when revealing the outlier is made not by 
test s, but by test Ryx. 

Let us consider another case. For the array (12) the noise array {e n } — {-2- x n - 
y n } = (-0.48; 0.27; 0.04; 0.44; -0.68; -12.0). We reduce by half the first 5 
magnitudes of the noise array { e „ }: {e*} D ew = (—0.24; 0.14; 0.02; 0.22; -0.34; —12.0); 
form a new array {y n } new = {-2- x„- (e„) aew } and determine again the number i of a 
reading, which elimination from the data array {y n , x„} ne w leads to the maximum 
absolute value of Ryx. In the described case Ryx reaches its maximum absolute value 
when the reading 6 (extremal outlier) is deleted from the array {y„, x„} new (Ryx — - 
0.989). If from the array (y„, x n } new we eliminate the reading 6, identified correctly 
by the test “the maximum absolute value of Ryx”, then by this test we are able to 
identify correctly the sequent outlier (the reading 5) in the discussed array. Thus, we 
obtain finally 

when a dependence between the analysed variables is to a certain extent close to 
a linear, one may use the correlation coefficient (10) for revealing outliers, 
presented in data arrays. 

It is known 15 -!7 , when the number of analysed variables K > 2, the structure 
and force of the connections between variables x\, xi_, ... , xk are determined by 
computing all possible pairs of correlation coefficients R x . x from (10). In this case 

all coefficients R x . x . are usually presented in the form of a square symmetric K by 

K matrix: 




which is called a correlation matrix (we note that in this matrix diagonal elements 
Ra= 1). Finding strong-interconnected pairs of variables x\, X2, ... , xk on the 
magnitudes of coefficients Rg from the matrix R is a traditional use of matrix (13) in 
data analysis. But, obviously, 

using the mentioned way, one should bear in mind all ideas presented above in 
outline concerning the correlation coefficient (10). 

3 . 2 Discriminant analysis 

Let a certain object W be characterised by a value of its vector parameter X w = (xi, 
X2 , ... , x K ) ; W u W 2 , ... , W p be p classes and the object W must be ranged in a class 
Wj on the value of its vector parameter X w . In discriminant analysis the formulated 
problem is the main one 18 ~ n . 

The accepted technique for solving the mentioned problem entails construction 
of a discriminant function D{ A, X). A form and coefficients {a,} 1 , 2 , values of 

this function are determined from the requirement, that values ofZ)(A, X) must have 
maximum dissimilarity, if parameters of objects, belonging to different populations 
W\, W 2 , ... , W p , are used as arguments of this function. 

It seems obvious that in a general case, firstly, D( A, X) may be either linear or 
non-linear function on {a,} and, secondly, must be some connection between the 
problem-solving techniques of discriminant and regression analyses. In particular, as 
stated 18 ~ 2 \ for solving problems of discriminant analysis one may use standard 
algorithms and programs of regression analysis. Thus, the similarity of techniques, 
used for solving problems of the regression and discriminant analyses, makes it 
possible in discriminant analysis to apply alternative algorithms and procedures of 
regression analysis and, consequently, 

if data analysis problems are solved by discriminant analysis techniques then in 
practice the researcher may meet the same difficulties which are discussed in 
Sect. 2. 

3.3 Regression analysis 

In the absence of theoretical models it is usual to employ regression analysis in 
order to express in a mathematical form the connections existing between variables 
under analysis. 

It happens with extreme frequency that researchers impose limitations on a type 
and form of approximative models or, in other words, approximative models are 
often chosen from a given set of ones. Evidently, in this case it is required to solve 
problem on finding the best approximative model from a given set of models. Let, 
for instance, it is required to find the best approximative multinomial with a minimal 
degree. With this in mind, in two examples below we consider some accepted 
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techniques, used for solving the mentioned problem in approximation and/or 
regression analysis theories. 

Example 4. Let 

{x„}= {-l+0.2(/i-l)}; / (x) = sinx; {y„} = { [ kf{x n ) ] Ik }, (14) 

where square brackets mean the integer part; n = 1 ,2, ... , 1 1 ; fix) is a given 
function, used for generating the array {y „ }; a factor k = 10 3 and its presence in (14) 
is necessary for “measuring” all values of_y„ within error e =10 -3 . It is required for 
the presented dependence {y „ , x„} to find the best approximative multinomial with a 
minimal degree. 

A. As well known in approximation theory 22 , if the type of the function /(x) is 
given and either (m+l)-derivative of this function weakly varies on the realisation 
{x„}, or on the x-interval [-1,1] the function /(x) is presented in form of even- 
converging power series, then the problem of finding the best even-approximating 
multinomial for the discrete dependence (y„, x„} offers no difficulty. Indeed, in the 
first case the solution of problem is an interpolative multinomial Pm( B, x) with a set 
of Chebyshev points (this multinomial is close to the best even-approximating one). 
In the second, case one may obtain the solution by the following economical 
procedure of an even-converging power series: 

1 . Choose the initial part of truncated Taylor series, approximating the function 
/(x) within error Zm < e, as the multinomial P m{ B, x) (the multinomial with the 
degree M and vector parameter B); 

2. Replace e M with Em - I b M \ /2 M_1 , where b M is a coefficient of the 
multinomial Pji/(B,x) at x^; 

3. If e m > 0, then replace the multinomial Pat (B, x) with the multinomial 

V,W = P M (x)-(b M 1 2 M ~ l )T M (x ), ( 15 ) 

where Tm(x) is Chebyshev multinomial: T 0 = 1, T = x and when M > 2 
T m = 2xT m _ { - T m _ 2 . Then decrement M by one and go to point 2. If £ w < 0, 
then go to point 4; 

4. End of computations: the multinomial P m (B, x) is the desirable one. 

By means of the foregoing economical procedure one may easy obtain that the 
multinomial with a minimal degree, even-approximating the function sinx, given 
within error e =10 -3 on the x-interval [-1, 1], has the following form 

Pi(x) = (383/384) x- (5/32) x 3 . (16) 

B. Let 1 < M < 9 and in the multinomial P \ (B, x) = 'L% = o‘k m b m x m all X m = 1 . By 
determining LS-estimates of vector parameter B for each value of M on the formed 
above array {y„,x„ }, we find that in all obtained approximative multinomials 
Px(B ', x), as well as in Taylor series of function sinx, the values of coefficients 
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b ' 21 = 0 at / = 0, 1, ... , 4. Thus, regression analysis of the discussed array {y„,x„} 
allows to determine a form of the best approximative multinomial: 

^ n+i opt (B,x) = Z>!X + bpr 3 + ... +b 2l+l x 2l+i + ... . (17) 

We note, if / = 1, then the values of parameters b\ and 63 , computed by regression 
analysis method (LS-method) for model (17), coincide with ones, shown in (16). 

We remind, that in classical variant of regression analysis theory the best 
approximative multinomial is chosen by the minimal value of the testsx = S\J(N-I\), 
where S\ is residual sum-of-squares, h = l N is total number of readings; X m 

is such characteristic number that X m = 0, if the approximative multinomial P\ (B, x) 
does not contain term b„x m , and Xm = 1, otherwise. For approximative models (17) 
and / = 0, 1, 2 , 3 the computation values of 5 are following 

/ 0 1 2 3 

s' 0.033 0.00063 0.00043 0.00054 

Since 5 has the minimal value at / = 2, for the discussed array in the frame of 
classical variant of regression analysis theory, the multinomial 

P 5 (x) = bix + bp? + bp 5 ( 1 8) 

is the best approximative one. 

C. Since, for each n in the data array (14), the difference of \y„ -y„' I must be 
kept within the limit of the error e = 10 -3 , the general solution of the discussed 

6 7 

problem has the following form ’ 

/£, x)=g a {(1.0012 - 0.0001^) x - (0.161200 - 0.000 1 27^) x 3 }, (19) 

where a = 0.001; 0 < ^ < 49 and, consequently, bi = 0.9987 ± 0.0025 and b 2 = - 
0.1582 ±0.0030. 

By analysing solutions (16), (18) and (19) we conclude that in the considered 
case 

the solution of the problem on finding the best fitting multinomial depends on 
the type of the used mathematical theory. 

Example 5. In some software products {for instance, in the different versions of 
software CURVE, designed by AISN} the solutions of problems on finding the best 
approximative models are found by the magnitude of a determination coefficient R, 
which value may be computed by a set of formulae 

*1 = Vi-Gr/G, Qr = X" 1 (y„ - y;> 2 , Q = X" 1 (y n ~ J ) 2 ( 20 ) 
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where y = 'L^ =l y n ! N , y n is /i-th reading of dependent variable, y' n is n-th value of 

dependent variable, computed on the fitting model; or by formula (firstly offered by 
K.Pearson) 

riUO’.-jW-n 

2 V2"=,(^-30 2 X",(:v„-r> 2 <21) 

(evidently, one may easy obtain formula (21) from formula (10)}. 



Table 5. 



The simulative data array to example 5 



n 


Yn 




n 


Yn 


*n 


1 


85 


11 


7 


205 


3.8136396 


2 


105 


5.6002132 


8 


225 


3.5774037 


3 


125 


5.0984022 


9 


245 


3.4193292 


4 


145 


4.7047836 


10 


265 


3.2903451 


5 


165 


4.3936608 


11 


285 


3.1026802 


6 


185 


4.0998636 


- 


- 


- 



There is a mathematical proof 23 of the equivalence of formulae (20) and (21). 
But, if the value of coefficient^ 2 is computed within error > 10 -8 , 

in actual practice, for some data arrays, firstly, R 2 2 ^Ri and, secondly, R 2 > 1. 
For instance, if one fits the simulative data array, presented in table 5, by the 
multinomial /m(B, x) with M= 8, then software CURVE-2.0 will give the value Ri 2 
= 1.00040. 



4 Problems of quantitative processing experimental dependences 

found for heterogeneous objects 

As it follows from the general consideration 24,25 , in practice at analysis of 
experimental dependences found for heterogeneous objects, three various situations 
can be realised: the heterogeneity of investigated objects causes a) no effect; b) a 
removable (local) inadequacy of postulated fitting model; c) an irremovable (global) 
inadequacy of the postulated model. In this section we discuss some computative 
difficulties which may occur at analysis of the mentioned experimental 
dependences. 

Example 6. As we know from Sect. 2.1 if F( A, x) is a linear model {/-"(A, x) = 
'5Li'=i a i h i( x )} ^en the value of A', minimising residual sum-of-squares S, is 
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computed by (1). Let rank H < L, or, in other words, there is a linear dependence 
between columns of matrix H: 

Ci h\ + Ci fi2 + . . . + cjJil = 0, (22) 

where at least one coefficient c/* 0. In this case matrix (H T H) -1 does not exist, that 
means one cannot find A' from (1). Such situation is known as strict 
multicollinearity. 

In the natural science investigations values of the independent variable X are 
always determined with a certain round-off error, although this error may be very 
small. Therefore, if even strict multicollinearity is present, in practice the equation 

(22) is satisfied only approximately and therefore rankH = L. In such situation 
application of equation (1) to find the estimate of vector parameter A gives A' 
values drastically deviating from true coefficients values ’ ’ . 

To correct this situation in regression on characteristics roots 26 it is suggested to 
obtain the information about the grade of matrix H T H conditioning from values of 
its eigennumbers Xj and first elements To y of its eigenvector Vy and to exclude from 
regression such /-components, whose eigennumbers Xj and elements Toy are small. 
Following values are recommended to use as critical ones: Xc r = 0.05 and V CT = 0.1. 
Let us demonstrate, that 

in some practical computations the difference between A' chr and A is of (1) can 
be explained not by the effects of multicollinearity, but by regression model 
inadequacy, which disappears simultaneously with the effects of multicollinearity 
after removing outliers. 

Indeed, let data array be following 

(y„, X„} ={1 + 0.5 n + 0.05 n 2 + 0.005 n 3 ; n, n 2 , « 3 }, (23) 

/t = 1,2, ... , 11 and we introduce two outliers in (23), by means of increasing values 
_y 3 and y% on 0.5. For this data array we obtain the following computation results: 

N step = 1, N= 11: A'ls = (1.017; 0.542; 0.0462; 0.0050), 

« a = {3,8} A / chr = (1.046; 0.516; 0.0516; 0.0047); 

N step = 2, N= 10: A'ls = (0.764; 0.801; -0.0169; 0.0090), 

«a= {3} A' CHR = (0.764; 0.801; -0.0169; 0.0090), 

where N slep is a number of the step in used computative procedure; « a is a vector to 
indicate the numbers of anomalous readings, contained in analysing array on the 
first and second steps of used computative procedure; N is the general quantity of 
analysing readings. In particular, after the first step of computative procedure from 

(23) the reading with number 8 is removed; after the second step — readings with 
numbers 3 and 8. And after the second step the values of A' are restored without any 
distortion by both examined algorithms. 
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By analysing obtained computation results one can conclude that 
the difference between A is and A 'chr may be caused not only by 
multicollinearity but also, for instance, by a set of outliers presented in the data 
array. 

Example 7. In table 6 we adduce an experimental data array, obtained by 
N.P.Bobrysheva , when she investigated magnetic susceptibility (x, relative units) 
of polycrystalline system V^Al^O^ (x = 0.078) on the temperature (T, K). Let us 

consider some computation results of quantitative processing this temperature 
dependence. 



Table 6. 



The experimental dependence of magnetic susceptibility of system 



V^Al^O^ s (x = 0.078) on temperature 



n 


T n 


In 


i 


80 


10.97 


2 


121 


8.06 


3 


144 


6.94 


4 


182 


5.56 


5 


202 


5.11 


6 


214 


4.75 


7 


220 


4.62 


8 


267 


4.00 



n 


T n 


In 


9 


292 


3.79 


10 


351 


3.31 


11 


360 


3.20 


12 


385 


3.06 


13 


401 


3.00 


14 


438 


2.78 


15 


464 


2.67 


16 


486 


2.56 



n 


T n 


X/i 


17 


501 


2.48 


18 


512 


2.46 


19 


523 


2.42 


20 


559 


2.28 


21 


601 


2.12 


22 


651 


2.02 


23 


668 


1.97 


- 


• 


_ 





a) b) 

Figure 4. Experimental (circles) and analytical (continuous curves) plots of dependences 
X(7)(a) and 1/(x + X 2 )~ T’(b) for system v x Al}_xOi 5 (x = 0.078). 



A. In figure 4 (a, b) for the discussed system experimental (circles) and analytical 
(continuous curves) plots of dependences % (7) and l/(x + % 2 ) ~ T are shown. For 
construction of analytical (continuous) curves we use modified Curie - Weiss 
law 6 ’ 7 ’ 24 ’ 25 
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X = Xo+ C/(T + 0), (24) 

where % is the experimental magnitude of specific magnetic susceptibility; T is 
absolute temperature, K; C, 0 and Xo are parameters: C— 988; 0 = 14, K; Xo = 0.54. 
From analysing graphical information, presented in figure 4 (a, b ), one can conclude, 
that 

the magnetic behaviour of system 'V x Al l _ x O l5 (x = 0.078) is well explained by 
the modified Curie - Weiss law (24). 

B. In figure 5(a) the dependence Ax = X - C/(T + 0) - Xo on T for system 
VxAli.xOi^ (x = 0.078) is shown. Since AXnu« = 0.2 » e = 0.01, where e is the 

measurement error in the diccussed experiment, we obtain 

in contradiction with the statement of point (A) in this case modified Curie - 
Weiss law (24) is an inadequate approximative model or, in other words, there is a 
set of outliers in the analysing experimental dependence. 

C. After deleting first 5 readings from the initial data array the parameters values 
of modified Curie - Weiss law (24) have magnitudes C— 1386; 0 = 89, K; Xo = 0.14. 
The plot of dependence Ax = X - C/(T + 0) - Xo with foregoing parameters values is 
shown in figure 5(b). 





a) b) 

Figure 5. Plots of A% = *- CI(T + 0) -Xo on 7 for system v x ai 1 _ x 0 1-5 (x = 0.078). 



Analysing plots A x(T), presented in figure 5(a, b ), and comparing with each 
other the parameters values of equation (24), mentioned in points (A) and (C), we 
conclude, that 

neglect of the local inadequacy of the approximative model in the discussed 
experiment leads to distortion of both form of function Ax(T) and parameters values 
of the modified Curie - Weiss law (24). 

Thus, if, for proving well-fitted properties of equation (24), researchers 27-30 
suggest to look at the graphic representation of dependences x (T) or 1/x (7), for the 
proof completeness one should ask these researchers to present information about 
the measurement error of values x and plots Ax(7) = x - C/(7+0) - X o- 
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For the sake of convenience, the main causes, given rise to computative 
difficulties at analysis of experimental dependences found for heterogeneous 
objects, and methods of their overcoming are adduced in table 7 together. In this 
table all methods, overcoming computative difficulties, are marked by the symbol 
0, if at present they are in the rough or absent in modem data analysis theory. 



Table 7. 

Main causes and overcoming methods of computative difficulties in 
modern data analysis theory 





Main causes 


Methods of overcoming 


1 . 


Impossibility to take heed of preset 
measurement accuracy of dependent 
variable values in the frame of 
accepted data analysis model 


0 Modification of data analysis model 


2. 


Limited accuracy of computations 


Increasing computation accuracy 


3. 


Point estimation of parameters 


Replacing the point estimation of 
parameters by interval one 


4. 


The deficient measurement accuracy 
of dependent variable values 


Increasing measurement accuracy 
of dependent variable values 


5. 


Ill-conditioning of estimation problem 


© Using alternative estimations methods; 
Increasing measurement accuracy of 
dependent variable values; 

© Revealing and removing outliers; 
Designing experiments 


6. 


Presence of outliers in analysing data 


0 Revealing and removing outliers; 




arrays 


0 Robust estimation of parameters 


7. 


Inadequacy of approximative model 


0 Eliminating inadequacy of 
approximative model; 

© Using advanced estimations 
methods 


8. 


Finding only single solution of the 
estimation problems for contaminated 
data array in the frame of modem data 
analysis theories 


© Finding a family of solutions 



Using information presented in table 7, let us clear up a question, whether one is 
able in the frame of modem data analysis theory to obtain reliable solutions for the 
problems of quantitative processing of experimental dependences, found for 
heterogeneous objects. 

Let, when an investigated object is homogeneous, a connection between 
characteristics y and X exist and it be close to functional one: y = F( A, X). As we 
said already in beginning of this section, in the discussed experiments three various 
situations can be realised: the structural heterogeneity 
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1) has no effect on the experimental dependence {y„, X„} or, in other words, in 
this case it is impossible to distinguish the homogeneous objects from 
heterogeneous ones on the dependence {y n , X„}; 

2) leads to a distortion of the dependence {y„, X„} in some small region 
{X nj } c {X„} {the approximative model F( A, X) has removable (local) 

inadequacy}. In this case for extracting effects, connected with the presence of a 
homogeneity in the investigated objects, one may use the following way 6, 7 

i) solve the problem on revealing outliers {y„, , X nj } ; 

ii ) determine the value A' on readings {y„, X„} \ (y ni ,X n) } {we remind, that 
a set of {y„, X„} \ {y B| , X n] } is to be well-fitted by the model F( A, X) } ; 

Hi) detect a type and degree of the effects, connected with the presence of a 
homogeneity in the investigated objects, on the data array {y n) -F(A', X B) ), X n| } . 

It follows from point 6 of table 7, that at solving problem (/) in actual practice 
some difficulties, which are unsurmountable in the frame of modem regression 
analysis theory, can be arisen; 

3) leads to a distortion of the dependence {y„ , X„} in a big region {X n| } £ {X„}: 

{the approximative model F( A, X) has irremovable (global) inadequacy}. 

It follows from point 7 of table 7, that in this case it is impossible to find a 
reliable solution of the discussed problem in the frame of modem data analysis 
theory. 

Summarising mentioned in points (1) - (3), we conclude 

since at present the methods, marked by the symbol 0 in table 7, are not effective 
for overcoming computative difficulties or absent in modem data analysis theory, 
one is not able to obtain reliable solutions for the problems of quantitative 
processing of experimental dependences found for heterogeneous objects. 

From our point of view, one of possible ways, overcoming computative 
difficulties in modem data analysis theory, is further development of this theory by 
means of translation of this theory from Aristotelian “binary logic” into more 
progressive “fuzzy logic” 6 ’ 7 ’ 24 ' 25, 31 ’ 32 . 

References 

1. Dumitrescu C and Seleacu V Some notions and questions in number theory (Erhus 

University Press, Vail, 1995) 

2. Smarandache F Paradoxist Mathematics (Pennsylvania, 1985) 

3. Linnik Y V Least squares method and foundations of mathematical-statistical theory of 

observations processing (Moscow, Publ. H. of phys.-math. lit., 1958) (In Russian) 

4. Rao C P Linear statistical inference and its applications. 2nd ed. (New York, Wiley, 

1973) 

5. Ermakov SM and Zhiglavskiy A A Mathematical theory of optimal experiment 

(Moscow, Nauka, 1987) (In Russian) 

6. Chebrakov Y V The parameters estimation theory in the measuring experiments (S.— 

Petersburg, S.-Petersburg State Univ. Press, 1997) (In Russian) 



79 




7. Chebrakov Y V and Shmagin V V Regression data analysis for physicists and chemists 

(S -Petersburg, S -Petersburg State Univ. Press, 1998) 

8. Huber P J Robust Statistics (New York, Wiley, 1981) 

9. Rousseeuw P J and Leroy A M Robust regression and outlier detection (New York, 

Wiley, 1987) 

10. Andrews D F Technometrics (1974) 16 523 

11. Draper KR and Smith H Applied regression analysis (New York, Wiley, 1981) 

12. Gill P E, Murray W and Wright M H Practical Optimization (New York - London, 
Academic Press, 1981) 

13. Mudrov V I and Kushko V P Methods of processing measurements: Quasi - believable 

estimates (Moscow, Radio i svyaz, 1983) (In Russian) 

14. Adler Y P The introduction to Russian edition of the book: Mosteller F and Tukey J W 
Data analysis and regression (Moscow, Finansy i statistika, 1982) (In Russian) 

15. Aivazyan S A, Yenyukov I S and Meshalkin L D Applied statistics. Analysis of 
dependences (Moscow, Finansy i statistika, 1985) (In Russian) 

16. Bendat J S and Piersol A G Random data: Analysis and measurement procedures 

(New York, Wiley, 1986) 

1 7. Rozanov Y A Probability theory, random processes and mathematical statistics: A 

textbook for higher education. (Moscow, Nauka, 1989) (In Russian) 

18. Afifi A A and Azen S P Stastistical analysis: A computer-oriented approach. 2nd ed. 

(New York - London, Academic Press, 1979) 

19. Lloyd E and Lederman W Handbook of applicable mathematics. Volume IV: Statistics. 

Parts A and B (New York, Wiley, 1984) 

20. Mandel I D Cluster analysis (Moscow, Nauka, 1988) (In Russian) 

21. Seber G A F Multivariate observations (New York, Wiley, 1984) 

22. Laurent P J Approximation et optimisation (Paris, Hermann, 1972) 

23. Vuchkov I, Boyadzieva L and Solakov E Applied linear regression analysis (Moscow, 
Finansy i statistika, 1981) (In Russian) 

24. Chebrakov Y V and Shmagin VVJ. Izv. Vysh. Uch. Zav. Fizika (1995) 5 40; (1998) 6 

1 10 (Engl. Transl. Russian Phys. J. (1995) 5 470; (1998) 6) 

25. Chebrakov Y V J. Izv. Vysh. Uch. Zav. Fizika (1997) 7 103 (Engl. Transl. Russian 
Phys. J. (1997) 7 687) 

26. Webster J, Gunst R and Mason R Technometrics (1974) 16 513 

27. Bobrysheva N P Magnetic dilution of a series of oxides containing elements of the first 

transitional period : Condensed version of the PhD thesis (Leningrad, 1974) (In 
Russian) 

28. Amine K, Tressand A, Imoto H et al. J. of Solid State Chemistry (1992) 96 287 

29. Miyata N, Kimishima Y, Akutsu N and Oguro I J. Magnetism and Magnetic Materials 

(1990) 90&91 337 

30. Zayachuk D M, Ivanchuk R D, Kempnik V I and Mikitjuk V I Fizika tverdogo tela 
(1996) 36 N° 8 2502 (In Russian) 

31. Chebrakov YV and Shmagin VV Correct solutions of fit problems in different 
experimental situations In: Advanced Mathematical Tools in Metrology III 
(Singapore, World Scientific, 1997)246 

32. Chebrakov Y V and Shmagin V V Some nontraditional approaches to solving a few 
data analysis problems In: Integral methods in science and engineering (Edinburg, 
Longman, 1997) 2 69. 



80 




THE NORMAL BEHAVIOR OF 
THE SMARANDACHE FUNCTION 



Kevin Ford 



Let S(n ) be the smelliest integer k so that n\k\. This is known as the Smarandache 
function and has been studied by many authors. If P(n) denotes the largest prime 
factor of n, it is clear that S(n) ^ P{n). In fact, S(n ) = P{n) for most n, as noted 
by Er dos [E]. This means that the number, N(x), of n ^ x for which S(n) ^ P(n) 
is o(x). In this note we prove an asymptotic formula for N(x). 

First, denote by p(u) the Dickman function, defined by 

p(u) = l (O^u^l), p(u ) = 1 - [ — ~ - dv (u> 1). 

J 1 v 



For u > 1 let £ = £(u) be defined by 



u = 



— 1 



€ 



It can be easily shown that 

£(u) = log u + log 2 u + O » 

where log fc x denotes the kth iterate of the logarithm function. Finally, let u 0 = 
uq(x) be defined by the equation 



log x = ul£{u 0 ). 



The function uq(x) may also be defined directly by 

logx -- uq (x 1/u ° — 1^ . 
It is straightforward to show that 



( 1 ) 



u o 



2 log x ' 
. log 2 X 



log 3 x + log 2 + Q 



2 log 2 x 2 log 2 x 



loggx' 



We can now state our main result. 
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Theorem 1. We have 



N ( x ) ~ ^ ( ^i° S2) (log x log 2 *) 3/ V 1/u V(wo)- 

There is no way to write the asymptotic formula in terms of “simple” functions, 
but we can get a rough approximation. 

Corollary 2 . We have 



N(x ) = xexp j— (\/2 4- o(l))v/logxlog 2 x| . 

The asymptotic formula can be made a bit simpler, without reference to the 
function p as follows. 

Corollary 3. We have 



N(x) ~ --- 1 2 ^ > ° - (log x) * (log 2 x)x* 





5 



where 7 = 0.5772 ... is the Euler-Mascheroni constant. 

This will follow from Theorem 1 using the formula in Lemma 2 which relates 
p(u ) and £(u). 

The distribution of S(n ) is very closely related to the distribution of the func- 
tion P(n). We begin with some standard estimates of the function \&(x,y), which 
denotes the number of integers n ^ x with P(n ) ^ y. 

Lemma 1 [HT, Theorem 1 . 1 ]. For every e > 0, 



^{x,y) = xp(u) 



(l+O 



/ log(u + l) \\ 

V lo g y ) ) ’ 



u = 



logx 

logy’ 



uniformly in 1 ^ u ^ exp{(logy ) 3 / 5 € }. 
Lemma 2 [HT, Theorem 2.1]. For u ^ 1, 




Lemma 3 [HT, Corollary 2.4]. If u > 2 , |u| ^ u/2, then 

p(u — v) = p(u) exp {v£(u) + 0((1 + v 2 )/u)}. 
Further, if u > 1 and 0 ^ v ^ u then 

p(u — v) p{u)e v ^ u \ 
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We will show that most of the numbers counted in N(x ) have 



P(n) « exp |^ilogxlog 2 x| 



Let 

Yi = exp > y 2 = Yi = exp {2 x /logxlog 2 | . 

Let N\ be the number of n counted by N(x) with P(n) ^ Yj, let N 2 be the number 
of n with P(n ) ^ Y 2 , and let N 3 = N(x) - Ni — N 2 . By Lemmas 1 and 2, 

Ni ^ ^(XjYi) = xexp{-(1.5 + o(l))v / logxlog 2 x}. 

For the remaining n ^ x counted by N(x), let p = Pin). Then either p 2 |n or for 
some prime q < p and 6 ^ 2 we have q b || n, q b \ p!. Since p! is divisible by 
and b ^ 21ogx, it follows that q > p/(31ogx) > p 1/2 . In all cases n is divisible by 
the square of a prime ^ Y 2 f (3 log x) and therefore 

v— ^ x 6xlogx f , „ /: 1 

iV 2 ^ ^ xexp |-1.9vlogxlog 2 x j . 

^ r 2 ” 2 

3 log x 



Since q > p 1 ^ 2 it follows that <p p / 9 l || pi. If n is counted by N 3 , there is a number 
6^2 and prime q G [p/ 6 , p] so that g 6 |n. For each 6 ^ 2, let iV^j/x) be the number 
of n counted in A 3 such that || n for some prime q ^ p/ 6 . We have 

^iV 3)fc <x r 31 ^ gx N j < xexp {—(5/3 + o(l))\/logx log 2 x| . 

6^6 V 1 / 



Next, using Lemma 1 and the fact that p is decreasing, for 3 ^ 6 ^ 5 we have 




By partial summation, the Prime Number Theorem, Lemma 2 and some algebra, 

IV 3, 6 < exp {-(1.5 + o(l))v / logxlog 2 x| . 
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The bulk of the contribution to N(x) will come from JV^. Using Lemma 1 we 
obtain 



( 2 ) 




By Lemma 3, we can write 



The contribution in (2) from p near Yf or Y 2 is negligible by previous analysis, and 
for fixed q G [Yi, Y 2 /2] the Prime Number Theorem implies 



E 

q<p<2q 



1 

p 



log 2 
log? 



+ 0((log?) 



2 ) 



log 2 
logp 




Reversing the roles of p, q in the second sum in (2), we obtain 




By partial summation, the Prime Number Theorem with error term, and the change 
of variable u = log x/ log p, 



(3) 1V 3; 2 = 1 + O 



(v^)Hr( 



p{u — 2) log 2 



u 



+ p(u — 3) ) x 1 / U du, 

log 2 # 



where 




/ 



logx 

log 2 X ’ 



U 2 = 6ui- 



The integrand attains its maximum value near u = uo and we next show that the 
most of the contribution of the integral comes from u close to uq. Let 



“o „ , — /l°g3z\ 1/2 

“ 0 = Ioo’ = ^ = - 

where A" is a large absolute constant. Let 7i be the contribution to the integral in 
(3) with |u — «o| > U/’o, let / 2 be the contribution from w\ < |u — uo| ^ u»o, let / 3 
be the contribution from ty 2 < |u — uol ^ ^l) and let I 4 be the contribution from 
Jxr — xi 0 | ^ W 2 - First, by Lemma 2, the integrand in (3) is 



exp {“(^“^ +0(1) ) Vlogxlog^}, 
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The function 1/c + c/2 has a minimum of \/2 at c = y/2, so it follows that 

h < exp (\/2 + 10 5 ) v /logilog 2 x| . 

Let u = uq — v. For w\ ^ |v| ^ w 0 , Lemma 2 and the definition (1) of uo imply 
that the integrand in (3) is 



^ p(uq) exp <{ t/f(u 0 ) - 



v v 2 v 3 \ ( v 2 

1 H 1 2 d 3 ) + O ( h log Uq 

Uq U 5 UlJ \Uq 



p(u 0 )x 1/Auo exp 



-O log X + O (— + log Uq 



< g; p(uo)x 1/,u ° exp < —O.9—3 logx 



for K large enough. It follows that 



l 2 <^u 0 p(uo)x 1 /“° exp{— 201og 2 x} <§: (logx) 10 p(u 0 )x 1//u °. 

For the remaining u, we first apply Lemma 3 with v = 2 and v = 3 to obtain 



h+h = ( 1 + 0 



p{u)x 



l/uf e2i(u) , lo g2 3£(u 



u logx 



3 ^ (u) du 



We will show that /3 r A » p(uo)x 1 / u °(logx) 3 / 2 , which imphes 



(4) N(x) =(l + 0 (vW)) /J*”' ^ 



u log X J 



This provides an asymptotic formula for N(x), but we can simplify the expression 
somewhat at the expense of weakening the error term. First, we use the formula 



£(u) = log u + log 2 u + O 



log2_u\ 
log u )' 



and then use u = uq + 0(uq ) and (1) to obtain 



^3 + I A — ( 1+0 



log 2 X J J 4 



[1 + log 2)x(log x) 2 (log 2 x)2 



r*tl 0 + Wi 



p(u)x 1//u du. 



By Lemma 3, when tu 2 ^ |u| ^ w\, where u = uq — v, we have 



p(uq — v)x “o » <g; p(u 0 )x "o exp < v£(u 0 ) 



logx ( V v 2 V 3 



1 1 2 ■* 3 

Uq \U 0 Ufi Ug 



<C p(uq)x u o exp < 3 log x 



p(uq)x “o exp < j logx 

l U Q 

< p{u 0 )x =? (log 2 x) 3 
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provided K is large enough. This gives 



/ p(u)x 1 / u du <£. p(uo)x 1/u °(loga;) 1/4 (log 2 :r) 3 ' 5 . 

JW2^\U lio|^T^l 

For the remaining v, Lemma 3 gives 

p(uo-v)x 1/(uo v) = (l + o(j^))p(«o)z 1/7X0 exp log xj. 

Therefore, 

pOo) 1 X ^ f p(u)x 1/u du= (l + offef)) f exp dv. 

Juo W2 ' ' ' ' J W2 V. U 0 J 

The extension of the limits of integration to (— oo, oo) introduces another factor 
1 + 0((log 2 x) 1 ), so we obtain 



J 3 + U = (l + O )) ^ {1 ^ 2) {logxlog, xf /A p{u 0 )x "o 

and Theorem 1 follows. Corollary 2 follows immediately from Theorem 1 and (1). 
To obtain Corollary 3, we first observe that £'(u) ~ u 1 and next use Lemma 2 to 
write 

e 7 f r uo 

By the definitions of £ and uq we then obtain 

r u o /•£(“ o) p v — 1 

/ £(t) dt = I e v dv 

J l Jo v 

= / £ J dv 

Jo v 

= logx_ r'-^t e v - 1 dv 

Uo Jo V 




Corollary 3 now follows from (1). 
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On values of arithmetical functions at factorials I 



J. Sandor 



Babe§-Bolyai University, 3400 Cluj-Napoca, Romania 



1. The Smarandache function is a characterization of factorials, since S(k\) = k, and 
is connected to values of other arithmetical functions at factorials. Indeed, the equation 



S(x) = k (k > 1 given) 



( 1 ) 



has d(k\) - d((k - 1)!) solutions, where d(n) denotes the number of divisors of n. This 
follows from {x : S(x) = k} = { x : z|fc!, x \ (k - 1)!}. Thus, equation (1) always has 
at least a solution, if d(k\) > d((k - 1)!) for k > 2. In what follows, we shall prove this 
inequality, and in fact we will consider the arithmetical functions ip, <7, d,u, 0 at factorials. 
Here <p(n) = Euler’s arithmetical function, cr(n) = sum of divisors of n, u/(n) = number 
of distinct prime factors of n. fl(n) = number of total divisors of n. As it is well known, 

r 

we have <^(1) — c/(l) = 1, while u;(l) = fi(l) = 0, and for 1 < JJp- 1 (a : > 1, p x distinct 
primes) one has 

r pV +1 - 1 



i=i 



Pi - i 



r (^) = n 

t — i 

u{n) = r, 

r 

C(n) = ^ a,, 

1=1 
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(2) 



^( n ) = m* + 

i=\ 

The functions p>^cr^d are multiplicative, lj is additive, while Cl is totally additive, i.e. 
c p,cr,d satisfy the functional equation f(mn) = f° r ( m ? n ) — ^ while u, Cl satisfy 

the equation g{rnn) = g(m) T g(n) for (m,n) = 1 in case of cj, and for all m,n is case of 
Cl (see [1]). 

r r 

2. Let m = 11?^% « = II P?' (a,-,/?,- > 0) be the canonical factorizations of m and n. 

t=i *=i 

(Here some a,- or /?,- can take the values 0, too). Then 

d(mn) — XI(o.‘ + Pi + 1) > H(A + 1) 

;=i «'=i 

with equality only if a,- = 0 for all i. Thus: 

d(mn) > d(n) (3) 

for all m, n, with equality only for m = 1. 

r r r 

Since n^+ft+ijsrKo.-njnw -f 1), we get the relation 

: = 1 t=l i=l 

J(mn) < d(m)d(n) (4) 

with equality only for (n,m) = 1. 

Let now m = k, n = (k — 1)! for k > 2. Then relation (3) gives 

d{kl) > d((k - 1)!) for all k > 2, (5) 

thus proving the assertion that equation (1) always has at least a solution (for k = 1 one 
can take x = 1). 

With the same substitutions, relation (4) yields 

d(k\) < d((k - 1 )l)d{k) for k > 2 (6) 
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Let k = p (prime) in (6). Since ((p - l)!,p) = 1, we have equality in (6): 



d{{p- 1)!) 



= 2 , p prime. 



3. Since S(k\)/k\ -4 0, 1 ^ k 00 ’ ° ne may ask the similar 

problems for such limits for other arithmetical functions. 



It is well known that 



cr(n!) 



— y oo as 7i — y oc. 



In fact, this follows from a(k) = d = , so 



d|fc d|A: 






as it is known. 



From the known inequality ([!]) <p(n)cr{n) < n 2 it follows 



” > <T ( W ) 

cp(n) — n 



oo, implying 



¥>(«!) 



-4 0 as n 4 oo. 



Since c p(n) > d(n ) for n > 30 (see [2]), we have <p(n!) > d(n\) for n! > 30 (i.e. n > 5), 



so, by (9) 



0 as n — > oc. 



In fact, much stronger relation is true, since > 0 for each 5 > 0 (n — > oc) (see 

[1]). From and the above remark on cr(n!) > nllogn. it follows that 

n! n! 



d(n\) 

lim sup — — log n < 1 . 

n— too n! 
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These relations are obtained by very elementary arguments. From the inequality 
<p{n)(u(n) + 1) > n (see [2]) we get 



u;(n!) — y oo as n — y co 



( 12 ) 



and, since fl(s) > u;(s), we have 



Q(n!) — >• oo as n — > oo. 



From the inequality nd(n ) > ip(n ) + cr(n) (see [2]), and (8), (9) we have 



d(n!) -> oo as n -4 oo. 



(13) 



(14) 



This follows also from the known inequality ip(n)d(n) > n and (9), by replacing n with 
n!. From cr(mn) > mcr(n) (see [3]) with n = (k — 1 )!, m = k we get 

<r(Jfe!) 



<r((fc - 1)!) - k (A: - 2) 
and, since cr(mn) < a(m)a(n), by the same argument 



(15) 



< „(*) (fc > 2). 



(16) 



Clearly, relation (15) implies 

ajkl) _ 
k™a{{k-l)\) 

From < jj(ran) < m<p(rz), we get, by the above remarks, that 

l o(k\) 

{k ~ 2) 

implying, by tp(k) -4- oo as A: — 4 oo (e.g. from ip(k) > \Jk for k > 6) that 

hm -,^.. =+ oo. 



(17) 



(18) 



k-+co <p((k — 1 )! 



(19) 
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a(kl) 

By writing cr(Ar!) — a((k — 1)!) = a((k — 1)!) — — — — — 1 , from (17) and 

.cr((K - 1)!) 

cr((k — 1)!) — >• oo as k -> oo, we trivially have: 



lim [er( A:!) — a((k — 1)!)] = +oc. 

k—HDO 



In completely analogous way, we can write: 



lim [<p(A!) — ¥?((£ — 1)')] = +oo. (21) 

k->oo 

4. Let us remark that for k = p (prime), clearly ((A — 1)1, A) = 1, while for k — 
composite, all prime factors of k are also prime factors of ( k — 1)!. Thus 



'(*0 = 



;((A — 1 )! A:) = u((k — 1)!) + cj(A) if k is prime 



'((*- 1)0 



if k is composite (k > 2). 



Thus 



;(&!) — u((k — 1)!) = 



1, for k — prime 
0, for k = composite 



Thus we have 



lim sup[u>(fc!) — uj{{k — 1)!)] = 1 

k—¥oo 

lim inf[u;(A!) — u;((A — 1)!)] = 0 

k-t oo 

Let p n be the nth prime number. From (22) we get 



"((*- 1)0 



n - 1 



if k = p n 



0, if A: = composite. 



Thus, we get 



lim = 1. 

k-t-co cj((A — 1)!) 



The function fl is totally additive, so 



fi(A!) = n((Jfc - 1 )!A) = ft((A - 1 )!) + ft(A), 
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giving 

«(*!)-«((*- !)!) = **(*)• 



This implies 



limsup[n(fc!) — Q,((k — 1)!)] — +oc 

k-t oo 

(take e.g. k = 2 m and let m — > oo), and 



lim inf[H(A:!) — — 1)!)] — 2 

k—*oo 



(take k = prime). 

For fi(fc!)/fi((A: - 1)!) we must evaluate 

fi(Jfc) fl(fc) 



fi((Jfe-l)!) fl(l) + fi(2) + ... + ft(fc-l) 

Since f l(k) < and by the theorem of Hardy and Ramanuj< 

log 2 



X fi ( n ) ~ x log log x (x — »• oo) 



n<r 



so, since 



log k 



{k - 1 ) log log(fc 1) 



— > 0 as k — > oo, we obtain 






oo n((jb - i)i) 



5. Inequality (18) applied for k = p (prime) implies 

1 y(p!) _ 

A-P>((P-1)!) ' 

This follows by ^(p) = p — 1- On the other hand, let A; > 4 
it is known (see [1]) that A:|(fc — 1)!- So <,?(&!) = ^({k — 1)!&) : 
<p(mn) = rrup(n) if m|n. In view of (28), we can write 

1 y(t!) 

9((fc — 1)!) 



(25) 

(26) 



a (see [1]) we have 

(27) 

( 28 ) 

be composite. Then. 
= k?((k — 1)!), since 

(29) 
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cr(p!) 

For the function a, by (15) and (16), we have for k = p (prime) that p < — — — < 

a (\P — 1)!) 



cr(p) = p+ 1, yielding 

1 <r{pl) _ . 

p ' <j((p - 1)!) 

In fact, in view of (15) this implies that 

. 1 ajkl) 

k ' a{{k- 1)!) 

By (6) and (7) we easily obtain 

d(k<) 

d(k)d((k - 1)!) ' 



(30) 



(31) 



(32) 



In fact, inequality (6) can be improved, if we remark that for k = p (prime) we have 
d(k\) = d((k — 1)!) • 2, while for k = composite, k > 4, it is known that k\(k — 1)!. We 
apply the following 

Lemma. 7/n|m, then 

d{mn) < d(n 2 ) ^ 



d(m) d(n) 



Proof. Let m = rb-n/.^n?” ( a ' < a) be the prime factorizations of m and 
n, where n|m. Then 

d(mn ) II(« + a ' + 1) IIC 5 + 1) _ yr f Q + a ' + l\ 

d{m) [](a + l)JJ(/5 + l) V a + l ) 

Now — — — - < — — A ^ o' < a as an easv calculations verifies. This immedi- 

a + 1 a' -r 1 

ately implies relation (33). 

By selecting now n = k, m = (k — 1)!, k > 4 composite we can deduce from (33): 



d(k\) < d(k 2 ) 



(34) 



d{{k - 1)!) - d(k) 

Bv (4) we can write d(k 2 ) < ( d(k )) 2 , so (34) represents indeed, a refinement of relation 



( 6 ). 
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THE AVERAGE VALUE OF THE SMARANDACHE FUNCTION 



Steven R. Finch 
MathSoft Inc. 

101 Main Street 
Cambridge, MA, USA 02142 
sfinch@mathsoft. com 



Given a positive integer n, let P(n) denote the largest prime factor of n and Sfn) denote 
the smallest integer m such that n divides ml 

The function S(n) is known as the Smarandache function and has been intensively studied 
[1]. Its behavior is quite erratic [2] and thus all we can reasonably hope for is a statistical 
approximation of its growth, e.g., an average. It appears that the sample mean E(S) 
satisfies [3] 



E(S(N)) = ^-f J S(n) = c{ 



N 



N ~ Un (N)J 

as A r approaches infinity, but I don't know of a rigorous proof. A natural question is if 
some other sense of average might be more amenable to analysis. 

Erdos [4,5] pointed out that P(n) = S(n) for almost all n, meaning 
|{n< V:P(«)<S(n)}| 



lim 

A— > oo 



N 



= 0 that is, { n < N: P(n ) < S(«)} = o{N) 



as N approaches infinity. Kastanas [5] proved this to be true, hence the following 
argument is valid. On one hand. 



/. = lim E 



r ln (/>(«)) N 



< lim E 






1 x- 

= 



InfS(nl) 



^ i / x — 11111 ^ , - — imi — * / 

"-* 30 V ln(^) ) v ln(«) ) n->od N ^ In (n) 

The above summation, on the other hand, breaks into two parts: 



r 



lim— -I T — — ■ " + 2 — s v " 1 

X^ccN U(«)=S( n ) ln («) P(nHS(n) ln («) J 



ln( P(n)) 



I 



ln(S(n))] 



The second part vanishes: 
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lim 77' 

A’-»oo A 



^ ln(5(/7)) 

\P(n)<S(n) ln(n) 



< lim — • 

A — >0C A 



Zi 

\P(n)<S(n)J 



«• n 

= lim — — = 0 

A —KC A 



while the first part is bounded from above: 



lim — 

N — >oc A 



^ ln(P(»)) 

\P(n)=S(n) l n ( W ) 



SHm > f;M£^ =lima 



A : — >oc A 



/r=l 



ln(«) 



l ln(«) j 



= A 



We deduce that 



lim E 

n — ►oo 



' ln^(n)) " 
ln(n) J 



= /i 



0.6243299885... 



where 2 is the famous Golomb-Dickman constant [6-9]. Therefore A • n is the asymptotic 
average number of digits in the output of S at an n-digit input, that is, 62.43% of the 
original number of digits. As far as I know, this result about the Smarandache function 
has not been published before. 

A closely related unsolved problem concerns estimating the variance of S'. 
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On the Irrationality of Certain Constants Related 
to the Smarandache Function 



J. Sandor 



Babe§-Bolyai University, 3400 Cluj-Napoca, Romania 



1. Let S(n) be the Smarandache function. Recently I. Cojocaru and S. Cojocaru [2] 

S(n) 



have proved the irrationality of ^ 



71=1 



n: 



The author of this note [5] showed that this is a consequence of an old irrationality 
criteria (which will be used here once again), and proved a result implying the irrationality 

offi-ir-sM 

71=1 U * 

°° SYfc) / 5 l\ 

E. Burton 111 has studied series of type Y — — , which has a value € ( e He 

flfi ( k + 1) ! V 2 2/ 

°° S(k ) 

showed that the series Y v - is convergent for all r € N. I. Cojocaru and S. Cojocaru 

k=2 ( k + r ) ! 

00 1 

[3] have introduced the "third constant of Smarandache” namely ^ . cTT’ 

n =2 j 

71 97 

which has a value between and . Our aim in the following is to prove that the 

100 100 

constants introduced by Burton and Cojocaru- Cojocaru are all irrational. 

2. The first result is in fact a refinement of an old irraionality criteria (see [4] p.5): 
Theorem 1 . Let ( x n ) be a sequence of nonnegative integers having the properties: 

( 1 ) there exists no E N’ such that x n <n for all n > n 0 ; 

(2) x n < n — 1 for infinitely many n; 

(3) x m > 0 for an infinity of m. 
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Then the series 



T —f is irrational. 
£l n\ 



Let now x n = S(n — 1). Then 



” S(k) ” x w 

k' 2 ( k + 1 ) 1 n^3 nl ' 



Here S(n — 1) < n — 1 < n for all n > 2; S(m — 1) < m — 2 for m > 3 composite, 



since by 5(m — 1) < -(m — 1) < m — 2 for m > 4 this holds true. (For the inequality 

3 

2 

£(&) < — for k > 3 composite, see [6]). Finally, S(m — 1) > 0 for all m > 1. This proves 
00 S(k) 

the irrationality of ^ . 

k = 2 \ k + 

Analogously, write 



~ S(k) ~ 5(m - r) 

i hi ( k + r ) ! _ m=r+2 m! 

Put = S(m — r). Here 5(m — r) < m — r < m, S(m — r)<m — r<m — 1 for 
r > 2, and 5(m — r) > 0 for m > r + 2. Thus, the above series is irrational for r > 2, too. 

3. The third constant of Smarandache will be studied with the following irrationality 
criterion (see [4], p.8): 

Theorem 2. Let (a n ), ( b n ) be two sequences of nonnegative integers satisfying the 
following conditions: 

(1) a n > 0 for an infinity of n; 

(2) b n > 2, 0 < a n < b n — 1 for all n > 1; 

(3) there exists an increasing sequence (i n ) of positive integers such that 



lim b ln = +oc, lim a, •„/&,„ = 0. 

n— *oc n— kx> 

oo 

TAen £Ae series is irrational. 

n=1 * * * b n 

OO 1 

Corollary. For b n > 2, (b n positive integers), ( b n ) unbounded the series V"] 

bib 2 ...b n 

is irrational. 
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Proof. Let a n = 1. Since limsup6 n = +oo, there exists a sequence (z„) such that 

n-*oo 

6, n oo. Then > 0, and the three conditions of Theorem 2 are verified. 

O.n 

By selecting b n = S(n ), we have b p = S(p) = p -4 oo for p a prime, so by the above 

S(l)S(2)...S(n) iS irrati0na1 ' 
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Smarandache Magic Squares 
Sabin Tabirca* 

•Bucks University College, Computing Department 



The objective of this article is to investigate the existence of magic squares made with 
Smarandache’s numbers [Tabirca, 1998], Magic squares have been studied intensively and 
many aspects concerning them have been found. Many interesting things about magic 
squares can be found at the following WEB page 
http://www.pse.che.tohoko.ac.jp/~msuzuki/MagicSquares.html. 



Definition 1. A Smarandache magic square is a square matrix a<=M n (N ) with the 
following properties: 



a) K, 1 i, j = 1, «}= {?(/) | i = 1,» 2 } 


(i) 


b) (vj = l,n)^a..=i 

i=l 


(2) 


o 

, 

II 

II 


(3) 



Therefore, a Smarandache magic square is a square matrix by order n that contains only 
the elements iS(l),£(2),...,£(« 2 ) [Smarandache, 1980] and satisfies the sum properties (2- 
3). According to these properties, the sum of elements on each row or column should be 
equal to the same number k. Obviously, this number satisfies the following equation 

k=^ . 

n 

n 2 

Theorem 1. If the equation n \ ^ S(i) does not hold, then there is not a Smarandache 

1=1 

magic square by order n. 

Proof 
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2 



!•?(/) 

This proof is obvious by using the simple remark k = — e N . If a &M n (N) is a 

n 

n 2 

Smarandache magic square, then the equation n \ ^S(i) should hold. Therefore, if this 

i=i 

equation does not hold, there is no a Smarandache magic square. 

* 

Theorem 1 provides simple criteria to find the non-existence of Smarandache magic 

n 2 

square. All the numbers l<n<101 that do not satisfy the equation n | y^ S(i) can be 

i=i 

found by using a simple computation [Ibstedt, 1997], They are {2, 3, ...., 100}\{6, 7, 9, 
58, 69}. Clearly, a Smarandache magic square does not exist for this numbers. If n is one 

* 2 

of the numbers 6, 7, 9, 58, 69 then the equation n | S(i ) holds [see Table 1], This does 

i=l 

not mean necessarily that there is a Smarandache magic square. In this case, a 
Smarandache magic square is found using other techniques such us detailed analysis or 
exhaustive computation. 



n 


S(n) 


h(o 

i=i 


6 


3 


330 


7 


7 


602 


9 


6 


1413 


58 


29 


1310162 


69 


23 


2506080 



Table 1. The values of n that satisfy n \ ^ S{i ) . 
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An algorithm to find a Smarandache magic square is proposed in the following. This 
algorithm uses Backtracking strategy to complete the matrix a that satisfies (1-3). The 
going trough matrix is done line by line from the left-hand side to the right-hand side. 



The algorithm computes: 

• Go trough the matrix 

♦ Find an unused element of the set S(l),S(2),...,S(n 2 ) for the current cell. 

♦ If there is no an unused element, then compute a step back. 

♦ If there is an used element, then 

■ Put this element in the current cell. 

■ Check the backtracking conditions. 

■ If they are verified and the matrix is full, then a Smarandache magic 
square has been found. 

■ If they are verified and the matrix is not full, then compute a step 
forward. 



procedure S mar_magi c_square( n) ; 


begin 


procedure forward(col, row); 


coI:=l; row:=l;a[col, row]:=0; 


begin 


while row>0 do begin 


col:=col+l; 


while a[col, row]<n*n do begin 


if col=n+l then begin 


a[col,row] :=a[col,row]+l ; 


col:=l;row:=row+l; 


call check(col,row,n,a,cont); 


end; 


if cont=0 then exit; 


end; 


end 


if cont =0 then call back(coI,row); 


procedure write_square(n,a); 


if cont=l and col=n and row=n 


begin 


then call write_square(n,a) 


for i:=l to n do begin 


else call forward(col,row); 


for j:=l to n do write (S(a[i j]),’ ‘); 


end; 


writeln; 


write(‘result negative’); 


end; 


end; 


stop; 

end; 


procedure back(col, row); 


begin 


procedure check(col,row,n,,k,a,cont); 


col:=col-l; 


begin 


if col=0 then begin 


cont:=l; sum:=0; 


col:=n;row:=row-l; 


for i:=l to col do sum:=sum+S(a[ij]); 


end; 


if (sum>k) or (col=n and sumok) then begin 


end; 


cont:=0; 
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return; 



cont:=0; 

return; 



end; 

sum:=0 end; 

for j:=l to row do sum:=sum+S(a[ij]); end; 

if (sum>k) or (row=n and sumok) then begin 

Figure 1 . Detailed algorithm for Smarandache magic squares. 



The backtracking conditions are the following: 



(v/=l,w)fX; and ( V -/ = 1 >«)lX> =k 
1=1 *=1 


(4a) 


(v/ = l,«)fV < k and (v/ = l,n)2X y . =k. 

7=1 7=1 


(4b) 


(V(/, j) < (row, col))a t ] * a rowcol 


(5) 



These conditions are checked by the procedure check. A detailed algorithm is presented in 
a pseudo-cod description in Figure 1. 

Theorem 2. If there is a Smarandache magic square by order n, then the procedure 
Smar_magic_square finds it. 

Proof 

This theorem establishes the correctness property of the procedure Smar_magic_square. 
The Backtracking conditions are computed correctly by the procedure check that verifies 
if the equations (4-5) hold. The correctness this algorithm is given by the correctness of 
the Backtracking strategy. Therefore, this procedure finds a Smarandache magic square. 

• 

Theorem 3. The complexity of the procedure complexity Smar_magic_square is 

0(n 2 " 2+l ). 

Proof 

The complexity of the procedure Smar_magic_square is studied in the worst case when 
there is not a Smarandache magic square. In this case, this procedure computes all the 
checking operations for the Backtracking strategy. Therefore, all the values 
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S(l),S(2),...,S(n 2 ) are gone through for each cell. For each value put into a cell, at most 
0(n) operations are computed by the procedure check. Therefore, the complexity is 

o(n-(n 2 Y^j = o((i 2nKl ). * 

Remark 1 . The complexity 0(n 2n +I ) is not polynomial. Moreover, because this is a very 
big complexity, the procedure Smar_magic_square can be applied only for small values of 
n. For example, this procedure computes at most 6 73 > 10 56 operations in the case n=6. 
The above procedure has been translated into a C program that has been run on a Pentium 
MMX 233 machine. The execution of this program has taken more than 4 hours for n=6. 
Unfortunately, there is not a Smarandache magic square for this value of n. The result of 
computation for n=l has not been provided by computer after a twelve hours execution. 
This reflects the huge number of operations that should be computed (7" > 10 83 ). 

According to these negative results, we believe that Smarandache magic squares do not 

n 1 

exist. If n is a big number that satisfy the equation n \ ^S(i), then we have many 

i=i 

possibilities to change, to permute and to arrange the numbers 5(1), S(2),. ..,S(n 2 ) into a 
square matrix. In spite of that Equations (2-3) cannot be satisfied. Therefore, we may 
conjecture the following: “ There are not Smarandache magic squares In order to 
confirm or infirm this conjecture, we need more powerful method than the above 
computation. Anyway, the computation looks for a particular solution, therefore it does 
not solve the problem. 

References 

Ibstedt, H. (1997) Surfing on the Ocean of Numbers - a few Smarandache Notions and 
Similar Topics, Erhus University Press. 

Smarandache, F (1980) A Function in Number Theory, Analele Univ. Timisoara, XVJLU. 
Tabirca, S (1998) Problem 24 - Smarandache magic squares, Smarandache Notion 
Journal, Vol 9 - nr 1-2. 



104 




Some inequalities concerning Smarandache’s function 
Sabin Tabirca* Tatiana Tabirca** 

* Bucks University College, Computing Department 
**Transilvania University of Brasov. Computer .Science Department 



The objectives of this article are to study the sum ^5(c/) and to find some upper 

d\n 

bounds for Smarandache’s function. This sum is proved to satisfy the inequality 
^S(d) < n at most all the composite numbers. Using this inequality, some new 

d\n 

upper bounds for Smarandache’s function are found. These bounds improve the well- 
known inequality S{n) < n . 

1. Introduction 

The object that is researched is Smarandache’s function. This function was 



introduced by Smarandache [1980] as follows: 

S:N* — » N defined by S(n) = min {A: e N\ k! =Mn} (Vn eA' *). 0) 

The following main properties are satisfied by 5 : 

(Va,Z> e N *) ( a,b ) = 1 => S(a-b) = max{S(a), S(A)} . (2) 

(Vcr e iV *) S(a ) < a and S(a) = a iif a is prim . (3) 

(\/p e A”*, p prime)(VA: e N *) S(p k ) < p-k . (4) 



Smarandache’s function has been researched for more than 20 years, and many 
properties have been found. Inequalities concerning the function S have a central 
place and many articles have been published [Smarandache, 1980], [Cojocaru, 1997], 
[Tabirca, 1997], [Tabirca, 1988], Two important directions can be identified among 
these inequalities. First direction and the most important is represented by the 
inequalities concerning directly the function S' such as upper and lower bounds. The 
second direction is given by the inequalities involving sums or products with the 
function S’. 
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2. About the sum ^5(J) 

d\n 



The aim of this section is to study the sum ^ S(d). 

d\n 

Let SS(n) = ^S(d) denote the above sum. Obviously, this sum satisfies 

dn 

SS(n) = 'Y J S{d). Table 1 presents the values of Sfn) and SS(n) for n < 50 [Ibstedt, 

\*dn 

1997], From this table, it can be seen that the inequality' SS(n) <n + 2 holds for all 

n= 1, 2, 50 and 12. Moreover, if n is a prim number, then the inequality becomes 
equality' SS(n ) = n . 

Remarks 1. 

a) If n is a prime number, then SS{n) = 5(1) + S(n) = n . 

b) If ri > 2 is a prim number, then 

55(2 •«) - 5(l) + 5(2) + 5(«) + 5(2-«) = 2 + n + n = 2-n + 2, 

c) SS(n 2 ) = 5(1) + S(n) + S(n 2 ) = n + 2 ■ n = 3 - n < n 1 . 



A' 


5 


55 




5 


55 


n 


5 


55 


n 


~”s~ 


55 


n 


5 


55 


f 


0 


0 


11 


1 1 


11 


21 


7 


17 


31 


31 


31 


41 


41 


41 


2 


2 


2 


12 


4 


16 


22 


11 


24 


32 


8 


24 


42 


7 


36 


j 


j 


j 


13 


13 


13 


23 


23 


23 




11 


25 


43 


43 


43 


4 


4 


6 


14 


7 


16 


24 


4 


24 


34 


17 


36 


44 


11 


39 


5 


5 


5 


15 


5 


13 


25 


10 


15 


35 


7 


19 


45 


6 


25 


6 




8 


16 


6 


16 


26 


13 


28 


36 


6 


34 


46 


23 


48 


7 


7 


7 


17 


17 


17 


27 


9 


18 


37 


37 


37 


47 


47 


47 


8 


4 


10 


18 


6 


20 


28 


7 


27 


38 


19 


40 


48 


6 


36 


9 


6 


9 


19 


19 


19 


29 


29 


29 


39 


13 


29 


49 


14 


21 


10 


5 


12 


20 


5 


21 


30 


5 


28 


40 


5 


30 


50 


10 


32 



Table 1 . The values of n, S, 55. 



106 




The inequality' SS(n) < n is proved to be true for the following particular values 
n = p k , 2 ■ p k , 3 • p k and 6 • p k . 

Lemma 1. If p>2 is a prime number and k> 1, then the inequality SS(p k ) < p k holds. 



Proof 



The following inequality holds according to inequality (4) and the definition of SS. 



k-(k + \) 



SS(p k ) = Y,S(P‘)< Y.P-1 = P 

i=\ /=! “ 

The inequality 

Y k-(k + 1) k 

Lp-‘ = p- — ^ ^ P 

/= i 

is proved to be true by analysing the following cases. 

• k= 2 => 3- p < p~ . 

• k = 3 => 6- p < p* . 

• k= 4 => 10 -p< p k . 

Inequalities (6-8) are true because p> 2. 

iz}k k - A 



(5) 

(6) 

(7) 

( 8 ) 



k ■ 4 => p k > p- p k 1 > p-2" 1 = 
of this sum are kept and it is found 



P k ^ P ■ 



-I v , j 

i=0 v ' ' 



The first and the last three terms 



f 


f£-l) 




(k-\) 




(k - lV 


1 J 


io J 


+ 2- 


vl J 


+ 2- 




V 






v2 Jj 



= p\k~ - k + 2) . The inequality' 



k~ - k + 2 > 



k-(k + 1) , k-(k + 1) 

holds because k> 4, therefore p > p- is true. 



Therefore, the inequality s(p k )< p k holds. * 

Remark 2. The inequality s(p k )< p k is still true for p= 2 and £>3 because (8) holds 
for these values. Table 1 shows that the inequality is not true for p= 2 and £=2,3. 



Lemma 2. If p>2 is a prime number and £>L then the inequality SS{2- p k ) <2 ■ p k 
holds. 
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Proof 

The definition of 55 gives the following equation 

SS(p t ) = S(2) + 'ts(p‘) + 'ts(2-p‘). 

1 = 1 1=1 

Applying the inequality 5(2 ■ p‘)< p-i and (4), we have 

k k 



SS(2 ■ p k ) <2 + ^p-i -Y jP' 1 = 2 + p-k-(k~\). (9) 

/ = i / = ! 

The inequality 

2 + p ■ k ■ (k + \) < 2 ■ p k (10) 

is proved to be true as before. 

• k= 2 =>2 + 6-p<2-p~. (11) 

• k=3 => 2 + 12-p < 2- p } . (12) 

• k=4 2 -r 20 • /7 < 2 • /? 4 . ( 13 ) 

• k=5 => 2 + 30 ■ p < 2 ■ p 5 . (14) 

• k=6 => 2 +42 -p < 2 ■ p 5 . (15) 

These above inequalities (1 1-15) are true because p> 2. 



• k>6 => p* > p- p k ] > p-2' 1 = /? • . I • The first and the last fourth terms 

/=<A 7 2 

of this sum are kept finding 




= p-(2- k 2 - 4- k + 4) > 2 + p-(k 2 + k) 

The last inequality holds because k>6, therefore 2 • p k > 2 + p - k ■ (k + 1) is true. 

The inequality 55(2 ■ p ) < 2 • p k holds because (10) has been found to be true. 

* 

Remark 3. Similarly, the inequality 55(3 ■ p k ) <3 ■ p' can be proved for all (p> 3 and 
^>1) or (p= 2 and k>3). 
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Lemma 3. If p > 3 is a prime number and £>1, then the inequality 55( 6 • p k ) < 6 • 
holds. 

Proof 

The starting point is given by the following equation (16) 

SS(6-p k ) = S(2) + S{3) + S(6) + f d S(p 1 ) + f^S(2-p l ) + j:S(3-p‘) + f^S(6-p ! ). 

i=i i=i i=i 1=1 

(16) 

The inequalities S(p‘),S(2- p‘ ),S(3 ■ p'),S{6- p')< p ■ i hold for all ;>1 because 
p>5. Therefore, the inequality' 

SS( 6 • p k ) < 8 + YjP ■ l + YjP' l + HP' 1 + YjP • / = 8 + 4 • X p - / (17) 

(=1 ;=1 /=! i =\ i =\ 

holds. The inequality 55(6 ■ p k ) < 8 + 4 • p k <6 • p k is found to be true by applying 
(5) in (17). 

* 



The following propositions give the main properties of the function 55. Let d(n) 
denote the number of divisors of n. 

Proposition 1. If a is natural numbers such that Sfa)> 4, then the inequality 
S{a)>2- d(a) holds. 

Proof 

The proof is made directly as follows: 

5(a) = £ . S(d) = £ S(d) + 5(a) > X 2 + S ( a ) = 2 • (J(a) - 2) + 5(a) = 

1 *d.a \.rt*d\a Ln*da 

= 2 • d(a) + 5(a) - 4 > 2 • d(a) . * 

Remark 4. The inequality S(a)> 4 is verified for all the numbers a>4 and o^6. 

Proposition 2. If a, b are two natural numbers such that then the inequality 

55(a • b ) < d(a) ■ SS(b ) + d{b) ■ 55(a) holds. 

Proof 

This proof is made by using (2) and the simple remark that 
a, b > 0=>max {a, b] < a + b . 
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The set of the divisors of ab is split into three sets as follows: 



(l*d a*b = Md} = 

(l*d a=Md}u{i^d |b=Md}u{d!d 2 



a =Mdi * 1 Ab =Md 2 * 1 A ( d, , d 2 ) = 1 } . 



08) 



The following transformations hold according to (18). 

SS(a b)= £SM)- 2>(</) + 2>(<0 + 2 = 

{l.dja-b = Md) (l.dla-Md} {l.d|b=Hd} { 1 *dj j a =Md 2 } { 1 .d 2 | b=Md 2 } 

= SS(a) + SS(b)+ ^ '£max{S(d l ),S(d 2 )} < 

{!*<>! j.. B d l} {l.d 2 |b=gd 2) 

< 55(a) + 55(A) + £ YX s W + S ( d 2 )]= 

( 1 ; a = Kd 1 } { 1 #d 2 ! b = Md 2 } 

= 55(a) + 55(A) + ^ '^ j S(d l )+ ^ S(d 2 ) = 

{l^d 1 ja=Md 1 } {l#d 2 j b=Md 2 > { 1-d], ( a = Md 1 > {l*d 2 ]b=Md 2 } 

= SS(a ) + SS(b) + 55(a ) • [d{b) - 1] + 55(A) - [d{a) - l] 

Therefore, the inequality SS(a ■ b ) < d(a) ■ SS(b) + d(b) ■ 55(a ) holds. * 

Proposition 3. If a, b are two natural numbers such that S(a), S(b)>\ and (a,bj= 1, then 
the inequality SS(a ■ b) < SS{a) ■ SS(b ) holds. 

Proof 

Proposition 1-2 are applied to prove this proposition as follows: 

S(a),S(b)>4=>S(a)>2-d(a)andS(b)>2-d(b) (19) 

(a,b) = l => SS(a - b) < d(a) ■ SS(b) + d(b) ■ SS(a) . (20) 

The proof is completed if the inequality d(a) ■ 55(A ) + d(b) - 55(a) < SS(a ) - SS(b ) is 
found to be true. This is given by the following equivalence 
d(a) ■ SS(b ) + d(b) ■ 55(a ) < SS(a ) - SS(b ) <=> 

d(a) - d{b) < [55(a) - d(a)\ [55(A) - d{b)\. 

This last inequality holds according to (19). 

Therefore, the inequality 55(a ■ b) < 55(a) - SS(b) is true. * 



Theorem 1. If n is a natural number such that n * 8, 12, 20 then 

a) SS(n) - n + 2 if (3 p prime) n = 2 • p . 

b) 55(n) < n , otherwise. 
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( 21 ) 

( 22 ) 




Proof 

The proof of this theorem is made by using the induction on n. 

Equation (21) is true according to Remark l.a. Table 1 shows that Equation (22) 
holds for «<51 and n 8, 12, 20. Let n>51 be a natural number. Let us suppose that 
Equation (9) is true for all the number k that satisfies k<n and k does not have the 
form k=2p. p prime. The following cases are analysed: 

• n is prime => SS(nj=n, therefore Equation (9) holds. 

• n-2p, p> 2 prime => SSfnj =n+ 2, therefore Equation (2 1 ) holds. 

• (n = 2 k and k>3) or ( n - p k and A>1) => SS{n) < n according to Lemma 1 

• n = 2 • p k , p> 2 prime number and A>1 => SS(n ) < n according to Lemma 2. 

• n = 3 • p k , (p>3 prime number and k> 1) or (p= 2 and k>2) => SS(n ) < n 

according to Remark 3. 

• n = 6/7*, p>3 prime number and k>\ => SS(n) < n according to Lemma 3. 

• Jc k lc 

• Otherwise => Let n-p l '-p 2 Z \..-p s s be the prime number decomposition of n 
with /7j < p 2 < ... < p s . We prove that there is a decomposition of n=ab, (a,b)=\ 
such that S(a), Sfb)>4. Let us select a = p\‘ and b = p\' ■ pp ■ ... ■ p k ’_]' . It is not 

difficult to see that this decomposition satisfies the above conditions. The 
induction’s hypotheses is applied for a,b<n and the inequalities SS(a)<a and 
SS(b)<b are obtained. Finally, Proposition 3 gives 
SS(n) = SS(a ■ b) < SS(b ) ■ SS(a ) <b-a = n. 

We can conclude that the inequality' SS(n)<h-2 holds for all the natural number 12. 

A 

Remark 5. The above analysis is necessary to be sure that the decomposition of n=ab, 
fa,bj=\. Ska), Sfb)> 4 exists. 

Theorem 1 has some interesting consequences that are presented in the following. 
These establish new upper bounds for Smarandache’s function. 

Consequence 1 . If n > 1 is a natural number, then the following inequality 

S(n)<n + 4-2-d{n) . (23) 



holds. 



Ill 




Proof 

The proof of this inequality is made by using Theorem 1. 

Obviously, (23) is true for n=p or n=2p,p prime number. 

Let n * 8, 12, 20 be a natural number. 

We have the following transformations: 

n > SS(n ) = d) = S{n ) + ]T S(<7) - 

\*dn \.n*dn 

> S(n) + 2 • j^/ = \,n Id ^ \,n/\d j «|= S(n) + 2-(d(n) - 2) = S(n ) + 2-d(n) — 4 
Inequality (23) is also satisfied for n= 8, 12, 20. 

Therefore, the inequality S(n ) < n + 4 - 2 • d{ri) holds. * 

Consequence 2. If n > 1 is a natural number, then the following inequality holds 

S(n ) < n + 4 - min{/?j p is prime and p\n} ■ d(n) . (24) 

Proof 

This proof is made similarly to the proof of the previous consequence by using the 
following strong inequality S(d) > min{/?| p is prime and pin) . 4» 

3. Final Remark 

Inequalities (23 - 24) give some generalisations of the well - known inequality 
S(n) < n . More important is the fact that these inequalities reflect. When n has many 
divisors, the value of n + 4 - min{/?j p is prime and p\n} ■ d{n) is small, therefore the 
value of S(n) is small as well according to Inequality (24). In spite of fact that 
Inequalities (23 - 24) reflect this situation, we could not say that the upper bounds are 
the lowest possible. Nevertheless, they offer a better upper bound than the inequality 
S(n ) < n . 
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Smarandache’s function applied to perfect 

numbers 

Sebastian Martin Ruiz 
5 August 1998 



Smarandache’s function may be defined as follows: 

S(n) = the smallest positive integer such that S(n)! is divisible by rz. [l] 

In this article we are going to see that the value this function takes when n is 
a perfect number of the. form n = 2 k ~' L • (2* - 1) , p = 2* — 1 being a prime 
number. 

Lemma 1 Let n = 2* ■ p when p is an odd prime number and i an integer such 
that: 

°- z - + *(£> + 

Where e 2 (p!) is the exponent of 2 in the prime number descomposition of pi. 
E(x) is the greatest integer less than or equal to x. 

One has that S(n) = p. 

D emonst rat ion : 

Given that gcd(2\p) = 1 ( gcd = greatest common divisor) one has that 
S(n) = max{s(2*), S(p)} > S(p ) — p. Therefore S(n) > p. 

If we prove that p! is divisible by n then one would have the equality. 






where pi is the i-th prime of the prime number descomposition of pi It is clear 
that pi = 2 , p* = p , e Pt (p\) = 1 for which: 



P- 



2 ea(p!) . 

Z P 2 P 3 — 1 'P 



From where one can deduce that: 

— _ 2 e =oo-' . . ■■p e, _r l(p,) 

n " 3 A 

is a positive integer since eo (p!) — i > 0. 

Therefore one has that S(n) = p 
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Proposition 1 If n a perfec number of the form n — 2*” 1 * (2* — 1) with k a 
positive integer 7 2 k — 1 = p prime , one has that S(n) = p. 

Demonstration: 

For the Lemma it is sufficient to prove that k — 1 < e 2 (p!). 

If we can prove that 

( 1 ) 

we will have proof of the proposition since: 



i - 1 < 2*~ i - - = 



P 

2 



2 2 

As k-Lis an integer one has that k — 1 < E{^) < e? (p!) 

Proving (1) is the same as proving k < 2 lc ~ 1 r | at the same time, since k 
is integer, is equivalent to proving k < 2 k ~ l (2). 

In order to prove (2) we may consider the function: f(x) = 2 r ~ 1 — x x £ R, 
This function may be derived and its derivate is /'(x) = 2 x ~ l ln 2 — 1. 

/ will be increasing when 2 r “ 1 /n 2 — 1 > 0 resolving x: 

In (In 2 ) 



> 1 - 



In 2 



T5287 



In particular / will be increasing V x > 2. 

Therefore V x > 2 /(x) > /( 2) = 0 that is to say 2 X_1 - x > 0 V x > 2 
therefore 



2* 1 > & y £ > 2 integer 



and thus is proved the proposition. 
EXAMPLES: 



6 = 23 


5(6) = 3 


28 = 2 2 • 7 


5(28) = 7 


496 = 2 4 • 31 


5(496) = 31 


8128 = 2 6 ■ 127 


5(8128) = 127 
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ON THE DIVERGENCE OF THE SMARANDACHE HARMONIC SERIES 



Florian Luca 



For any positive integer n let S(n) be the minimal positive integer m such that 
n | ml. In [3], the authors showed that 

§ ^ (1) 

is divergent and attempted, with limited succes, to gain information about the 
behaviour of the partial sum 



A ^ x ) = E 

TL<X 



1 



by comparing it with both logs and logx -F log log x. 

In this note we show that none of these two functions is a suitable candidate 
for the order of magnitude of -4(x). 

Here is the result. For any S > 0 and x > 1 denote by 

M x ) = E s^n) 7, (2) 

n<x K } 



Then, 

Theorem. 

For any S > 0 , 



Bs(x) = log As{x). 



Bs(x) > log 2 • 



logx 
log log X 




log I 
log logx 



)' 



( 3 ) 



( 4 ) 



What the above theorem basically says is that for fixed S and for arbitrary 
e > 0 . there exists some constant C (depending on both 6 and e), such that 

A 5 (x) > 2 ^”^ lo s° lS 4* for x > C. (5) 



Notice that equation (5) asserts that A$(x) grows much faster than any polynomial 
in log(x). so one certainly shouldn’t try to approximate it by a linear in logx. 

The Proof. 



In [ 1 ], we showed that 



E 



1 



( 6 ) 



diverges for all 6 >0. Since the argument employed in the proof is relevant for our 
purposes, we reproduce it here. 

Let t > 1 be an integer and pi < P 2 < -** < Pt be the first t prime numbers. 
Notice that any integer n — p t m w r here m is squarefree and all the prime factors of 
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m belong to {pi, p 2 , Pt-i) will certainly satisfy S(n) = p t . Since there are at 
least 2 l ~ l such m’s (the power of the set {p i7 pt-i}), it follows that series (6) 
is bounded below' by 

2*— 1— <5 log 2 p t * ^ 

t> i ^ t> i 

The argument ends noticing that since 



, ... k ot-i ^ 

EV = E 



Pt 



lim 
t-*oo tlogt 



= 1, 



it follows that the exponent £ — 1 — <51og 2 Pt is always positive for t large enough. 
This proves the divergence of the series (6). 

For the present theorem, the only new thing is the fact that we do not work 
with the whole series (6) but only with its partial sum In particular, the 

parameter t from the above argument is precisely the maximal value of s for which 
P\P2---Ps < In order to prove our theorem, we need to come up with a good lower 
bound on t. 

We show- that for all t > 0 one has 



t > (1 - «) 



log x 

log log x 



(8) 



provided that x is enough large. Assume that this is not so. It then follows that 
there exists some e > 0 such that 



t < (1 - e) 



logs 

log log X 



( 9 ) 



for arbitrarily large values of x. Since t was the value of the maximal s such that 
Pi P'2 ---Ps < x, it follows that 

P\P 2 —Pt+l > X. (10) 

From a formula in [2], it follows easily that 



Pi < 2ilogi for i > 3. (11) 

It now follows, by taking logarithms in (10) and using (11), that 



(-1 t-r 1 

iogx < log Pi < Ci -r (t — i) log 2 + yy^og* + i°giogi) < 

1=1 1=3 

/■‘+ 2 

Cl + (t - 1) log 2 + J (log y + log log y)dy < 

Cl + (t - 1 ) log 2 + (: t + 2)(log (t + 2 ) + loglogif + 2)), (12) 

where C. = log 6. Since t can be arbitrarily large (because x is arbitrarily large), it 
follows that one can just work with 

logx < f(logf + 21oglogf). (13) 

Indeed, the amount (t-r2)(log(t + 2) + loglog(t+2)) appearing in the right hand side 
of (12) can be replaced by t(logt + loglogt) + /(f) where f(t) = 2 log t -f 2 log log f -f- 
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0(1) and then the sum of /(£) with the linear term from the right hand side of (12) 
can certainly be bounded above by t log log t for t large enough. Hence, 

logs < £(log£ + 2 log log t). (14) 



Using inequality (9) to bound the factor t appearing in (14) in terms of x and the 
obvious inequality 



t<( 1 - e) 



logx 
log log X 



< logx 



to bound the Vs appearing inside the logs in (14), one gets 



logI < (l-^i^^(log I og^ + 21 o g 1 og 1 °g I ) = 



or. after some immediate simplifications, 



log logx < 



2(1 -e) 

e 



log log logx. 



(15) 



Since e was fixed, it follows that inequality (15) cannot happen for arbitrarily large 
values of x. This proves that indeed (8) holds for any e provided that x is large 
enough. We are now done. Indeed, going back to formula (7), it follows that 

A{ (x) > 2‘- 1 - 51 o sp- 



or 



B 6 (x) = log A* (x) > log2(£ - 1 - <51ogpt) > log2(f - 1 - S log(2i log*)), (16) 



where the last inequality in (16) follows from (11). By inequalities (8) and (16), we 
get 

Bs(x) = £ log 2 - o(f) = log 2(1 - e) lQgX o{t). 

log logx 

Since e could, in fact, be chosen arbitrarily small, we get 



B s (x) 



log 2 



logx 
log logx 




logx \ 
log logx / ' 



which concludes the proof. 



Remark. 

We conjecture that the exact order of B$(x) 



logx / logx \ 
log log X V (log log x) 2 ) 



16) 
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A Generalisation of Euler’s function 
Tatiana Tabirca* Sabin Tabirca** 

•Transilvania University of Brasov, Computer Science Department 
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The aim of this article is to propose a generalisation for Euler’s function. This function 
is <p:N-*N defined as follows (VweAf)^(tt) = j^ = l,«](£,«) = lj. Perhaps, this is 

the most important function in number theory having many properties in number theory', 
combinatorics, etc. The main properties [Hardy & Wright, 1979] of this function are 
enumerated in the following: 

a, b g N)(a,b) = 1 => <p{a • b) — <p(a) ■ <p(b) - the multiplicative property (1) 

(VfleAf)^(i) = a. (3) 

da 

More properties concerning this function can be found in [Hardy & Wright, 1979], [Jones 
& Jones, 1998] or [Rosen, 1993], 



1. Euler’s Function by order k 

In the following, we shall see how this function is generalised such that the above 
properties are still kept. The way that will be used to introduce Euler’s generalised 
function is from the function’s formula to the function’s properties. 



Definition 1. Euler’s function by order k e N is <p k : N — » N defined by 



(Vo - P7' • PT- • - • P? Wa) = •( 1 - V . ] ■ fl - J 



Pi J 



Pi 



1 ~/p:} 



Remarks 1. 

1 . Let us assume that <p k (1) = 1 . 
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2. Euler function by order 1 is Euler’s function. Obviously, Euler’s function by order 0 
is the constant function 1. 

In the following, the main properties of Euler’s function by order k are proposed. 
Theorem 1. Euler’s function by order k is multiplicative 

(V a, b e N ) (a, b) = 1 => cp k ( a-b) = (p k ( a ) ■ cp k (b) . (4) 

Proof 

This proof is obvious from the definition. * 



Theorem 2. (Vae N) y (p k ( d ) = a k 



( 5 ) 



da 



Proof 



The function <p k {a ) = ^<p k (cl) is multiplicative because cp k is a multiplicative function. 



d\a 



If a - p m , then the following transformation proves (5) 



d p m 



f 



;=0 



/=! 



1 



A 



K~P k ) 



= 1 + ]T (p ki - p Hi ~ X) )= 1 + p km - 1 = a k 

7 = 1 

If a - p™' ■ pp • ... - pp then the multiplicative property is applied as follows: 

VM - Mp? ) • V k {p? MpP ) - P k r ■ pT ■ - Ps m ■ =a k . * 



Definition 2. A natural number n is said to be £-power free if there is not a prim number 
p such that p k | n . 

Remarks 2. 

1 There is not a 0-power free number. 

2. Assume that 1 is the only 1 -power free number. 
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The combinatorial property of Euler’s function by order k is given by the following 
theorem. This property is introduced by using the £-power free notion. 



Theorem 3. (v n e N^)(p k {n) = jjr = \,n k [(/,«*) is£ - power free] (6) 

Proof 

This proof is made using the Inclusion-Exclusion theorem. 

Let a = p™' ■ pp •.. • p™’ be the prime number decomposition of a. 

If dn, then the set S d = | = 1 ,a k \ d k |/j contains all the numbers that have the divisor 
d k . This set satisfies the following properties: 



S d =<d ,2- d 



i^y'i <Ji = 5 



J7 rj i r J2 



I id a ) is k - power ffee}= { = \,a k }- (S Pi n S Pz n ... n^J. 
The Inclusion-Exclusion theorem and (7-9) give the following transformations: 
= | (/,a*) is £- power freely* -j(S A nS f; n...nS p J| = 



= a k - S' S + Y S„ nS„ -... + (-1 ) s+1 Y S n nS„ n...nS, | = 

^ Pj ^ Pj\ Pjl V ; ^ Pjx Pjl Pjs 



1 —Jl <Jz 







Therefore, the equation (6) holds. 



* 
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3. Conclusion 

Euler’s function by order k represents a successful way to generalise Euler’s function. 
Firstly, because the main properties of Euler’s function (1-3) have been extended for 
Euler’s function by order k. Secondly and more important, because a combinatorial 
property has been found for this generalised function. Obviously, many other properties 
can be deduced for Euler’s function by order k. Unfortunately, a similar property with 
Euler’s theorem a 9(n) = 1 mod n has not been found so far. 
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A result obtained using Smarandache 
Function 

Sebastian Martin Ruiz 
21 November 1998 



Smarandache Function is defined as followed: 

S(m) = The smallest positive integer so that 5(m)! is divisible by m. [1] 

Let's see the value which such function takes for m = p?" with n integer, n > 2 
and p prime number. To do so a Lemma is required. 

Lemma IV m, n e N m, n > 2 



m n = E 



n^-rl 



\+E 



m n + ml 



•+E 



m n+1 — m n -I- m 
77r B P°s»*( m " +: ] 



ITAere iT(r) gives the greatest integer less than or equal to x. 

D emons t rat ion : 

Let's see in the first place the value taken by £pog m (m" +1 — m" + m)]. 

If n > 2: m m+1 — m n + m < m" +1 and therefore log m (m n+1 — m n + m) < 
l°g m m " +1 = n 4- 1. As a result £'[log Tn (m ,l+1 — m" 4- m)] <n|I. 

And if m > 2: mm" > 2m" => m" +1 > 2m" => m" +1 4- m > 2m" => 
m" +1 - m" + m > m" => log m (rn" +1 - m" + m) > log m m" = n => 
•^P°gm( mn+1 - m" + m)] > n 

As a result: n < £'[log m (m" +1 - m" + m)j < n + 1 therefore: 



•£’P°gm( mn+1 ~ + m )] = n if n.m > 2 



Now let’s see the value which it takes for 1 < k < n: E I 

L m * 

1 



m""*" 1 — m" -f m 



= E 



If k = 1: E I 



If 1 < * < n: £ far” 



m" +1 -* _ 



1 



= m"-*- 1-4 - m" 
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Let’s see what is the value of the sum: 



k = 1 
4 = 2 
k = 3 



n 1 



m —m 



-hi 



m n 1 —m n ~ 2 



m n 2 — m n “ 3 



/: = n — 1 
k = n 



m~ —m 

m -1 



Therefore: 



X> 

jfe=l 



f m ri+1 — m n -f m 



— m n m, n > 2 



Proposition 1 V p prime number Vn>2; 

^‘)=p n+ 1 -p n +p 

Demonstration: 

Having e p (fc) = exponent of the prime number p in the prime numbers 
descomposition of k. 

We get: 

e p (k\) = E& + £(4) + £(4) + • • ' + S(- * .. ) 

P p 3 ' p £ (l°S P *) ' 



And using the Lemma we have: 



e P l(p n+1 -p n +p)'} = E 



P n+1 -P n + P 



+E 



P n+1 - p n + p 



+■ ■ -+E 



p n+1 - p n + p 



m E l l °s ,( p “ +1 - p “+ p )] 



Therefore: 

(p" +1 ~ P n +p)! 



P pn 



£ N and 



(p r>+1 ~ P n + p — 1)! 



P"” 



£ N 



And: 



S(p p “)=p" +1 -p"+i 
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On an inequality for the Smarandache function 



J. Sandor 

Babe§-Bolyai University, 3400 Cluj-Napoca, Romania 



1. In paper [2] the author proved among others the inequality S(ab) < aS(b) for all 
a, b positive integers. This was refined to 

S(ab) < S(a) + S(b) (1) 

in [1]. Our aim is to show that certain results from our recent paper [3] can be obtained 
in a simpler way from a generalization of relation (1). On the other hand, by the method 
of Le [l] we can deduce similar, more complicated inequalities of type (1). 

2. By mathematical induction we have from (1) immediately: 

S(aiCi 2 . . . a n ) < 5(^1 ) + S(a 2 ) + . . . + S(a n ) (2) 

for all integers a ; > 1 (i = 1, . . . , n). When a\ — . . . = a n = n we obtain 

S(a n ) < nS(a). (3) 

For three applications of this inequality, remark that 

5((m!) n ) < nS(m\) ~ nm (4) 

since 5(m!) = m. This is inequality 3) part 1. from [3]. By the same way, 5((n!)^ n_1 ^ ! ) < 
(n — 1)15(72!) = (n — l)!n = n!, i.e. 

5((n!) (n - 1)! ) < nl (5) 
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Inequality (5) has been obtained in [3] by other arguments (see 4) part 1.). 

Finally, by 5( 77 2 ) < 25(n) < n for n even (see [3], inequality 1), n > 4, we have 
obtained a refinement of 5(n 2 ) < n: 

S{n 2 ) < 25(n) < n (6) 

for n > 4, even. 

3. Let m be a divisor of n. i.e. n = km. Then (1) gives S{n ) = S{km.) < S(m) + S(k), 
so we obtain: 

If 77? 1 77, then 

S{n) - 5 ( 777 ) < 5 (£-j . (7) 

As an application of (7), let d(n) be the number of divisors of n. Since JJ it = n^"^ 2 , 

&jn 

and k = nl (see [3]), and by n*in k, from (7) we can deduce that 

k<n k\n k<n 

S(77 rf( ">/ 2 ) + S{nl/n d(n)/2 ) > n. (8) 

This improves our relation (10) from [3]. 

4. Let 5(a) = u. S(b ) = v. Then 6ju! and zf!|or(ar — 1 ) . . . (x — 77 + 1) for all integers x > u. 
But from a|u! we have a|.r(x — 1 ) . . . (x — u + 1) for all x > u. Let x = u + v + k (k > 1). 
Then, clearly ab(v + 1 ) . . . (v + k)\(u + v + k)\, so we have 5[a6(u + l) . . . (u + fc)] < u + v + k. 
Here v = 5(6), so we have obtained that 

5[a6(5(6) + 1) . . . (5(6) + k)] < 5(a) + 5(6) + k. (9) 

For example, for k = 1 one has 

5[a6(5(6) + 1)1 < 5(a) + 5(6) + 1. (10) 

This is not a consequence of (2) for 77 = 3, since 5[5(6) + 1] may be much larger than 1. 
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ON A SERIES INVOLVING 5(1) - 5(2)... - 5(n) 



Florian Luca 



For any positive integer n let 5(n) be the minimal positive integer m such that 
n | m!. It is known that for any a > 0, the series 



E 



n Q 

5(1) • 5(2) ■ ... • 5(n) 



(1) 



is convergent, although we do not know who was the first to prove the above state- 
ment (for example, the authors of [4] credit the paper [1] appeared in 1997, while 
the result appears also as Proposition 1.6.12 in [2] which was wTitten in 1996). 

In this paper we show T that, in fact: 

Theorem. 

The series 

'7* n 

V (2) 



converges absolutely for every x. 

Proof 

Write 

i*r 

S(l)-S(2)-...-S(n)- 



Then 



Qn-rl _ \x\ 

a n S(n + 1) 



(3) 

(4) 



But for |x| fixed, the ratio |x|/5(n 4- 1) tends to zero. Indeed, to see this, choose 
any positive real number m, and let n m = |_m|x| 1 j ! . When n > n m , it follows 

that 5(n + l) > |_m]x| + lj > m|x|, or 5(n H- l)/|x| > m. Since m w r as arbitrary, it 
follow’s that the sequence 5(n + l)/|x| tends to infinity. 

Remarks. 

1. The convergence of (2) is certainly better than the convergence of (1). 
Indeed, if one fixes any x > 1 and any a, then certainly x n > n a for n large enough. 

2. The convergence of (2) combined with the root test imply that 



(S(l)-5(2)-...-5(n)) 1 /" 



diverges to infinity. This is equivalent to the fact that the average function of the 
logs of S, namely 

LS{x) = - log 5(n) for x > 1 

n< x 

tends to infinity with x. It would be of interest to study the order of magnitude of 
the function LS(x). We conjecture that 

LS(x) = logx — loglogx -f 0(1). (5) 
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The fact that LS{x) cannot be larger than what shows up in the right side of (5) 
follows from a result from [3]. Indeed, in [3], we showed that 

1 x X 

.4(x) = - > S(n) < 2- for x > 64. (6) 

X ^ logx “ 

Now the fact that LS(x) - logx -b log logx is bounded above follows from (6) and 
from Jensen's inequality for the log function (or the logarithmic form of the AGM 
inequality). It seems to be considerably harder to prove that LS(x) — log x+log log x 
is bounded below. 

3. As a fun application we mention that for every integer k > 1, the series 

\ x n - 

) ‘ 5(1) • 5(2) - ... • S{n) (,) 

is absolutely convergent. Indeed, it is a straightforward computation to verify that 
if one denotes by C(x) the sum of the series (2), then the series (7) is precisely 

x* d k C 

kl dx k ' ^ 

When k = x = 1 series (7) becomes precisely series (1) for a = 1. 

4. It could be of interest to study the rationality of (2) for integer values of 
x. Indeed, if the function S is replaced with the identity in formula (2), then one 
obtains the more familiar e x whose value is irrational (in fact, transcendental) at 
all integer values of x. Is that still true for series (2)? 
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A Congruence with Smarandache’s Function 

Sebastian Martin Ruiz 
25 September 1998 



Smarandache’s function is defined thus: 

S(n ) = is the smallest integer such that S(n)! is divisible by n. [1] 

In this article we are going to look at the value that has S(2 k — 1) ( mod k ) for 

all k integer, 2 < k < 97. 

One can observe in the following table that gives the continuation 
S( 2* - 1) = 1 ( mod k ) in the majority of cases, there are only 4 exeptions for 

2 < k < 97. 



k 


S( 2 “ - 1 ) 


5(2* - 1) ( mod k) 


2 


3 


1 


3 


7 


1 


4 


5 


1 


5 


31 


1 


6 


7 


1 


7 


127 


1 f 


S 


17 


1 


9 


73 


1 


10 


31 


1 


11 


89 


1 


12 


13 


1 


13 


8191 


1 


14 


127 


1 


15 


151 


1 


16 


257 


1 


17 


131071 


1 


18 


73 


1 


19 


524287 


1 
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m 




5(2* - 1) ( mod k ) 


k 


5(2* - 1) 




20 


41 


1 


59 


3203431780337 


l 


21 


337 


1 


60 


1321 


l 


22 


683 


1 


61 


2305843009213693951 


1 


23 


178481 


1 


62 


2147483647 


1 


24 


241 


1 


63 


649657 


l 


25 


1801 


1 


64 


6700417 


l 


26 


8191 


1 


65 


145295143558111 


l 


27 


262657 


1 


66 


599479 


l 


28 


127 


15 


67 


761838257287 


l 


29 


2089 


1 


68 


131071 


35 


30 


331 


1 


69 


10052678938039 


1 


31 


2147483647 


1 


70 


122921 


1 


32 


65537 


1 


71 


212885833 


1 


33 


599479 


1 


72 


38737 


1 


34 


131071 


1 


73 


9361973132609 


1 


35 


122921 


1 


74 


616318177 


1 


36 


109 


1 


75 


10567201 


1 


37 


616318177 


1 


76 


525313 


1 


38 


524287 


1 


77 


581283643249112959 


1 


39 


121369 


1 


78 


22366891 


1 


40 


61681 


1 


79 


1113491139767 


1 


41 


164511353 


1 


80 


4278255361 


1 


42 


5419 


1 


81 


97685839 


1 


43 


2099863 


1 


82 


8831418697 


1 


44 


2113 


1 


83 


57912614113275649087721 


1 


45 


23311 


1 


84 


14449 


1 


46 


2796203 


1 


85 


9520972806333758431 


1 


47 


13264529 


1 


86 


2932031007403 


1 


48 


673 


1 


87 


9857737155463 


1 


49 


4432676798593 


1 


88 


2931542417 


1 


50 


4051 


1 


89 


618970019642690137449562111 


1 


51 


131071 


1 


90 


18837001 


1 


52 


8191 


27 


91 


23140471537 


1 


53 


20394401 


1 


92 


2796203 


47 


54 


262657 


1 


93 


658812288653553079 


1 


55 


201961 


1 


94 


165768537521 


1 


56 


15790321 


1 


95 


30327152671 


1 


57 


1212847 


1 


96 


22253377 


1 


58 


3033169 


1 


97 


13842607235828485645766393 


1 1 



One can see from the table that there are only 4 exeptions for 2 < k < 97. 








We can see in detail the 4 exeptions in a table: 



L = 28 = 2 2 • 7 


5( 2 2A - 1) = 15 (mod 28) 


k = 52 = 2 2 • 13 


5(2 52 - 1) = 27 ( mod 52) 


it = 68 = 2 2 • 17 


5(2 68 - 1) = 35 ( mod 68) 


k = 92 = 2 2 • 23 


5(2 92 - 1) = 47 ( mod 92) 



One can observe in these 4 cases that k = 2 2 • p with p prime and moreover 
5(2* - 1) = |+ 1 ( modk ) 

Unsolved Question: 

One can obtain a general formula that gives us, in function of k the value 
5(2* — 1) [mod k) for ail pisitive integer values of k ?. 
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An Integer as a Sum of Consecutive Integers 

Henry Ibstedt 



Abstract: This is a simple study of expressions of positive integers as sums of 
consecutive integers. In the first part proof is given for the fact that N can be 
expressed in exactly d(L>l ways as a sum of consecutive integers, L is the 
largest odd factor of N and d(L) is the number of divisors of L. In the second 
part answer is given to the question: Which is the smallest integer that can be 
expressed as a sum of consecutive integers in n ways. 

Introduction: There is a remarkable similarity between the four definitions given below. 
The first is the well known Smarandache Function The second function was defined by K. 
Kashihara and was elaborated on in his book Comments and Topics on Smarandache 
Notions and Problems 1 . This function and the Smarandache Ceil Function were also 
examined in the author’s book Surfing on the Ocean of Numbers 2 . These three functions 
have in common that they aim to answer the question which is the smallest positive 
integer N which possesses a certain property pertaining to a given integer n. It is 
possible to pose a large number of questions of this nature. 

1. The Smarandache Function S(n): 

S(n)=N where N is the smallest positive integer which divides n!. 

2. The Pseudo-Smarandache Function Z(n): 

Z(n)=N where N is the smallest positive integer such that 1+2+...+N is 
divisible by n. 

3. The Smarandache Ceil Function of order k, S k (n): 

Sk(n)=N where N is he smallest positive integer for which n divides N* 

4 . The n-way consecutive integer representation R(n): 

R(n)=N where N is the smallest positive integer which can be represented as a 
sum of consecutive integer is n ways. 

There may be many positive integers which can be represented as a sum of positive integers 
in n distinct ways - but which is the smallest of them? This article gives the answer to this 
question. In the study of R(n) it is found that the arithmetic function d(n), the number of 
divisors of n, plays an important role. 



1 Erhus University Press, 1996 

2 Erhus University Press, 1997 
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Question 1: In how many ways n can a given positive integer N be expressed as the sum of 
consecutive positive integers? 

Let the first term in a sequence of consecutive integers be Q and the number terms in the 
sequence be M. We have N=Q+(Q+1)+ ... +(Q+M-1) where M>1. 



or 



M(2Q + M-1) 
2 




M-l 

2 



( 1 ) 

( 2 ) 



For a given positive integer N the number of sequences n is equal to the number of positive 
integer solutions to (2) in respect of Q. Let us write N=L*2 S and M=m*2 k where L and m are 
odd integers. Furthermore we express L as a product of any of its factors L=mi m 2 . We will 
now consider the following cases: 



1. 


s=0, k=0 


2. 


s=0, krO 


3. 


s*0, k=0 


4. 


s?K), k^O 


s=0. 


k=0. 



Equation (2) takes the form 



Q = 



m 1 m 2 m-l 



m 2 

Obviously we must have m*l and m^L (=N). 



(3) 



For m=mi we have Q>0 when m 2 -(mi-l)/2>0 or mi<2m 2 +l. Since mi and m 2 are odd the 
latter inequality is equivalent to m, <2m 2 or, since m 2 =N/mi, m 1 < V2N . 

We conclude that a factor m (*1 and *N) of N (odd) for which m < V2N gives a solution 
for Q when M=m is inserted in equation (2). 

Case 2. s=Q, k^O. 

Since N is odd we see form (2) that we must have k=l. With M=2m equation (2) takes the 
form 



m 1 m 2 2m- 1 
2m 2 



(4) 
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For m=l (M=2) we find Q=(N-1)2 which corresponds to the obvious solution 
N-l N + l. V7 
2 2 

Since we can have no solution for m=N we now consider m=m 2 (*1, *N). We find Q=(m r 
2m 2 +l)/2. Q>0 when mi>2m 2 -l or, since mi and m 2 are odd, mi>2m 2 Since miin 2 =N, 

m 2 =N/mi we find m > V2N . 

We conclude that a factor m (*1 and *N) of N (odd) for which m > V2N gives a solution 
for Q when M=2m is inserted in equation (2). 

The number of divisors of N is known as the function d(N). Since all factors of N except 1 
and N provide solutions to (2) while M=2, which is not a factor of N, also provides a 
solution (2) we find that the number of solutions n to (2) when N is odd is 



n=d(N)-l 



(5) 



Case 3. s^O, k=0. 

Equation (2) takes the form 

Q = i ( J& 2 s+1 -m+l) (6) 

2 m 

Q>1 requires m 2 <L*2^V We distinguish three cases: 

Case 3.1. k=0, m= 1 . There is no solution. 

Case 3.2. k=0, m=m! . Q>1 for m\<m 2 2* r} with a solution for Q when M=mi . 

Case 3.3. k=0, m=m 1 m 2 . Q>1 for L<2^ ! with a solution for Q when M=L. 

Case 4, s^O. fc*Q. 

Equation (2) takes the form 

Q = A(J?^2 s ~ k+1 -m-2 k +l) (7) 

2 m 

Q is an integer if and only if m divides L and 2 S ~ M =1. The latter gives k=s+l. Q>1 gives 

Q = I(£hEi + i)_ m .2 s s 1 (8) 

2 m 

Again we distinguish three cases: 
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Case 4.1. 


k=s+l, m=l. 


Q>1 for L>2**' with a solution for Q when M=: 


Case 4.2. 


k=s+l, m=m 2 


Q^l for mi>m 2 2^ 1 with a solution for Q when 






M=m 2 2~ l 


Case 4.3. 


k=s+l, m=L 


Q>1 for l-L-2^1. No solution 



Since all factors of L except 1 provide solutions to (2) we find that the number of solutions 
n to (2) when N is even is 

n=d(L)-l ( 9 ) 



Conclusions: 



• The number of sequences of consecutive positive integers by which a positive integer 
N=L-2 5 where L=1 (mod 2) can be represented is n=d(L)-l. 

• We see that the number of integer sequences is the same for N=2’L and N=L no matter 
how large we make s. 

• When L<2 5 the values of M which produce integer values of Q are odd, i.e. N can in 
this case only be represented by sequences of consecutive integers with an odd number 
of terms. 

• There are solutions for all positive integers L except for L=l, which means that N=2 5 
are the only positive integers which cannot be expressed as the sum of consecutive 
integers. 

• N=P-2 5 has only one representation which has a different number of terms (<p) for 
different s until 2^>P when the number of terms will be p and remain constant for all 
larger s. 

A few examples are given in table 1. 

Table 1 . The number of sequences for L=105 is 7 and is independenf of s in N=L-2 S . 



N=105 


s=0 


N=210 


01 

It 


N=3360 


s=5 

L>2 S+1 


N=6720 


s=6 

L<2 s * 1 


Q 


M 


Q 


M 


Q 


M 


Q 


M 


34 


3 


69 


3 


1119 


3 


2239 


3 


19 


5 


40 


5 


670 


5 


1342 


5 


12 


7 


27 


7 


477 


7 


957 


7 


1 


14 


7 


15 


217 


15 


441 


15 


6 


10 


1 


20 


150 


21 


310 


21 


15 


6 


12 


12 


79 


35 


175 


35 


52 


2 


51 


4 


21 


64 


12 


105 
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Question 2: Which is the smallest positive integer N which can be represented as a sum of 
consecutive positive integers in n different ways. 

We can now construct the smallest positive integer R(n)=N which can be represented in n 
ways as the sum of consecutive integers. As we have already seen this smallest integer is 
necessarily odd and satisfies n=d(N)-l. 

With the representation N = 2 -pJ j we have 

d(N)=(a 1 +l)(a 2 +l)...(a j +l) 

or 

n+l=(a 1 +l)(a 2 +l)...(cxj + l) ( i0 ) 

The first step is therefore to factorize n +1 and arrange the factors (oq+ 1 ), (a 2 +l) ... ( 04 + 1) 
in descending order. Let us first assume that ai>a 2 > ... >04 then, remembering that N must 
be odd, the smallest N with the largest number of divisors is 

R(n)= N = 3 ttl 5 az l ai ... p 

where the primes are assigned to the exponents in ascending order starting with pi=3. Even' 
factor in ( 10 ) corresponds to a different prime even if there are factors which are equal. 

Example: n =269 

n+l= 2*3 3 *5 = 5 *3 *3 -3 -2 
R(n) = 3 4 -5 2 *7 2 *ll 2 * 13= 156080925 

When n is even it is seen from (10) that oq, a 2 , ... cxj must all be even. In other words the 
smallest positive integer which can be represented as a sum of consecutive integers in a 
given number of ways must be a square. It is therefore not surprising that even values of n 
in general generate larger smallest N than odd values of n. For example, the smallest 
integer that can be represented as a sum of integers in 100 ways is N=3 100 , which is a 48- 
digit integer, while the smallest integer that can be expressed as a sum of integer in 99 ways 
is only a 7-digit integer, namely 3898125. 



Conclusions: 



• 3 is always a factor of the smallest integer that can be represented as a sum of 
consecutive integers in n ways. 

• The smallest positive integer which can be represented as a sum of consecutive integers 
in given even number of ways must be a square. 
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Table 2. The smallest integer R(n) which can be represented in n ways as a sum of 

consecutive positive integers. 



n 


R (n) 


R(n) in factor form 


1 


3 


3 


2 


9 


3 2 


3 


15 


3-5 


4 


81 


3 4 


5 


45 


3 2 -5 


6 


729 


3 s 


7 


105 


3-5-7 


8 


225 


3 2 5 2 


9 


405 


3 4 5 


10 


59049 


3 10 


11 


315 


3 2 5-7 


12 


531441 


3 12 
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A NEW INEQUALITY FOR THE SMARANDACHE FUNCTION 



Mihaly Bencze 
Str. Harmanului 6 
2212 Sacele 3, jud. Brasov 
Romania. 



Theorem. Let S(m) — min{ks N: ml k!} be the Smarandache Function, 
and a k , b k 6 N* (k=l,2,...,n), then we have the following inequalitv 

n b k n 

S(IT ( a k !) )<Ya k b k 



Proof: 






Z( a k bQ! 

fc-1 


i(a k b,J! 

fc-1 

* 


Il( a k bt)! 

k-1 


b k 

ri( a t -) 

k-1 


IT( a t bk)! 

k-1 


. b k 

]I(a k !) 

k-1 


~ a 2 b 2 + ... + a m 


i b m w a 2 b 2 + 

l f 


••• + a n b n , 


i 


/ ' a 2 b 2 




a k b k) ( 3a k \ 

a k /.Ja J 


/ 2a k U 
\ a t // e N* 




From this result 


b k 

s cn ( a k o 


> £ i >ii>i 






k-1 
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A FORMULA OF THE SMARANDACHE FUNCTION 



Maohua Le 

Department of Mathematics, Zhanjiang Normal College 
Zhanjiang, Guangdong, P.R. China. 



Abstract. In this paper we give a formula expressing the 
Smarandache function S(n) by means of n without using the 
factorization of n. 



For any positive integer n, let S (n) denote the Smarandache 
function of n. Then we have 

(1) S(n) = mini a I ae N, n I a!}, 

(See [1]). In this paper we give a formula of S(n) without 
using the factorization of n as follows: 

Theorem. For any positive integer n, we have 

f „ "( n sin (k! n / n)) 2 ] 



(1) S(n) = n+1 - IE n 


J 


lc=l 




Proof. Let a = S(n). 


It is an obvious fact that 1 ^ a ^ n 


We see from (1) that 




(2) n 1 k!, k = a, a+1, .. 


- , n. 


It implies that 
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- (n sin (k! 7t / n)) 2 0 

(4) n = n = 1, k = a, a+1, , n. 

On the other hand, since n / k! for k ~ 1, ..., a-1, we 
have sin (klit/n) * 0 and 

k! n tt 

(5) (nsin ) 2 £ (n sin — ) 2 > 1, k=l, ... ,a-l. 

n n 



Hense, by (5), we get 

-(n sin (k! re/ n)) 2 

(6) 0< n <l/n, k=l, ... , a-1. 



Therefore, by (4) and (6), we obtain 



n -(n sin (k! Tt/n)) 2 

(7) n+l-a < E n < n +1 - a + (a-1) / n < n +2 - a. 

k=l 



Thus, by (7), we get (1) immediately. The theorem is proved. 



Reference 

1. F Smarandache, A function in the number theory, Smarandache 
function J. 1 (1990), No. 1, 3 - 17. 
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ON THE DIOPHANTINE EQUATION S(n) = n 



Maohua Le 

Department of Mathematics* Zhanjiang Normal College 
Zhanjiang, Guangdong, P.R. China. 



Abstract. Let S(n) denote the Smarandache function of n. 

In this paper we prove that S(n) = n if and only if 
n = 1, 4 or p, where p is a prime. 

Let N be the set of all positive integers. For any 
positive integer n, let S(n) denote the Smarandache function of 
n (see[l]). It is an obvious fact that S(n) ^ n. In 
this paper we consider the diophantine equation 

(1) S(n) = n, ne N. 

We prove a general result as follows: 

Theorem. The equation (1) has only the solutions n = 1,4 
or p, where p is a prime. 

Proof. If n = 1,4 or p, then (1) holds. Let n be 
an another solution of (1). Then n must be a composite integer 
with n > 4. Since n is a composite integer, we have n = uv, 
where u,v are integers satisfying us vs 2. If u* v, 
then we get n I u!. It implies that S(n) <u = n/ v<n, a contradiction. 
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If u = v, then we have n = u 2 and n| (2u)! 

It implies that S(n) £ 2u. Since n > 4, we get u > 2 and 
S(n) s 2u < u 2 = n, a contradiction. Thus, (1) has only the 
solution n = 1,4 or p. The theorem is proved. 



Reference 



1. F Smarandache, A function in the number theory, Smarandache 
function J. 1 (1990), No.l, 3-17. 
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ON SMARANDACHE DIVISOR PRODUCTS 



Maohua Le 

Department of Mathematics, Zhanjiang Normal College 
Zhanjiang, Guangdong, PR. China. 



Abstract. In this paper we give a formula for Smarandache 
divisor products. 

Let n be a positive integer. In [1, Notion 20], the 
product of all positive divisors of n is called the Smarandache 
divisor product of n and denoted by P d (n). In this paper 
we give a formula of P d (n) as follows: 

r k 

be the factorization of n, 

if n is not a square, 

1), if n is a square. 



r i 





Theorem. Let 


n - 


Pi ... 


Pk 


and 


let 










f V4( r,+ 


i)~ 


• o-k+ 


i). 


0) 


r (n )= > 










l ^((r, 4 


-i) 


- ( r k+ 


i)- 



Then we have P d (n) = n r(n) . 

Proof. Let f (n) denote the number of distinct positive 
divisors of n. It is a well known fact that 

(2) f ( n ) = (r, + 1) ... (r k + 1), 



144 




(See [2, Theorem 273]). If n is not a square and d is 
a positive divisor of n, then n/d is also a positive divisor 
of n with n / d * d. It implies that 

(3) Pj (n) - n f(n)/2 . 

Hence, by (1), (2) and (3), we get P d (n) = n r(n) . 

If n is a square and d is a positive divisor of n with 
d * \/n, then n/d is also a positive divisor of n with n / d * d. 

So we have 

n f(n)/2 

(4) P d (n) = = n (f( "- 1)/2 . 

•Jn 

Therefore, by (1), (2) and (4), we get P d (n) = n r(n) too. The 
theorem is proved. 



References. 

1 . Dumitrescu and V. Seleacu, Some Notions and Questions 
In Number Theory, Erhus Univ. Press, Glendale, 1994. 

2. G.H.Hardy and e.M. Wright, An Introduction to the 
Theory of numbers, Oxford Univ. Press, Oxford, 1938. 
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ON THE SMARANDACHE N-ARY SIEVE 



Maohua Le 

Department of Mathematics, Zhanjiang Normal College 
Zhanjiang, Guangdong, P.R. China. 



Abstract. Let n be a positive integer with n >1 . 
In this paper we prove that the remaining sequence of 
Smarandache n-ary sieve contains infinitely many composite 
numbers. 



Let n be a positive integer with n > 1 . Let S n denote 
the sequence of Smarandache n-ary sieve (see [1, Notions 
29-31]). For example: 

S 2 = {1,3, 5,9,1 1,13,17,21,25,27, ...} , 

S 3 ={ 1,2,4,5, 7,8,10,11, 14,16,1 7,1 9,20, ...} 

In [1], Dumitrescu and Seleacu conjectured that S n contains 
infinitely many composite numbers. In this paper we verify 
the above conjecture as follows: 

Theorem. For any positive integer n with n>l , 

S n contains infinitely many composite numbers. 

Proof. By the definition of Smarandache n -ary sieve 
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(see [1, Notions 29-3 1]), the sequence S„ contains the 
numbers n k +1 for any positive integer k . If k is an odd 
integer with k > 1, then we have 

(1) n k +l= (n+lXn 1 - 1 - n k ’ 2 + ... + 1). 

We see from (1) that (n+l)| (n k +l ) and n k +l is a composite 
number. Notice that there exist infinitely many odd integers 
k with k>l.Thus, S n contains infinitely many composite 
numbers n k + 1 . The theorem is proved. 



References. 

1 . Dumitrescu and V. Seleacu, Some Notions and Questions 
In Number Theory, Erhus Univ. Press, Glendale, 1994. 
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PERFECT POWERS IN THE SMARANDACHE PERMUTATION SEQUENCE 



Maohua Le 

Department of Mathematics, Zhanjiang Normal College 
Zhanjiang, Guangdong, P.R. China. 

Abstract. In this paper we prove that the Smarandache 
permutation sequience does not contain perfect powers. 

Let S = { Sn) n=1 be the Smarandache permutation sequence. 
Then we have 

(1) Sj =12, ^ = 1342, Sj = 135642, s 4 = 13578642, ... . 

In [1, Notion 6], Dumitrescu and Seleacu posed the following 
quiestion: 

Question . Is there any perfect power belonging to S? 

In this respect, Smarandache [2] conjectured: no! In this 
paper we verify the above conjecture as follows: 

Theorem. The sequence S does not contain powers. 

Proof. Let s n be a perfect power. Since 2 [ s „ by (1), 
then we have 

( 2 ) 4 | S „. 

Since s t = 12 is not a perfect power, we get n > 1. Then 
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from (1) we get 



(3) s n = 10 2 a + 42, 

where a is a positive integer. Notice that 4| 10 2 and 4 / 42. 
We find from (3) that 4 / s n , which contradicts (2). Thus, 
the theorem is proved. 



References 

1 . Dumitrescu and Seleacu, Some Notions and Questions 
In Number Theory, Erhus Univ. Press, Glendale, 1994. 

2. F.Smarandache, Only Problems, not Solutions! Xiquan 
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ON SMARANDACHE PSEUDO - POWERS OF THIRD KIND 



Maohua Le 

Department of Mathematics, Zhanjiang Normal College 
Zhanjiang, Guangdong, P.R. China. 



Abstract. Let m be a positive integer with m>l. In 
this paper we prove that there exist infinitely many m* perfect powers 
which are Smarandache pseudo - m * powers of third kind. 



Let m he a positive integer with m > 1 . For a positive 
integer a , if some nontrivial permutation of the digits is an 
m* power, then a is called a Smarandache pseudo - m* power. 

There were many questions concerning the number of Smarandache 
pseudo - m* powers (see [ 1, Notions 71, 74 and 77]). In 
general, Smarandache [2] posed the following 

Conjecture. For any positive integer m with m>l, there 
exist infinitely many m powers which are Smarandache pseudo-m th 
powers of third kind. 

In this paper we verify the above conjecture as follows. 

Theorem. For any positive integer m with m > 1, there 
exist infinitely many m th powers are Smarandache pseudo-m * 
powers of third kind. 

Proof. For any positive integer k , the positive integer 
is an m* power. Notice that 0 ...01 is a nontrivial permutation 
of the digits of 10 km and 1 is also an m* power. 

It implies that there exist infinitely many Smarandache pseudo - m * 
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powers of third kind. The theorem is proved. 
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AN IMPROVEMENT ON THE SMARANDACHE DIVISIBILITY THEOREM 



Maohua Le 

Department of Mathematics, Zhanjiang Normal College 
Zhanjiang, Guangdong, P.R. China. 



Abstract. Let a, n be positive integers. In this paper 
we prove that n I (a n - a)[n / 2]! 



For any positive integer a and n , Smarandache [3] proved that 

(1) n I (a n - a)(n - 1)!. 

The above division relation is the Smarandache divisibility 
theorem (see [1, Notions 126]). In this paper we give an 
improvement on (1) as follows: 

Theorem. For any positive integers a and n, we have 

(2) n I (a n - a)[ n / 2]!, 

where [n / 2] is the largest integer which does not exceed n / 2. 

Proof. The division relation (2) holds for n 5 9, we may 
assume that n>9. By Fermat’s theorem (see [ 2, Theorem 71]), 
if n is a prime, then we have 

(3) n I (a n - a). 
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for any a . We see from (3) that (2) is true. 

If n is a composite number, then we have n = pd, where 
p,d are integers satisfying p * q £ 2. Further, if p * q, then 
we have n|p! It implies that n|(n / q)! Since q>2, we get 

(4) n | In/21! 

If p = q , Then n = p 2 and 

(5) n | (2p)! 

Since n > 9 , we have n £ 4 2 , p £ 4 and 2p <. n / 2. Hence, 
we see from (5) that (4) is also true in this case. The 
combination of (3) and (4) , the theorem is proved. 
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ON PRIMES IN THE SMARANDACHE PIERCED CHAIN 



Maohua Le 

Department of Mathematics, Zhanjiang Normal College 
Zhanjiang, Guangdong, P.R. China. 



Abstract. Let C = (c „} n=1 be the Smarandache pierced 
chain. In this paper we prove that if n > 2, then c„ / 101 
is not a prime. 

For any positive integer n, let 
(1) c n = 101* 100010001 ...0001. 

n -I times 



Then the sequence C = {c „} n=1 is colled the Smarandache 
percied chain (see[2. Notion 19]). In [3], Smarandache 
asked the following question: 

Question. How many c n /101 are primes? 

In this paper we give a complete anser as follows: 

Theorem. If n > 2, then c „ / 101 is not a prime. 

Proof. Let C„ = e 2rJ ~ 1 / n be a primitive roof of unity with 

the degree n , and let 

f n (x) = Il ( X -C n k )- 

Is ksn 
g cd ( k, n) =1 



Then f „(x) is a polynomial with integer coefficients. Further, 

it is a well known fact that if x>2,then f„(x)>l (see 

[1]). This implies that if x is an integer with x > 2 , then f n (x) is an 

integer with f n (x) > 1 . On the other hand, we have 



(2) x n - 1 = n f d ( x ) 

d / n 

We see from (1) that if n > 1 , then 
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